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Experiments on two-dimensional flow over a 
normal wall 


By MIKIO ARIE 
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SUMMARY 

Measurements of the velocity, pressure, and turbulence behind 
a series of normal plates in the uniform test section of an air tunnel 
are described, the oscillation of the wake being prevented in all but 
ye test through use of symmetrically located tail plates. By a 
combination of experimental and computational techniques, 
details of the pattern of flow over a wall on a plane boundary in an 
infinite fluid are closely approximated. A significant difference 
is indicated between the characteristics of such a flow and those of 
the flow past an isolated plate with oscillating wake. 


INTRODUCTION 

For want of experimental information about the characteristics of flow 
past a sufficiently great variety of boundary shapes, one is forced either to 
perform specific tests upon each body design that comes into question, or 
to estimate the expected characteristics from what is known about flow past 
the individual elements of which the body consists. Under some conditions 
the latter method is quite effective; often the purely qualitative indication 
that it provides is far better than none at all, and sometimes the combination 

‘f available data with mathematical analysis permits the indication to be of 
quantitative value as well. Under other conditions the information at hand 
is merely sufficient to throw serious doubt upon even the qualitative value of 
such a procedure. 

A case in point is the utilization of experimental data for two-dimensional 
flow past a cylindrical body to predict the conditions of motion that would 
prevail if a boundary were to exist at the plane of symmetry. Assume, for 
example, that it was desired to estimate the force exerted by the wind upon 
a wall, or the effectiveness of such a wall in either shielding something behind 
it or in promoting the diffusion of some substance in its wake. Drag co- 
efficients and eddy patterns are well known for rectangular plates (Fage & 
Johansen 1927, Goldstein 1938) but it is also known (Roshko 1955) that 
interference with an oscillating wake is likely to change the primary nature 
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of the flow. Just how great this change might be in any particular case can 
be determined as yet only by experiment. 

Whereas the wake behind an immersed body is of basic importance in 
such professions as aeronautics, ballistics, and marine engineering, it is the 
wake behind a boundary irregularity, such as the wall, that is of primary 
interest in hydraulics. Whether the wall typifies roughness elements or part 
of an orifice or valve makes little difference; and the standing eddy in its 
wake also has much in common with those which form at any abrupt changes 
in conduit section. Yet so little is known about their geometrical pro- 
portions, about the dynamics of the mean flow within them, and about the 
generation of turbulence along their borders, that the experimental analysis 
of such standing eddies has received the attention of the Lowa institute of 
Hydraulic Research for nearly a decade. ‘The most recent in this series of 
investigations has been a study of the wake behind a simple wall (Arie 1955), 
the purpose of the present paper being to make generally available the most 
essential of the results. 

Flow past a plate that is fully surrounded by fluid ditfers in several ways 


from flow past one that is in contact with a continuous boundary. — First of 


all, the former may be considered to lie in the path of essentially irrotational 


1 


flow, whereas the latter will be immersed in a boundary laver of indefinite 


thickness. Secondly, as has already been noted, the flow in the wake ot an 
isolated plate is free to oscillate about the plane of symmetry, whereas such 
oscillation is wholly prevented in the vicinity of a rigid longitudinal 


boundary. Now the lateral freedom or restraint of the wake marks a 
dynamically basic distinction between the two cases, whereas the boundary 
layer is a secondary factor that is not determined by local conditions. 
Hence, in the experiments under discussion, boundary laver development 
both upstream and downstream was obviated by placing the plate at the 
axis, rather than at the floor, of an air tunnel, the lateral boundary restraint 
at the rear being effected by means of a tail plate which was only slightly 
longer than the standing eddy. 

Whether a plate fully surrounded by the fluid or one in contact with a 
boundary is simulated in a tunnel, the presence of the test section tloor and 
ceiling will introduce an additional threefold effect. First, these surfaces 
will correspond geometrically to planes of symmetry between the given 
plate and standing eddy, on the one hand, and their first images in a series 
having a finite lateral spacing, on the other; the velocity of the surrounding 
field being increased thereby. Second, there will be an additional aug- 
mentation of velocity as a result of the development of the boundary laver 
along the tunnel surfaces. Finally, despite the use of a foreshortened 
downstream boundary, the conversion of energy in the wake of the test plate 
will produce something akin to a boundary layer effect, which will cause 
the velocity of the flow outside the wake to become still higher than it would 


be in an infinite fluid. 
Three courses were thus open to the investigator. ‘The pattern of flow 
could be determined as a function of the relative spacing of a series of plates, 
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with the double aim of permitting extrapolation to the infinite case and of 
providing supplementary information about flow between opposite walls 
in a passage of finite size. Secondly, by means of various experimental 
artifices, an effort could be made to simulate the infinite case itself without 
recourse to extrapolation. Finally, the data obtained for any spacing 
could be recomputed for infinite spacing by means of an analytical approxi- 
mation. In actuality, these three courses were followed in sequence, the 
first serving as a guide to the partial refinements of the second, and the 
final correction proceeding analytically. 


EXPERIMENTAL PROCEDURI 
lhe air tunnel used for the experiments in question was of the closed 
circuit tvpe, with a uniform test section 3 feet high, 3 feet wide, and 10 feet 


long. Air speeds could be varied at will from 10 to 50 feet per second, 
but for ease in measurement all tests were made in the upper part oi 
this range. ‘The main air speed was indicated by the change in pressure 


at the bell entrance to the test section; the corresponding pressure-velocity 
relationship had been determined through test-section traverses with a 
Prandtl-type Pitot tube, which in turn had been calibrated in the irrota- 
tional core of a jet from a 6-inch rounded orifice. No effort was made t 
reduce the turbulence of the flow, beyond the use of vanes at the turns and 
a fourfold reduction in cross-sectional area between the plenum chambet 


and the test section; the root-mean-square \ elocity 


fluctuation was about 
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Figure 1. Definition sketch of plate arrangement. 


All test plates were introduced into the tunnel approximately midway 
between the ends of the test section, completely spanning its 3-foot width. 
‘They were invariably 3 inches in effective height, 1 4 inch in thickness, and 
bevelled at an angle of 60° along the edges on the rear side. All but one 
were of tempered hardboard, the one being of carefully machined brass. 
Each was guyed with piano wire, at its one-third points, to the floor and 
ceiling of the tunnel to maintain its alignment and eliminate all possibility 
of vibration. ‘The tail plates (see figure 1), which were of 1/16-inch sheet 
aluminium, were attached to the test plates at mid-height, and likewise 
guyed in place. They extended a distance of 30 inches to the rear, which 


12 
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had been found -by observation to be longer than the eddy by a sufficient 
amount for the size of the latter to be unaffected by a further lengthening of 
yf the plate. By successively introducing dummy plate assemblies at 
planes of symmetry on either side of the centrally located test plate, the 
ratio Bb of spacing to width could be varied through the consecutive 
values of 12, 6, 3, etc. 

Because of the two-dimensional nature of the flow, it was possible to 
make measurements of the magnitude and direction of the mean velocity at 
any point by means of one of two alternative Pitot cylinders. Each was 
constructed of 1 8-inch brass tubing mounted horizontally on a fork-shaped 
support with bearings 6 inches on centre. One cylinder contained a 
15-inch hole on each side of a central plug (see figure 2), the two portions 











f the tube having been rotated differentially about their common axis tll 
the holes lav at an angle of 30° to each other in the normal plane. A pro- 
tractor at either bearing permitted the common angular position of the 
holes to be read to within about 0-2. ‘The other cylinder differed from the 
first only in the angle—now 180 —between the holes. After careful 
calibration of each cylinder, the magnitude of the velocity could be com- 
puted from the differential pressure between the respective holes, the 180 
cvlinder being oriented to yield the maximum reading, and the 30° cylinder 
being turned so that first one and then the other hole lay directly upstream 
(in the latter instance the direction of the flow was first determined by noting 
the angle at which the differential reading was exactly zero). The 30 
cylinder was used so long as the velocity gradient was not excessive; since, 
in zones of high gradient, this cylinder was subject to a pronounced error 
in both direction and magnitude of the velocity, the 180° cylinder was then 
employed, and the direction determined by successive steps in plotting the 
streamlines. 

In the zones in which both the inclination and the magnitude of the 
velocity vector were directly measurable, it was a simple matter to deter- 
mine and plot, from vertical traverses with the 30° cylinder, the distribution 
of the longitudinal component of velocity over a series of normal sections. 


Graphical calculations of the areas % = [udy enclosed by these curves 
(figure 3), as a function of distance normal to the tail plate, then permitted 
successive streamlines ys = C to be located for any selected increment Ai. 
In those zones for which only the magnitude of the velocity could be 
measured, this was first assumed to be indicative of the longitudinal com- 
ponent. Evaluation of the streamline location according to the foregoing 
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method y.elded the inclination of the velocity vectors as a first approximation, 
and a second approximation was usually sufficient to determine the flow 
pattern with acceptable precision. Figure + shows the pattern obtained 
trom the single run made without tail plate for purposes of reference. 

This method of determining the flow pattern was followed for B / equal 
to 12 and to 6 without undue experimental difficulty. All attempts t 


repeat the measurements for Bb equal to 3 met with no success, however, 











Figure 3. Determination of streamline locatior 


as the pattern was then not only unsymmetrical about the axis of any plate 
but unstable as well. ‘This situation is in approximate accord with pre- 
vious experience, flow through lattice screens being known to become 
unstable as the solidity ratio approaches the magnitude 0-5. For this 
reason, no attempt was made to carry the study any farther in the direction 
of decreasing spacing, which meant that a sufficient number of cases to 


« 











Figure +. Approximate pattern of mean streamlines for flow without tail plate 


Bh 12. 


permit extrapolation to greater spacing could not thus be realized. A 
ditferent tack was hence followed in approximating the desired conditions 
of infinite spacing. 

Since the profile of the standing eddy behind the test plate was roughly 
elliptical in form (figure 5), it was assumed that the streamlines at a con- 
siderable distance would approximate those of potential flow ‘past either an 
ellipse or a Rankine oval of comparable dimensions. ‘The latter being 
somewhat easier to describe mathematically, the location and strength of 
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the source and sink were determined so as to yield a profile coinciding most 
closely with the measured eddy profile for B 6 = 12. ‘The resulting stream 
tunction had the form 


/, 0-564 1 y 1 V 
ib Vy + 0°56 bu,| tan ! —~—— — tan-! —~—— ], 
sili . v—3-45 y+3-46 


the profile of the oval corresponding to % = 0. From this relationship it 
was possible to evaluate the coordinates of any desired streamline in the 


4 


vicinity of the test-section wall; the maximum deviation of the streamline 
at the approximate lateral distance of 1-5 feet was found to be about 1-5 inches. 
Because of the boundary layer development along these walls, however, a 
still further change in their shape would be necessary to eliminate this 
additional wall effect. ‘The displacement thickness of the boundary layer 
as a function of distance along this tunnel was hence evaluated by a com- 
bination of velocity measurement with extrapolation, in accordance with 
known boundary layer relationships; its entire variation was found to be 
slightly less than 1 4 inch over the 10-foot length of test section. False 
boundaries of hardboard were then introduced along the ceiling and floor 
of the test section; they were curved according to the stream function of 
the eddy profile, flared by an amount equal to the displacement thickness 
of the boundary layer, and smoothly joined to the bell entrance of the test 


section. 
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Figure 5. Derivation of stream function for Rankine oval having proportions 


of standing eddy. 


With the falsework in position, the velocity distribution over successive 
vertical sections throughout the eddy region was again measured with the 
Pitot cylinder. In addition to such velocity traverses, measurements were 
also made of the pressure distribution at the same sections, and of the pressure 
over the faces of the test plate and tail plate. Pressure measurements at 
the plates themselves required simply the introduction of piezometer 
orifices along the one side, the velocity measurements being restricted to the 
same side in order to avoid the disturbing effect of the piezometer connec- 
tions on the other. ‘The pressure traverses in a plane normal to the two 
plates involved the use of a third plate, 1 8 inch thick, 5 inches wide (in the 
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longitudinal direction), and about 18 inches high, sharply bevelled on 
the leading edge, and having a series of piezometer orifices on one 
side in a vertical row, which could be placed at any section. Finally, 
measurements weré made with the constant-temperature type of hot-wire 
anemometer (Hubbard 1954) at the same vertical sections, the wires and 
circuit being so arranged as to yield the root-mean-square values of the 
three components of fluctuation and the cross product of the components 
in the plane normal to the two test plates. 


ANALYSIS OF MEASURED DATA 

Had the revision in tunnel profile actually been sufficient to simulate 
infinite-fluid conditions, the experimental results would have yielded 
directly the patterns of velocity, pressure, and turbulence produced by flow 
over a normal wall. However, it will be recalled that the boundary pro- 
file was determined on the assumption of irrotational flow around a Rankine 
oval, whereas flow over a wall—though it may be essentially irrotational at 
the outset—departs from this state more and more with distance down- 
stream, as a result of the pronounced shear and generation of turbulence 
at the border of the standing eddy. In other words, although a fairly suc- 
cessful effort had been made to compensate for the presence of the eddy and 
the development of the boundary layer on the tunnel walls, there remained 
uncompensated the effect of the wake of the plate itself, the importance of 
which was not at once appreciated. Finally, rather than modify the tunnel 
profile and repeat the measurements a second time, it was decided to 
determine by analytical means the magnitude of the corrections still 
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Figure 6. Displacement etfect of standing eddy and turbulent wake 


With reference to figure 6 (in which the wall height A is introduced as 
reference length), the total effect of the wall upon the surrounding flow 
could be considered equivalent to that of the combination of the Rankine 
body (shown by double cross-hatching) with an additional tail section 
(shown by single cross-hatching) corresponding to the displacement thick- 
ness of the flow in the wake. ‘The total thickness of the combination was 


assumed to be given to a first approximation by 
{ u—t ; 


nar 


} —<—(}/, 
all Umax 


the quantity 0 representing the value of y, just outside the wake, at which 
the velocity becomes a maximum. 
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Just as the profile of the eddy was simulated by the combination of a 
source and a sink with a uniform flow, the composite profile of the eddy 
and the tail section (or the profile of either one alone) could be realized by 
a distribution of sources and sinks, or, more conveniently, by a distribution 
of doublets (a method analogous to that of Landweber 1951), in place of the 
single source and sink of the Rankine body. For the case of an infinite 
fluid the stream function would then have the form 

m(t) y ; 
y= —YUy a (xt? +2 at, 
in which m(t) is the doublet strength per unit length as a function of distance 
t along the line of symmetry and a is the abscissa of the initial point of the 
distribution which will yield the desired profile for the condition “ = 0 
For the case of flow that is confined between parallel walls (i.e., for a series 
of similar bodies) with the spacing #, on the other hand, the stream function 
would have to be written as 
7 m(t)sin 2ry/B 


— | it, 
BJ, cosh2z(x—t) B—cos2zv B 





y Vulg + 


the doublets being distributed symmetrically along a series of parallel axes 
midway between the walls. 

Because the tunnel boundaries in the present experiments were formed 
according to the approximate stream function of the eddy alone, the 
differences between the velocity distributions indicated by the proper 
space derivatives of these two equations must be due almost entirely to 
the tail section or wake. Evaluation of the difference, in the form of an 
approximate correction factor, hence depended upon determination of 
the doublet distribution m(t) for the tail section alone. Fortunately, the 
elements shown in figure 6 were sufficiently slender for a considerable 
simplification (again by analogy with the method of Landweber 
(1951)) to be made. Thus, upon replacing m(t) by m(x) with ys = 0 in the 
stream function for the infinite fluid, it could be shown that, with accep- 


F | | ¥—@ 
mx) -Uyg V(X)| 57+ tan -—— 
\ 


‘The computed distribution of m(x) tor the tail section (and as well for the 


table accuracy, 


composite profile, with a value of a that would yield a satisfactory approxi- 
mation in the vicinity of the plate) was that shown in figure 7. Use of 
the values indicated by the full curve thus permitted calculation of the 
factor by which the measured velocity outside the wake should be re- 
duced to correspond with that of the infinite case. (Use of the composite 
values, it should be noted, will lead to an approximate stream function for 
the infinite case.) 


‘The magnitude of the resulting corrections may be judged from figure 8 a, 


in which are superposed the longitudinal components of the measured 


velocities (the points themselves) and the corrected values (the distribution 
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curves) at every other measured section. Also shown are the streamlines 
determined from the corrected distribution curves in the region of the eddy 
From these it will be seen that it was actually of littke moment whether or 
not the correction (based upon the assumption of irrotational flow) was made 
in the rotational region, where the velocity was already quite low. As a 
matter of fact, the velocity correction was so small in all zones as to have 
almost no effect upon the flow pattern itself. Its effect upon the pressure 


distribution, on the other hand, was considerable 








Figure 7. Distribution of doublet strength for eddy and wake displacement 


‘ ot hgure 6 


Figure 86 reproduces the actual measurements of the pressure intensity 
as points and the corrected values as continuous curves plotted to the left 
of the respective sections (i.e., negatively). Once again the correction 
was strictly applicable only outside the wake, for it was based upon the 
assumption of irrotational flow and constancy of the sum }p(u?+7?)+p. 
Unlike the velocity correction, however, that of the pressure could by no 
means be ignored within the wake, since there it was evidently a maximum. 
For want of a more precise method of estimating it, the customary boundary 
layer assumption ot negligible pressure change over the distance 6 had to be 
made; in other words, the correction at any point in the eddy was assumed 

be the same as that just outside at the same cross-section. ‘That the 
probable error in this assumption was small may be concluded from the fact 
that the measurements themselves varied only slightly across the eddy. 

With the corrected magnitudes of the velocity and pressure at hand 
throughout the zone in question, it was a simple matter to compute and 
plot the distribution of total pressure (i.e., as would be read at the stag- 
nation opening of a Pitot tube) relative to that in the ambient fluid, which is 
also shown in figure 84. Comparison of this sequence of curves with the 
locus of points of maximum velocity will indicate that the latter corresponds 
rather closely to the approximate limit of the zone of irrotational flow. 
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Figure 8. Experimental data and corrected distributions of (a) velocity, 


) pressure, and (c) turbulence, behind test plate 


(he summary plots in figure 8c of the three components of the turbu- 
lent velocity and of the cross product of the two components in the plane 
‘f motion (obviously uncorrected for the wake effect) quite logically show 








Fil Fif owe | 
3 44, © 


@ ' 


ot*— ——-+ 





Experiments on two-dimensional flow over a normal wall 139 


the region of maximum turbulence and maximum shear to coincide at the 
outset with the mean streamline dividing the main flow from that in the 
zone of separation, remaining at about twice the wall height as the eddy 
terminates. Were data for additional characteristics of the turbulence at 
hand—in particular, the microscale—it would be readily possible to trace 
the transformation of the energy of the primary or mean flow through its 
secondary or turbulent stage to its final dissipation as heat (Hsu 1950). 
For the present, it must suffice to note that in the zone of maximum turbu- 
lence level the energy of the secondary motion is some 15”, of that of the 


primary 


Se - a 
. Saar 

ae: Sas f | ' 

08 -— =. a eee KH, —— + _ _{—__+__4 





* ; 
in oe ee ee It aesnd 
¥ UC 7 | | 
/ nstr m ace ownistrear Tace | 
h i Upstrea e } L ‘ 
: a aan aa cn aa ena ——— rl ar 
| 1 alith ++ +o | 
2 —_—_—+—_- SSS ee ee ee 4 4 ; 
Uc ¥ i 4 | 
' 
| 
0 4 +— - + balks = 
~ 10 8 06 04 02 Y, A r ? 


Figure 9. Distribution of pressure on test plate, with and without tail plate 


There remains to be considered here only the evaluation of the force F 
exerted by the flow upon the wall. From the direct measurements plotted 
in figure 9, the integral of the variable intensity on the front and the corrected 
constant intensity on the rear vields for the drag coefficient the magnitude 


FAL 


5 pur 


~ 


Phat this value is in accord with the change in flow pattern that the plate 
produced was demonstrated by evaluating the momentum equation 


(1 ue" )a(? tf? a(t | “a(z)+ 
' u>}) \h Jo pur h nx 4s h 


 Y [? A(u/ug) (7) _# AL «30, 


h pu> 


uh! xh) 
between a section far upstream (i.e. where u ity 1) and an arbitrary section 
downstream. For the latter the section x 4 = 4 was chosen, and the inte- 
gration was carried out to—in lieu of infinity— the limit y h = 100 (extra- 


polation beyond the zone of measurement being facilitated by use of the 
stream function based upon the distribution of m(t) for the composite 
profile plotted in figure 7). The successive integrals were found to have 
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the values — 2-29, 3-04, 0-05, and 0-00, respectively, with the result that 
Ce 1-40), which, as a residual value, is in unexpectedly good agreement 
with 1-38 


INTERPRETATION OF RESULTS 
()ne must conclude from this investigation that the use of measure- 
ments made on isolated cylindrical bodies to approximate conditions i: 
hich the wake is not free to oscillate can result in rather large errors. ‘lhe 
actual drag coefficient of the wall (C, 1-+), for example, is only two-thirds 


) 
that of the isolated plate (C,, = 2-1), and the relative pressure change 
\p \p= beh nd the wall (—0-57) is less than half the mean value (— 1-36) 
behind the plate. ‘That this corresponds to a considerable difference in the 
length of the mean eddy zone (in accordance with the Riabouchinsky and 
related theories (Birkhotf 1950) that the free-streamline curvature must 
increase to produce a pressure decrease in a wake or cavity) is verified by the 
comparative measurements made in the course of the present study, the 
standing eddy behind the wall being more than three times as long (/ h 17) 
as the homologous portion of the temporal mean pattern of flow past an 


2-3). 


isolated plate (Lb 

Such reduction in drag and pressure drop, and such lengthening of the 
average eddy zone, point to the fact that a major difference between the two 
cases lies in the considerable amount of mean-flow energy that is required 
to produce the vortices that are shed alternately into the oscillating wake. 
At low Reynolds numbers these vortices are known to persist for a con- 
siderable distance prior to their disintegration into turbulence; in such 
circumstances, the wake behind the isolated plate should be narrower than 
that behind the wall. At high Reynolds numbers, on the other hand, the 
vortices tend to disintegrate quite rapidly; the rate of lateral spread of the 
turbulent zone behind the plate then probably exceeds that behind the wall. 
Since at all Reynolds numbers beyond a few thousand the standing eddy 
behind the wall appears to maintain an essentially constant form, the same 
can reasonably be expected of its wake as a whole; and since both the scale 
and the energy of the turbulence generated at the periphery of the standing 
eddy are smaller than the corresponding quantities in the oscillating wake, 
the disturbance caused by the wall should invariably dissipate more rapidly. 

How great an effect the existence of boundary shear upstream from the 
well (deliberately set aside at the outset because of its undefined magnitude) 
would actually have upon these characteristics can only be surmised. ‘The 
complexity of its experimental evaluation is seen from the fact that under 
even the simplest probable conditions of velocity variation—for example, 
the semi-logarithmic distribution wv 1+ Blog y—at least two variables 
must be considered: in brief, a magnitude and a rate of change. Fortu- 
nately, a phenomenon analogous to that under discussion is found in flow 
through a pipe orifice, for which it is both logical and experimentally 
demonstrable that maintaining a constant mean velocity and varying the 


transverse distribution by roughening the pipe surface will reduce both the 
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pressure change and the jet contraction. ‘There is little basic difference 
between the axisymmetric orifice and the two-dimensional slot, and the 
slot differs even less from the boundary arrangement investigated in the 
present study. One may therefore conclude that an increase in the velo- 
city variation over a given vertical distance (say 2h, the velocity at level / 
thereby remaining constant) could be expected to cause a decrease in the 
drag of a wall, in the pressure reduction behind it, and in the size of the 
standing eddy which it produces. Until the extent of this effect has been 
determined experimentally, however, the prediction of such related phe- 
nomena as Cavitation in the eddies formed by boundary projections must 
continue according to the rule of thumb that the significant velocity is that 
at a boundary distance equal to the height of the projection. 


This investigation was undertaken by the first author under the 
direction of the second, in part as a thesis project for the degree of Master 
of Science at the State University of Iowa and in part as a sponsored project 
under Contract NS8onr—500 between the Iowa Institute of Hydraulic 
Research and the Office of Naval Research. In the course of the hydro- 
dynamic analysis of the flow pattern, considerable assistance was provided 
by Dr Louis Landweber of the Institute staff. 
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SUMMARY 


(his paper contains a_ theoretical investigation of the 
displacement etfect of a pitot tube in a shear flow. Viscosity is 
neglected throughout so that the vorticity field alone is considered 

It is first shown that a two-dimensional approach does not 
produce a large enough displacement effect because it does not 


include the stretching of vortex tubes that takes place around 


a three-dimensional pitot tube. ‘Then the three-dimensional 
problem is considered. A solution is obtained in the plane of 
symmetry ror a sphere Ina she 1] How "This solut mis found 
by making an assumption about the rate of stretching of vortex 
tubes perpendicular to the plane of symmetry and then considering 
the shear flow as a smail perturbation of a uniform flow. A 
solution in the plane of symmetry is sufficient to obtain the displace- 
ment ettect, which is found to be of the same order as the 
xpel rental result obtained DY Young & Maas (1936) for a 
conventional pitot tube. ‘The sphere may be considered to 
represent a conventional pitot tube (of slightly smaller diameter), 


so it is concluded that a large part of the displacement ettect 
of a pitot tube may be accounted for without the inclusion of 
viscosity, i.e. by consideration of the vorticity field alon 

lo a first approximation, the vorticity in the plane of symmetry 
is found to depend only on the distance from the centre of the 
sphere ° 

\n outline of shear flows past some two-dimensional bodies 
is given in an appendix. ‘The bodies considered are a circular 
cylinder and a two-dimensional ‘pitot-tube’ consisting of two 
parallel semi-infinite plates. 


1. INTRODUCTION 
Correct interpretation of the pressure measured by a pitot tube in 
a shear flow is important in the experimental study of boundary layers and 
wakes, whose velocity profiles are generally obtained from pitot traverses. 
he problem has been investigated experimentally by Young & Maas 
(1936) and their approach, outlined below, shows the nature of the problem. 


‘The corresponding theoretical problem is then discussed, and an approximate 
solution is obtained for a spherical ‘ pitot tube’. 
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The variation of total pressure across a wake, measured by a pitot tube, 
is shown in figure 1, together with the hypothetical curve which would 
be obtained using a tube of zero diameter. In effect, the pitot tube with 
its axis at B measures the total pressure at 4. Hence, if the distance 5 
were known, the actual distribution could be obtained from the measured 
one. The distance 5 must be a function of the shape of the pitot tube, its 
diameter D, the velocity and its derivatives, and the kinematic viscosity. 











f 
Ww 
a 
> 
7) 
”) 
Ww 
a 
a 
il 
< 
ro 
oO 
we 
——> § «—__ 
AB “ 
DISTANCE ACROSS WAKE 
Fig | Dhe mn t ta I iC Ss ng displ ’ 


Suppose now that the second and higher derivatives of the velocity are 
zero, or so small that they can be neglected; i.e. the velocity of the oncoming 
flow is 

U+ Ay’, (1.1) 


where U is the velocity on the axis of the pitot tube, and \’ is a rectangular 
coordinate perpendicular to the axis. Then, if we consider only 
geometrically similar tubes, in an inviscid fluid, 


6 = function of D, U, A, (1.2) 


or, trom dimensional considerations, 


D = function of A, (1.3) 


where K = ~—-: (1.4) 
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Using flat-nosed pitot tubes with internal diameters approximately 
equal to 0-6 of their external diameters, Young & Maas found that 
= 0-18sgn AK (1.5) 
i.e. the displacement is towards the region of higher velocity. ‘Vhis result 
contains a difficulty at A 0, where there is a discontinuity in the displace- 
ment effect. ‘This is not physically plausible since it implies that there is 
1 discontinuous change in the flow pattern when the pitot tube is moved 
a small distance from a point of zero shear. However, the result obtained 
by Young & Maas is not necessarily valid near A = 0. ‘Their values of 
8 D show increasing scatter as K decreases, and no values were obtained for 
K less than 0-025. Also, when AK is small, it is probably not always justified 
to neglect second and higher derivatives of the velocity. ‘Therefore, it is 
likely that the actual variation of 6/D with A is rapid, but firmite, when A is 
very small, and that 6/D tends towards some value, constant for a particular 
pitot tube geometry (0-18 for the shape of tube used by Young & Maas), 
at larger values of A. During recent work at Cambridge University, 
further experimental values have been obtained for the displacement etfect 
vhich are in fair agreement with the Young & Maas result. 
So far, no theoretical treatment of a pitot tube in a shear flow has been 
made. ‘Iwo classes of problems have been solved in which the oncoming 
flow has a uniform rate of shear: 
(a) those in which the body is cylindrical and the whole tlow is two- 
dimensional (Lamb 1932, p. 233, Lighthill, Reichardt 1954, ‘l'sien 
1943); some of these solutions are outlined in an appendix ; 

(b) those in which the body is cylindrical with generators perpendicular 
to the oncoming vorticity (Hawthorne 1954, Squire & Winter 1951). 

Other solutions involving bodies in non-uniform streams have been 
obtained, but these are not relevant to the pitot tube problem (e.g. Nagamatsu 
1951, Hawthorne & Martin 1955), 

In the completely two-dimensional problem, the displacement effect 
can be found easily, but it is considerably smaller than the experimental 
values obtained with an axially-symmetric pitottube. ‘This is not necessarily 
due to the omission of viscosity (Reichardt’s experiments show very good 
agreement with inviscid theory) but is probably due to an essential ditierence 
between the vorticity fields in shear flows past two and three-dimensional 
bodies. ‘Thus, consider a vortex tube moving towards the pitot tube. 
\t a large distance upstream, the vorticity is everywhere perpendicular 
to both the velocity and the direction of the shear, so that a vortex tube is 
initially a cylinder with its generators perpendicular to the stream direction. 
When the flow is completely two-dimensional, the generators remain 
perpendicular to the stream, so that the cross-section of a vortex tube is 
deformed but its area is unaltered. Hence the vorticity is constant along 
any streamline. In front of a three-dimensional object, however, there is 
a retardation of the fluid, while at large lateral distances the flow is almost 
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unattected by the body. ‘Therefore, the vortex tube must be stretched so 
that there is a decrease in cross-sectional area (by continuity) and hence an 
increase in vorticity (by Kelvin’s theorem). ‘The net result is an accumula- 
tion of vorticity in front of the body and, in consequence, larger streamline 
curvatures than occur in the two-dimensional case. ‘This makes possible 
a much larger displacement effect. For example, it is shown in this paper 
that, when the vorticity of the oncoming flow is small, a sphere produces, 
according to certain approximate assumptions, a displacement approxi- 
mately five times as large as that produced by a cylinder of the same radius. 

\nother feature introduced by the stretching of the vortex tubes around 
a three-dimensional object is a streamwise component of vorticity. In the 
plane of symmetry, the velocity field of this component seems to be 
predominantly downwards, which would tend to increase the displacement 
effect. This feature is not accounted for explicity in the analysis of this 
paper, and would have to be taken into account in any attempt to improve 
the present theory. 

The three-dimensional problem is much more difficult, because the 
full equations are complicated and non-linear. However, the streamline 
reaching the stagnation point must always lie in the plane of symmetry, 
x = 0, so that a solution in this plane is sufficient to determine the displace- 
ment effect. ‘The flow pattern in this plane is not independent of the rest 
of the flow pattern, but it can be shown to depend on one quantity only, 
namely (Cw/0z), 9, where w is the velocity in the 3-direction. As a first 
ittack on the problem, the equations of motion in the plane of symmetry 
are solved below by assigning an approximate value to (¢w/0s), 9. Consider 
again a vortex tube moving towards the body. It is clear that the behaviour 
of the flow pattern in the plane of symmetry is governed by the rate of 
stretching of the vortex tube perpendicular to this plane. ‘This quantity 
is the relevant component of the rate of strain tensor, i.e. 0w/éz. ‘Therefore, 
it is necessary to know (¢w/éz), _ y in order to solve the problem in the plane 
of symmetry. ‘This is not possible without a full solution of the three- 
dimensional problem, so an approximate value is sought. Since the 
vorticity in the plane of symmetry is in the z-direction, it seems reasonable 
to assume that, in this plane, dw/dz is independent of the vorticity. On 
this assumption, (cew/¢ez), _, must take its value in the irrotational flow. 

At this stage the problem is still non-linear : in fact, it involves the solution 
of three simultaneous partial differential equations. However, in practice 
the oncoming shear is often small, which suggests that linearization may 
be effected by considering the shear flow as a perturbation of the uniform 
How. Since small shear corresponds to small K, a solution is sought in 
the form of a series in ascending powers of K. ‘To be able to use this 
approach, it is necessary to know the solution when K is zero. ‘This solution 
is not available for a pitot tube of conventional shape but only for bodies 
with closed, rounded noses. Therefore, since the solution for a sphere 
in a uniform flow is simple and well known (e.g. Milne-Thomson 1949, 
p. +13), an attempt has been made to obtain a solution for a sphere in a 
shear flow. 


F.M. K 
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If the sphere itself were used as a pitot tube, by making a small pressure 
hole on the axis, the definition of 6 as the displacement of the stagnation 
streamline would not have the desired physical meaning. The pressure 
measured would be the static pressure at the orifice, which is less than the 
total pressure on the stagnation streamline. However, a pitot tube of this 
type would be unsatisfactory because of its sensitivity to yaw. Instead, 
one may think of the sphere as replacing a conventional pitot tube of slightly 
smaller outside diameter, so that the streamline reaching the stagnation 
point is the same as the one reaching the stagnation point inside the pitot 
tube. ‘Then the displacement effect of the sphere is defined in the same 
way as for an open body; 1.e. the asymptotic distance of the stagnation 
streamline from the axis. ‘lhe above argument is illustrated by the two- 
dimensional case (see the appendix) in which the displacement effects of 
a cylinder and a parallel plate pitot tube are equal when the diameter of the 
cylinder is \ 2 times the spacing of the plates. 

By expanding the full equations of motion in powers of K, it can be 
shown that the assumption for (ew/dz),_»9 gives the vorticity correct to 
the first order in A, but omits an irrotational contribution to the velocity 
components. Unfortunately, this cannot be remedied without solving 
the problem for the entire three-dimensional flow pattern*. However, 
in order to obtain other than a linear variation for the displacement effect, 
results have been obtained up to the order of K* using the above assumption 
for (Cw/dz), 9. “The cubic term produces a maximum value of 6/D which 
is of the same order as the asymptotic values found experimentally and 
which occurs within the possible working range of values of K. 

‘lhe result for the first approximation to the vorticity is of some interest 
since it shows that it is a function of the distance from the centre of 
the sphere only, and that it has a square root infinity on the surface. In fact, 
the first approximation is (1 —7~*)-"? times the oncoming vorticity, where r 
is the distance from the centre of the sphere divided by its radius. 

The mathematical formulation of the problem, using the above 
assumptions, is given in § 2, and the equations are solvedin§3. ‘lhe solution 
is summarized in $4, and the velocity components are obtained. In $5, 
differential equations for the streamlines are found and solved for the 
stagnation streamline. Outlines of some two-dimensional solutions are 


given in an appendix. 


2. MATHEMATICAL FORMULATION OF THE PROBLEM 


Consider a cylindrical coordinate system r’, 0, 2’ with its origin at the 
centre of a sphere of radius a. ‘The undisturbed flow is independent of 2’, 
so that s’ = 0 is the plane of symmetry of the flow pattern, and is in the 
direction 6 = 7 (figure 2). Non-dimensional coordinates r, 3, are detined 
asx ld. 2% 1a. 

Since this work was compieted, Lighthill has solved the first order problem 
over the entire flow field. He obtains the same order of magnitude for the displace- 


ment effect and for the downwash velocity ahead of the sphere. This increases 


one’s confidence in the cubic term obtained in this paper 
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In the plane of symmetry, the vorticity components are (0, 0, ©) and the 
velocity is (w,v,0) where w and v are even functions of x. Therefore, in 


this plane, Helmhoitz’s equations of motion reduce to 


bo 
~ 


where = e wa ae 
‘The equation of continuity is 
nfo i {2 (2.3 
ro roé es) o a 
Equations (2.1) to (2.3) are the same as those for a two-dimensional rotational 


flow except for the additional terms containing (dw/és),_,. If a value is 
inserted for (Cw/dz), _ 9, the problem is reduced, in effect, toa two-dimensional 








Figure 2. The plane of symmetry of the flow pattern. 


one, and the above equations are sufficient to determine u, v, and ¢, provided 
that the boundary conditions are suitable. In $1, it was shown that a good 


approximation to (dw/Cz). _ 9 is its value in the irrotational flow. ‘Therefore 


(Milne-Thomson 1949), 
ow 3U cos 
3s), =0 . —— * (4.4) 


ie) 


: CG o ‘ us lod ?. ous 
ae (+x ze 
} ; A) 
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bo 
bo 


equations ( 1 (2.3) become 


- Lay a 1 3 
where V2? = -~(r—)+=. (2.8) 
rer\. er r? cA 


Ihe oncoming flow is parallel and has uniform rate of shear. Hence, 


the boundary conditions are 
(a) atr 0, C A, 
—u = Ucosé#+ Aarsin#@ cos 6, 
= Usin6@+ Aarsin?9, 
(b) ati L w=: 


However, when 4 is zero, the solution is known (Milne-Thomson 1949, 
p. +13). Thus, ¢ satisfies the equations 


“dh : l 

: cos 4 ( 1——}, 
) \ r 
ma | 

— = sind (1. —}. 
ri 7p 


it this solution for ¢ satisfies (2.6) when 4 is not zero.* Hence the 


? Y) 


2) 
rd 


problem is reduced to that of solving (2.1) and (2.7) with boundary conditions 


(a) atr = 0, {= A, 
ys = —Ir* sin*s, + (2.10) 
(6) atr=1, ws = constant. 


The solution is obtained by inserting expansions for ¢ and % in powers 
of K. These are 


(2.11) 


5 = OtKGtKa+.... 


On substituting (2.11) in (2.1) and (2.7) two series of differential equations 
ire obtained. ‘These are 


dy cl Cd el 3 COS 0 ‘ ! ; . - >) >) 
or } rene oe 2" 9 (Yos- ++ Baas Sore +e Sua) (2.12) 


V4b, = —Las (2.13) 


* It would not do so if the term of order K in (¢z/ 22), were included. 
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where the functions g,, are given by 


g,=W 


Sv ’ 

cus cl cus cl.\ 
.- ¥ (3.2 
; cr rc@ reO en 


‘The boundary conditions are 


) 

) 
to 
— 
+ 

mee 


(a) atr Ly Co = I, by — fr? sin?6, 
t () i; () (7 0), (2.15) 
(b) atr l, ub constant. 


In $3, equations (2.12) and (2.13) are solved for m = 0, 1, and 2. First, 
equation (2.12) with 7 = 0 is solved to obtain the first approximation to the 
vorticity. ‘This is substituted in (2.13) to obtain the first approximation 

) 


tow. ‘Then, substituting back in (2.12) with n = 1, the second approxima- 
tion to the vorticitv is found This is inserted in (2.13), and so on. 


Physically, this corresponds to finding first the vorticity field obtained by 


convecting initially constant vorticity along the irrotational streamlines, 
and then the velocity perturbations caused by this vorticity field. On 
substituting back into (2.12), the vorticity field that is convected along the 
perturbed streamlines is found. ‘Then the further velocity perturbation 
is obtained, and the process is continued until the desired approximation is 
reached. 

3. SOLUTION OF THE DIFFERENTIAL EQUATIONS 


Ihe general solutions of equations (2.12) and (2.13) are found below 
and then the special solutions satisfying the boundary conditions (2.15) 
are determined. 


Substituting from (2.9), equation (2.12) is 


[ 1 ec { | cl 3cos4 
“os 4 amy age ‘a i Ret ee Ee Se ag Lg 3 
cos 4 | | 4 = sin 4 | 3) <a, +5 (3.1) 


The general solution of (3.1) is the sum of a particular integral, which 


depends on g,, and the complementary function 7, which is the solution 


ot 


I\cn  . 1 \e 3 cos 4 
cosé{1— =.) > +sin@{ 1+ =)— =—— 7. (3.2) 
rye 2r3 }rcd _ 





Changing the independent variables from r, 6 tor, y, (3.2) becomes 


-eo(1-3)(2), +{ ene(1- 3) 


j 
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lherefore, if y 1s an integral of 


dr rdé 7 
——_———_-_—___ —————— (3.4) 
cos (1 — r-) sin 6(1 + 31-3)’ 
rP—1\12 , es 
g . ( — ) sin @, 5) 
\ } 


equation (3.3) reduces to 


whose solution 1s 


p3 12 om 
(x{=—) (i 
KX (= i) 


where / is an arbitrary function, determined by the boundary condition. 
Equation (2.13) is 

Vb, = —C,,. (3.8) 
Chis is a Poisson-type equation whose general solution is the sum ofa 
particular integral, depending on ¢,, and the complementary function, 
which is independent of ¢,. ‘The latter is a solution of Laplace’s equation 
in two dimensions, and is 
ay > [(a, 7° +b, r-")cos nb +(c, rr" +d, r-")sin n§], (3.9) 


where a,, 4,, c,, d, are constants determined by the boundary conditions. 


First approximation to the vorticity 
gy is zero. ‘Therefore, when n=O, the solution of (3.1) is the 
complementary function alone, given by (3.7), i.e. 





: ( r3 12 
r ae . (3.10 
>0 HX) r3 i) 1) 
Since l at? OQ, 
Hy) = I (3.11) 
lheretore, 
rs \12 ue 
So = (= i) ° (382) 


This equation shows that, to a first approximation, ¢ is independent of 6. 
This is because, for the particular case of irrotational flow past a sphere, 
u and (dw/ez),_, have the same variation with 6. ‘This can only occur 
when the body can be represented by a doublet at the origin, since this is 
the only case in which the velocity potential is the product of two terms, 
one of which is a function of @ only and the other is independent of #. 

Equation (3.12) also shows that the vorticity becomes infinite on the 
sphere. This is to be expected from the previous considerations regarding 
the stretching of a vortex tube. In fact, the vorticity must become infinite 


on any body with a stagnation point. 
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First approximation to ud 


Substituting (3.12) in (3.8) gives 





A particular integral is %, where 


l1d/ d¥\ re \12 . 
= A oo do ~(=3) le 





Hence, 


JW Ey , : a 
rs = — | x°2(x3—1)-)? dx. (3.15) 
ar ; 


To ensure convergence of the integral as r — x, (3.15) is rewritten as 


jy z 
a at ess ae > 
i by? + | {x92(x3—1)-12— x} dx. (3.16) 
ar 
Integrating again, 
Y= —1,-O, (3.17) 
dae ; — 
where O=; —| {x°%(x3-—1)-)2-—-x} dx, (3.18) 
§ 
r - 2 


The solution of (3.13) is the sum of Y and the complementary function 
(3.9). Hence, applying the boundary conditions (2.15), 
l 1 \ 
by = 4 r— 3 cos 20 — 37° O. (3.21) 
r- ~ 


Second approximation to ¢ 


Substituting (3.12) and (3.21) in (2.14) gives 


g, = £;(r)sin # cos 6, (3.22 
where a, = —37- 2A -— 1) — 1)°32, (3.23) 


Therefore, when n = 1, a particular integral of (3.1) is Z(r)sin @, where 


dZ Z ( Sa . 
ee ee ee 3.24 
Ar r ( I 3) £ ’ {9 ) 
A ee ; 3 
i Z=—s) x '%x'-1)(x3-1)-°? dx (3.25) 





152 


The solution of (3.1) is 


r3 12 : Mey 
G=i(a-q) +20nsina, 3 
r? — 
where y is given by (3.5). Then, since ¢, = Oatr = %, 
(x) = 6, (Ss: 
and C¢, = Zr)sin 6. (3. 
Second approximation to ¢ 
Substituting (3.29) in (3.8) gives 
lo Cw 1 cd ; as os 
=e ') - — —— = —Zipain’, (3 
ror cr y? op" 
A particular integral is R(r)sin 4, where 
1 d/ dR Rk : Z 
2 <-2 ( 
rdr\ d r 


It is convenient, both f 
consider the integrals R 
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j 


i 


or algebraic manipulation an¢ 


, and R,, where 


/R 1? 
R.=r— -~R=- Z(x) dr 
di ¥ 
dR 
ate - R= -—r| Zx) dx. 
= adi 


The double integrals o 


reduced to single integer: 


R, 


where d(Z x) dx is known trom (3.24). 


Therefore, subtracting 


adding 3r/2 outside the 


n the right hand sides of (3.32 


> 


ils. Integrating by parts, (3.3 
d (Z\ 
= —1rPZ+ vi - ( dx, 
; 4) aN x 


From (3.26), 
d 


= 


3/2 from the integrand to ensu 


ZA \ 


\ 


integral, 


l | zt ae 3 H 
4) Y a ) a ss 4y’ 
where H= >| {x7%(x4—1)(x3-1)52-1) 
Similarly, 
R, Ra IZ lrG 
where G : | x3 2(x4-—1)(x3 — 1)-°? dy, 


for computation, 


(3.53 


) and (3.33) can 
) becomes 


+ 


re convergence, 





0) 


and 


.36) 
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As r—> ©, 
G ~ 3r}, (3.40) 
H ~ Br, (3.41) 
‘Therefore, as r > ~, 
R= a ~-}. (3.42) 


The function ys, is the sum of Rsin#@ and the complementary function 
(3.9). Since ys must be constant on r = 1, 


oe 
ys, =<R- —— Sing 4a, + = \a,(r" —r-")cosn6+c,(r” —1r—")sin n6}. 


(3.43) 
This solution cannot satisfy the boundary condition at infinity for any 
values of a,, ¢,. However, the part of the flow pattern in which deviations 
from the upstream conditions would be expected to be important is a strip 
|v] = |rsiné| <k, (3.44) 
where & is finite. In principle, the complete flow pattern could be obtained 
by joining three solutions; the present one in the strip |v| < k, and two 
solutions for |y| > & which would have to be found by another method. 
However, since the outer flows would be very nearly parallel, and would 
have little effect on the inner flow near y = 0,if k were taken to be large enough, 
it should be sufficient to satisfy the boundary condition at infinity in the 
strip |vy|< k only. In this strip, the boundary condition is satisfied if 
a,. c, are zero for all n, 1.e. if 


R11). 
fy =<R- - ; sin 4 (3.45) 


Third approximation to € 


The equation for g, may be written in the form 


Yo = go(r)cos 6 + g5(r)sin?6 cos 4, (3.46) 

” ._ 40 ; ae : say Tis 

where g, = $742 —9rZ ate 1/273 — 1) -3/2, (3.47) 
a a3 


g=(r-4)% -(+4)2 
> ie \ xr?) dr r =) 7 ” 


When a = 2, a particular integral of (3.1) is 
X(r) + Y(r)sin?é, 


dr 2r 


where -(1 : ) i ini X +25, (3.48) 
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1\ da) joe) id : 
At. -}— + (1+ 55 —=~—Y+83. (3.49) 
r?} dr Ly +r ~§=2F - 
Hence, solving these equations, 
rs \1/2 ‘ 
X = ( s—} X;,(r), (3.50 
oe | 1 ) ) 
Y x ply? — 1)"*Y (1), (3.51) 
where : moc) ae Naa, = 
; Ce (3-7) 2o(x) dx, (3.52) 
Nee ; 
Y,= | #°0F—Ip"e"(e) &. (3.53) 
Asr—> o, 
X ~ —4r-, (3.54) 
Y ~ }r-}, (3.55) 
Therefore, the solution of (3.1) for ¢,, satisfying ¢, = 0 atr = 2%, is 
Co = A(r)+ Y(r)sin76. (3.56) 
Third approximation to ys 
Substituting (3.56) in (3.8) gives 
| C ous, ] Cub, 7 . P “ - fa gata 
-=-(|r <a -—— = —X—]V+1Ycos20. (3.57) 
ror\ of Y? 6g" . 5 
A particular integral is 
S(r)+ T(r)cos 20, 
where ] | dS Ae 
: r = —X-1Y, 3.58 
y dr | ; \ ) 
| = dT 47° Ly 3.59 
-(r—-)—-—> = }Y. 3.5 
r dr — r* < \ ) 
Integration of (3.58) gives, after manipulation similar to the previous cases, 
: fe : ¥ 
S=—fr—| ((X+4Y)x+} Hog (5) dx. (3.60) 
Jr = ae 


The integration of (3.59) is carried out in the same manner as the integration 
of (3.31). ‘Two integrals, 7, and T,, are defined by 


dT } 
p 4! _op 
‘a wr as 
(3.61) 
dT 
1 _ at on 
, dr af. } 
‘These are found to be 
T, = 57+ gr? logr—r? fSx8Y —dx— $271} dx, (3.62) 
Ty = -3P | xVdx (3.63) 


~-}r asr>o. 
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Hence, asr> ©, 

 —— = ax 2 

f yee <a (3.64) 
Again, it is not possible to satisfy the boundary condition at infinity for 
all values of 6. However, the solution satisfying the condition y. = 0 
at infinity in the strip |y| <A, as well as the condition 4. = constant on 
y = 1, is 


os | 
bo = S+47T a) cos 24 + (> - Jos @ (3.65) 
r } 


4. ‘THE VELOCITY FIELD 


I'he results of $3 are summarised below: 
= = (y+ Ke, sind + Klay + A259 cos 26 + O( K?), (4.1) 


| = Po + Poo cos 20+ K¥,, sin + K2P s+ 
K2Y,, cos 0+ K2¥,o cos 20+ O(K%), (4.2) 


1\-12 
where t.4 (1. 2] a ee 








2 
bog = X+4Y, fen oo 
| 1 \ 
wy ee O w = 7 \ i 
no ? ’ 2 } > / é 
, ? - * 7) F (4.3) 
R(1) 
P= R-— P= S, 
} 
| I FC) 
os . Boa 2e 
a 3 } FE 
The integrals QO, R, S, etc. have been evaluated and are given in table 1. 
aaa tent : 7 | 
ir }° QO | R | oS | Ff 4 =H i Fs —Z .| 
i— — — = 
0-0 | 0-0000 0:5000 | x XL | 0:0000 | 0-0000 0-0000 | 
O-1 | 0-0500 04099 2-2005 | 00-8358 | 0-0000 | 0-0501 0-0050 | 
0-2 01003 | 0-4994 | 1-2540 | 0-4329 | 0-0003 | 0-1011 0-0202 | 
| o-3 | 6-250 0-4980 | 0-8373 | 0-2908 | 0-0014 00-1556 0-0464 
()-4 0-2006 0-4953 | 0-6320 0)-2274 (0)-0047 0-2177 | 0-0859 | 
0-5 2515 0-4909 0-5108 0-1917 0-0125 00-2947 | 00-1436 | 
0-6 0-3032 -4841 | 00-4319 0-1691 | 0O-0341 00-3999 | 0-2295 | 
0-7 0:3562 | 0:4745 | 0-3773 0-1532 00-0821 00-5617 00-3659 
O-S 00-4114 0-4609 00-3378 00-1401 0-2038 | O-8573 0-6133 
0-9 0-470] 0-4413 | 0-3081 0-1268 0:6450 | 16112 | 1-2219 
1-0 | 05377 | 0-4089 | 0-2818 | 0-1077 L L | 


Table 1 
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‘The velocity components wu and @ are obtained from (2.5), (2.9) and (4.2), 
in the forms 


= = —(1—r-4)cos 4 — 2Kr-¥,, sin 20 + KY, cos 6 - 
- K*{r t¥,, sin @+ 2r-!¥,,, sin 20+ O(K4), (4.4) 
i , _ {dP oo 
; = 4 beth = fe se oe 
7 (1+4r-*)siné—K e% 7 Cos 
ly’ iv, ‘&¥, i. : r 
~K*= sin ~ K*)——2 + — 2 ent 4 ore cos 26 >+O(K4). (4.5) 
di dr dr dr 


The velocity is always finite, and this suggests that the series (4.4) and 
(4.5) converge when K is small enough. 
The position of the stagnation point is easily found from (4.5). Onyr = 1 


; = 1:5000sin #— A'0-3505 + cos 26! + 0:9667 K2sin 6 + 
- K3/()-5262 — 0-6667 cos 4 + 0-5524 cos 26? + O( K4). (4.6) 
Put y = rsin#d., (4.7) 


‘Then, if v, is the value of y at the stagnation point, 
vy, = Ky,,+ K3v,, + O(K®), (4.8) 
and equation (4.6) gives 
() K(1-5000y,, — 1-3505) + 


+ K8(1-5000y,, + 2-0000y,, + 0-9667y,, + 0-4119) + OCK®), (4.9) 


whence Vie = 00-9004, 
(4.10) 
Ve. = — 1°9357 
5. STREAMLINES AND DISPLACEMENT EFFECT 
The streamlines cre integrals of the equation 
dr rdé rs 
= ——, (5.1) 
u v 


Substituting from (2.5), this can be written 


C us C us ‘ C d mm 
K (3 dr+-—= de) = —dr— =, dé, (5.2) 
or 0 ar 


C 
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or, using (2.9), 


1 \ 
where P= (1 -= )” py*. (5.4) 


Using y and ¢ as independent variables, (5.3) can be written 


ous _ ous cP | oP : 
K--dr+K=- dy = = = G+ = @, (5:5) 
r oy * 2y or 2y oy - 
1 
ey 1 oP) dy oY 1 oP 
~-CWw a\ _Oow " 
‘Ke —-— =} 5++Ka-5 =; =9, (5.6) 
cy 2y dy) dr Or 2y or 


which is a differential equation for y as a function of ¢ along a streamline. 
Now, ¢% can be expanded in the form 
| Apo . 2 rah ’ 5 7 
b = don t+ Von v2 + Ady, yt..., (5.7) 
where the ,,’s can be found by comparison with (4.2). “hen, assuming 
that on a streamline 


y = vo(r)+ Ky, (r) +... , (5.8) 


a series of linear differential equations for y,(r) is obtained. ‘These equations 


ire 
dy l dp 5 
— +- —y = AAP, Vo,---5 Vi-4)s (5.9) 
Par 2 drtt — A Yorr Jew) 
where ie =O, 
pe 1 ’ 
cw Woo o ay ~ 
Ny = tt 2 + Qo Vos, + (5.10) 
cr oF * seins? ( 


lhe solutions of (5.9) are 


Vo = Cy p??, (5.11) 
y= apt? |" pte, dr Cp (5.12) 


It can be seen that (5.11) is the equation for the streamlines in the irrotational 
flow, as would be expected from (5.8). 
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The streamline reaching the stagnation point is of particular interest 
since it gives the displacement effect. On this streamline, 

















Vo = Vo = My = =O. (5.13) 
Hence 
ny 0 
} 
CW 9 Oss) OP )) : dy, (5.14 
nN —_ \ a Vs — ( 2s no Vy Y41) 
e or eC) : a 
Since y, must be finite when r = 1, the coefficient of p~'? in (5.12) must be 
zero, 1.€. 
C= | pn. (5.15) 
1 
Hence y=Pu one. (5.16) 
= 5 20 
Hus 
—— ‘$2 casesiahgete 4, —4ps 
* nesses 
ollie nal. 
4 05 
ene j i i 4 
= v s 4 3 = { 
~ 


Figure 3, The stagnation steamline, (7) Ky, (r)— K*y,(7) — OUK®) 


Using (5.14), the functions y, may now be evaluated. ‘The results for y, 
and y, are shown graphically in figure 3. 
The displacement, 5, of the stagnation streamline at infinity is given by 


Ky,(«) + Kys(2%) + O(K?) (5.17) 


o 


1-2400K — 1:1752K3+ O(K?). 
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The variation of 5/2a with K is shown in figure 4. It can be seen that the 
displacement effect found in this investigation is of the same order of 
magnitude as the experimental result of Young & Maas. ‘The gradient of 
8/2a with K, when K is small, is somewhat less than would be desired for 
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Figure +. The displacement effect of a sphere: (a) first approximation neglecting 


O.K*); (6) next approximation neglecting O(K°); (c) experiment, flat-nosed 
pitot tube 


good agreement with experiment. However, as indicated in §1, the 
displacement effect of a sphere is probably the same as that of a pitot tube 
of smaller diameter. Therefore, the initial gradient of 6/2a would be 
expected to be larger for a conventional pitot tube than for a sphere. 

It may be concluded from the above results that the major part of the 
displacement effect of a pitot tube in a shear flow may be accounted for 
without the inclusion of viscosity; i.e. the vorticity field is of prime 
importance. 

This work was done while the author was at the Fluid Motion Laboratory, 
University of Manchester. The author wishes to thank Professor M. J. 
Lighthill for suggesting the problem and for his continued interest and 
advice during the preparation of this paper. 
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APPENDIX. SHEAR FLOW PAST TWO-DIMENSIONAL BODIES 


‘The flow pattern around a two dimensional body in a flow with constant 
rate of shear 4 is defined by a stream function ¢ which is a solution of 
V4 = —A, (A. 1) 


subject to the condition that 4% is constant on the surface of the body 
(e.g. Milne-Thomson 1949, p. 108). Since & is constant on any streamline 
the displacement 6 of the stagnation streamline is obtained directly from 
the value of 4% on the surface of the body. 

Solutions of (A. 1), and hence 6, are given below for a circular cylinder 


and for a parallel plate pitot tube. 


(a) The circular cylinder 
The centre of a circular cylinder of radius a is taken as the origin of 
coordinates 7, #, and the undisturbed flow has velocity U' + Arsin@ in the 
direction 4 z. ‘lhe stream function for this flow is (Lamb 1932, Tsien 
1943, Reichardt 1954) 
d a’ a’. 
h | 4(r? — a?) 4 }A( 2 “z eos 20 l (, sin A. (A. 2) 


? 
Since ys is zero on the cylinder, the displacement 5 is the value of rsin# 
since wZw 1s Zero on the evinder, the aisp acement 0 1s the value of rsin 


for which ¢ is zero at infinity, so that 


| 4a? — |} Ad? — US = 0. (A. 3) 
Hence 
> = K-Mf—144/(14+ }K%} (A. 4) 
d 
~1K asK-0. (A. 5) 


(he negative square root has been discarded in (A. +) because 6 must be 
zero When AK is zero. ‘The variation of 6/2a with K is shown in figure 5. 
An interesting deduction about the stagnation streamline can be made 
rom (A. 2). The term — |} Aa*r-? cos 24 represents the ettects of the image 
vorticity inside the cylinder, which tends to deflect the oncoming flow 
lownwards. But the doublet term Ua?r-!sin @ tends to deflect the flow 
iway from the axis. Hence, the stagnation streamline is the one on which 
these effects balance, i.e. the one on which 
2 1 


— 
— sin@ — cos 20 = 0. (A. 6) 





lheretore, if 6 is always small on the stagnation streamline, the latter is 
straight and parallel to the axis at a distance r sin # from it, where 


m 5 ‘ ~ 
: sin? = 1K, (A. /) 


in agreement with (A. 5). This result can also be obtained mathematically 
from (A. 2). 
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(6) The parallel-plate pitot tube 

This idealized pitot tube consists of two semi-infinite flat plates at 
y = +c, stretching from x = 0 to x = o, with no net mass flew between 
the plates (cf. Lamb 1932, p. 75). The undisturbed flow has velocity in 
the positive x-direction. 





8 
26 bs 
0-3 
L 
O2h- , 
0-1 = 











Figure 5. The displacement effect of two-dimensional bodies: (a) parallel plate 
pitot tube, distance 2a apart; (6) circular cylinder; (c) experimental result 
for flat-nosed axially symmetric pitot tube. 


Let the stream function be % when A is zero, where yf, is zero on the 
surface of the body. At infinity, Y,~ Uy. ‘Then, in the general case, 
the stream function is 


b= f+ 4A. (A. 8) 
The displacement 6 must satisfy the equation 
Us + 3A(5*—c?) = 0, (A. 9) 
whence 
5 =~ rs 
= K-'-14++/(1+K)} (A. 10) 
~tiK aK-—0O. (A. 11) 


Comparison of (A. 9) with (A. 3) shows that, for identical oncoming flows, 
the displacement effect of the cylinder is equal to that of the parallel plate 
pitot tube when 


a 
c¢ = — A. 12 
wy (A. 12) 


i.e. when the diameter of the cylinder is 1/2 times the spacing of the plates. 


F.M. L 








162 I. M. Hall 


Figure 5 shows that 6/2c always increases with K and only reaches 
0-18 when K = 0-83. Since the highest value of K used in the experiments 
of Young & Maas was 0-3, this solution does not afford an explanation of the 
displacement effect. 

The above analysis is only slightly modified when the two-dimensional 
pitot tube is of arbitrary shape, provided that it consists of two parallel plates 
at infinity. ‘These are at a distance 2c apart as before. In this case, the 
stream function is 


b = y+ FAY? +4, (A. 13) 


where ys, is the stream function when A is zero, and y, is an irrotational 
stream function included to satisfy the boundary conditions. But, on the 
body, x, can differ from —}Ac? over a finite region only. | Hence it must 
always take this value on the body, and the displacement effect is the same 


as before. 
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SUMMARY 


This paper considers the changes that occur in the character 
of short-crested sea waves when they are refracted by a 
shallowing depth of water. Besides a change in mean wave- 
length and direction there is also a change (usually an increase) in 
the mean length of the crests. If the waves approach obliquely 
they become skew, that is, the crests become staggered one behind 
another. 

When a short-crested sea is superposed on a long-crested 
swell, refraction tends to amplify the longer waves more than the 
shorter ones. This also produces an increase in the mean length 
of the crests. 

Numerical examples are given. 


INTRODUCTION 

It is often noticed that sea waves tend to become more regular as they 
approach a coast, and that the length of the crests perpendicular to the 
direction of propagation appears to increase. Jeffreys (1924) considered 
two possible explanations: (1) non-linear effects reduce the short-crested 
waves by breaking, while the long-crested waves, being of smaller amplitude, 
do not break until later ; (2) the reduction in depth has a greater magnifying 
effect on the long-crested waves than on the short-crested waves of the same 
period. On the basis of theoretical analysis, Jeffreys decided in favour of 
the first alternative. 

There are other possible causes, however. The most important of these, 
now to be considered, is the refraction of the waves on entering shallow water. 

It is well known that when waves enter shallow water their velocity and 
wavelength decrease ; also, if they approach obliquely, the direction of the 
crests changes so as to become nearly parallel to the shore. But in general 
the crest-length also changes. If the water surface were a simple sine-wave, 
the waves would be infinitely long-crested and would remain so after re- 
fraction. But in fact the sea surface must be thought of as the sum of a 
whole spectrum of harmonic components having different wavelengths and 
directions ; the ‘long-crestedness’ depends on the angular deviation of the 
wave components from the mean direction: the smaller the deviation, the 
greater is the ratio of mean crest-length to mean wavelength. Refraction 
thus affects the crest-length in the following two ways. 

L2 
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(1) Wave components of the same wavelength and different direction are 
‘collimated’ so as to become more nearly parallel to the coast and to one 
another ; this tends to diminish the angular deviation of the wave components 
and to make the waves more long-crested (§ 2(a)). 

(2) Wave components in the same direction but of different wavelength 
are refracted through different angles and so separated (as in the separation 
of white light by a prism); this tends to increase the angular deviation of the 
wave components and to make the waves more short-crested (§ 2 (6)). 

The two effects (collimation and separation) are illustrated in §2 by 
considering a sea made up of just two harmonic components. The full 
representation of the sea surface, which requires a continuous spectrum of 
harmonic components, is discussed in §3. In general, the mean wave- 
length, mean crest-length, etc., must be defined statistically (see Longuet- 
Higgins 1956a*). ‘The changes in these quantities are deduced in § 4+ and 
§5, and a numerical example is given in § 6. 

One of the more interesting consequences of refraction is that waves 
which were originally symmetrical about their mean direction become 
skew, i.e. the principal direction of the wave envelope is in a different 
direction to that of the crests. ‘The waves then tend to line up in an échelon 
pattern, each wave crest being displaced (on the average) sideways relative 
to the preceding crest. 

The change in the amplitude of the waves is also studied, and it is shown 
that the increase in amplitude may itself be considered as a result of refraction. 

Finally, in $7, we consider what happens when two distinct bands of 
swell, of different wavelength, are present in the spectrum. For example, 
swell from a distant storm may be present simultaneously with shorter waves 
due to local winds. In general, the longer waves are amplified more than the 
shorter waves, and so, since the longer waves are usually also more long- 
crested, there is an increase in the average crest-length. We calculate the 
extent of this increase and give a numerical example. 


1. REFRACTION OF A SIMPLE WAVE TRAIN 
Consider a single, infinitely long-crested train of waves of wavelength 
A= 27, w and period += 27, o approaching a straight coast at an angle @ to the 
normal (see figure 1(a)). If the depth of water varies gradually the wave- 
length will be related to the local depth h by the equation (see Lamb 1932, 


ch. 9 2 
) wtanhwh = - =w’, (1.1) 


where g denotes gravity and a prime (’) denotes the value in deep water. 


The direction @ will depend on the wavelength according to Snell’s law: 
sind A we’ 


=y>=-—, (1.2) 


sin@’ A Ww 


which may also be written 
wsin@ = w’ sin 6’ = constant, (1.3) 


that is, the wave-number parallel to the coast is a constant. 


* This paper will subsequently be referred to as I. 
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We shall find it useful to represent the waves by a point P in a wave- 
number diagram as in figure 1(b): the direction of OP is opposite to the 
direction of propagation of the waves, and the length of OP is w, the wave- 
number in the direction of propagation. It may easily be shown that the 
wave-number in any other direction is the projection of OP in that direction ; 
in particular the coordinates 

(u,v) = (wecos8, wsin@) (1.4) 


are the wave-numbers perpendicular and parallel to the coast. 
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(a) 
Figure 1. (a) A regular, long-crested train of waves approaching a straight coast ; 
(b) its representation in the wave-number diagram. 


Let P’=(u',v’), represent the wave train in deep water, and P=(u, v) 
the refracted waves. Then, by (1.1), 


cothwh > 1, (1.5) 





II 


and, by (1.3) 
v=0' (1.6) 


’ 


so that the éffect of refraction is to displace the representative point P away 
from the origin and parallel to the axis of u. 


2. ‘TWO LONG-CRESTED WAVE TRAINS 
Two effects of refraction can be simply seen by considering the sum of a 
pair of long-crested wave trains in the following two cases. 
(a) Wave trains having the same wave-number w and slightly different directions 
6 and 6+ 66 
The sum of two such trains is a short-crested pattern whose wave-number 
is w in a direction perpendicular to the crests and wéé along the crests (see 
figure 2(a)). The ratio of the wavelength to the crest-length is 60. Let 
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the two wave trains be represented by points P, and P, in the wave-number 
diagram (figure 2(b)). Then OP, = OP, =w and / P,OP, = 66. Also 
P,P, = w6, the wave-number along the crests. Now refraction displaces 
P, and P, parallel to the u-axis, and 

P,P, dvsec#  cos@’ 





1, (2.1) 


P}P; — dvsec@’ cos 8 i 


\ 


\ \ 
SAR 


\\\ 
a ~ 
0 





bu’ Pat oe 





SS) 








AAAS 


(a) (4) 


Figure 2. (a) The pattern of crests formed by two wave trains of equal wavelength 


and slightly different direction ; (6) their representation in the wave-number 
diagram. 


since 6< 6’. Therefore, refraction increases the length of the crests, 
though the wavelength is diminished. The angle between the component 
wave trains, which is equal to P,P, OP,, is also diminished. 
(b) Wave trains having the same direction 6 but different wave-numbers w and 
w+ dw 

The sum of two such waves is a pattern which is long-crested but of 
varying amplitude, the waves coming in ‘groups’ of length 7/5w. In the 
wave-number diagram the unrefracted waves are represented by two points 
P; and P on the same line through the origin (figure 3 (6)) and with 
P}P; = dw’. However, (if P, and P, represent the refracted waves) P,P, 
does not in general pass through O; for in that case we should have 
Sw/w = dv/v = dw’/w’ and so dw’ dw = w’ w, whereas in fact we have, on 
differentiating equation (1.1), 





ow ; w’ 
— = tanhwh+whsech*wh > tanh wh = —. (2.2) 
ow w 


Hence P,P, subtends a positive angle at the origin, and the refracted wave 
trains have differing directions. Thus the refracted pattern is short-crested, 
i.e, the crest-length is reduced. 
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Since the refracted waves are also of different wavelength, the resulting 
pattern is ‘staggered’ (see I). A similar effect has been pointed out by 
Stoneley (1935) in connection with seismic waves refracted at a discontinuity. 
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Figure 3. (a) The pattern of crests formed by two wave trains initially in the same 
direction and of slightly different wavelength ; (5) their representation in the 
wave-number diagram. 


3. GENERAL REPRESENTATION 

To describe actual, irregular sea waves we use the representation 
described in I. It is assumed that the waves in any given locality are 
statistically uniform, and that they contain wave components of every possible 
wavelength and direction. Corresponding to each wave number (u,v) we 
define the energy density E(u,v) such that Edudv is the energy associated 
with wave-numbers in a small region dudv of the (u,v)-plane. Thus the 
total energy per unit of the sea surface, for example, is given by 


p= [” [ E(u, v) dudz. (3.1) 


~—-D-+~—o 
Many of the frequency characteristics of the sea surface may be derived 
in terms of E (see I). Here we shall make the assumption that the spectrum 
is narrow, i.e. that the energy E is mainly clustered near the mean wave- 
number (u,v), defined by 
-_) -2 ea a) 
| (u—u)Edudv = | | (v—v)Edudv = 0. (3.2) 
+~—-D-~—-®@ - - a0 


The waves are then more or less long-crested, having a mean wave-number 
zw in the direction @ perpendicular to the crests, where 


(u, @) = (wecos§, wsin§) (3.3) 
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(see figure 4). Let coordinates (€,») be taken with origin at the mean point 
and with € andy measured along and perpendicular to OP, so that, approxi- 
mately, 

















&€=w-—w = bu, ] 
ee (3.4) 
7 = w(6—6) = w88, ; 
VU 
/ "4 
v . T ) 
aA Ee 
ey oS 
es 
ee 
6 > U 
O 
Figure 4. The effect of refraction on a continuous spectrum. 
say. We define the moments 

Hoo = | | EE dédn, | 

r °*¢ 
Hay = | | EnE dédn, + (3.5) 

- my —« j 

“oO f xX " 2 

Loo = | 7° E dédn, | 


an interpretation of which is as follows (see I). First, 

Hoo = pwr, (3.6) 
where v is the r.m.s. proportional spread in wave-number in a direction 
perpendicular to the crests ; v is inversely proportional to the average number 
of waves ina ‘group’. Secondly, 


=. © _ 
Hog = pewy?, (3.7) 
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where y is the r.m.s. angular deviation of the waves from their mean 
direction ; it may be shown that 


y= Ninin (3.8) 
N max 
where N,,,,, is the average number of times per unit distance that the surface 


crosses a line at mean level perpendicular to the crests, and N,,;, is the 
corresponding number for a line parallel to the crests. ‘Thus if A is the mean 
wavelength and « the mean crest-length, 


A = 2 | kK = 2 a Di r K. (3.9) 
Thirdly, 
Hy = pw(v?— yg, (3.10) 
where q is a measure of the skewness of the waves; it may be shown that 
q = }tan 2A, (3.11) 


where f is the angle between the mean direction of the waves and the principal 
direction of their envelope (see 1). 


4. CHANGES DUE TO REFRACTION 

If the wave amplitude is sufficiently small, the energy corresponding to 
each wave-number is refracted independently. Further, it is assumed that 
the energy reflected from the coast is negligible. Under these conditions, 
it may be shown (see Longuet-Higgins 1956b) that if E’(u,v) denotes the 
energy density before refraction and E(u,v) the energy density after re- 
fraction, then E is related simply to E’ by 

E(u,v) = E'(u',v’), (4.1) 
where (u,v) is related to (u’,v’) asin §1. Hence a contour of energy density 
E’ = constant is transformed into a contour E = constant (see figure 4), 
though the area enclosed by the contour may be changed, as also may be its 
shape. 

We shall consider how the parameters 1, v, y, g defined above will change 
on refraction of the waves. Let (£’, 7’) be coordinates defined for the un- 
refracted waves similarly to (€’,7’) for the refracted waves (equations (3.4)). 
Now, to the first order of small quantities, 





. dw . 
p ae 4.2 
Bw = 7 bw’, (4.2) 
where from (1.1) 
dw ] 
— = P 4.3 
du’ tanh wh+ wh sech?wh oe 
Also, on differentiating (1.2) logarithmically, we have 
= a = dw—dw’, (4.4) 





tan@ tan @’ 
and so, using (4.2) and (1.2), 
w tané 


r w dw Lemp t poe = 
wod = | — — —— |tan@dw’ + — —,w' d0’. (4.5) 
Ww dw w’ tandé 
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Equations (4.2) and (4.5) may be written 


0, , 
E = 06" +2407) 











‘ > (4.6) 
N= O91 F" +99 77's | 
where 
0&€ dw 0€ 0 ] 
oa = = Ce ay oy 
“ 0€ dw . On (4 7) 
7 
on w dw On w tané 
is ee ee NM, agg ere es | 
= vy = dw ““ én’ sw’ tan” } 
Therefore, 


(én) _ 
a(".7') = 211 %aq- (4.8) 


Over the narrow region of integration, 11, %12)%1,% 23 may be treated as 
constants. Hence, on changing the variable in (3.1) and remembering (4.1), 
we have 


Pe = yy Bae fe 5 (4.9) 
similarly, from (3.5) we find 
ee ) 
H29 = %11 %20(%11 M20)» | 
Heyy = %1 %oa(%11 %o1 Mao + % 11% Mia)» (4.10) 


- 2 , 9, ee. - , 
Mor = %11 %a0(%21 Mo9 + 2091 Xoo My1 +%22 Hoa)» | 
where p’ is the initial value of 1,and soon. Thus », g, y are transformed as 
follows : 





eae ae } 
7 = a) a"; j 
¢ 
, 
w* , 
(v? — y?)q = =z [411% Y 7 +11 Haa(v ?—y *)q'], f (4.11) 
qy’2 
PF a ee eee a ee, ee ee 
yo = —z [4g VP + 2091 Aaa(v °— 7 )q’ +222 y *], | 
w* 


(the bar in w has been dropped). We may assume that q’, the skewness of 
the unrefracted waves, is zero. Then, on inserting the values of «11, %19; %29 
in (4.11), we have 


eo. 

— ie 
(v?—y?)g = r(1—r)tanév? (4.12) 

tan?0 ' 

a8 a (1 ~7 tan yp? + fa I 

, (1—r)?tan?@1 cantg | 

where 

Mi (4.13) 


w dw’ 


and @ is given by Snell’s law (1.2). 
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5. Discussion 

The first of equations (4.12) gives 

v=rv’, (5.1) 

that is to say, the average number of waves in a group is divided by 7 on 
refraction. Now, from (4.13), the factor r may be written 


aw oa’ /w’ (5.2) 
= st a 
da/dw ~ do’ dw’ 





since, if the frequencies are preserved, og = o’. But o/wand do/dw are the 
phase-velocity and group-velocity respectively, whose ratio varies trom 
2 in deep water to 1 in shallow water (see Lamb 1932, ch. 9). Thus in 
shallow water r = 1/2, and hence the number of waves in a group is doubled 
when the waves reach shallow water. (This effect was noticed by Barber 
(1950) in the case when the waves approach the shore normally.) However, 
if the waves were recorded in time at a fixed point, the groups in the record 
would appear to have the same average number of waves, regardless of the 
depth of water. 

Equations (4.12) show that v and g depend only on v’, the initial spread 
in wave-number, but that y has a contribution from both v’ and y’. By 
setting y’ = 0 and v’ =0 in turn, we see that y may be either greater or less 
than y’. 

The average length of the crests is given by 


K = Aly = 2n/yw. (5.3) 
So from the third of equations (4.12) we have 
1 1 ey ee cos?@’ r 
=3 = ga (l—r?)w v *tan?6+ cost zB? (5.4) 


which shows that the mean crest-length, also, may be either increased or 
decreased by refraction. 

A particularly simple case is when 6’ is small, or when the wave crests 
are nearly parallel to the coast. ‘Then from (1.3), 


wb = w’'é’" (5.5) 
approximately, and from (4.12) and (5.4), 
9 w’? ’9 oe reel e 
Peay, tek (5.6) 


(assuming v6’ <y’). Thus the average crest-length remains unchanged, 
although the average wavelength and the angular deviation are both 
diminished. 
The wave amplitude 

The change in amplitude of the waves may be found as follows. If a 
denotes the r.m.s. wave amplitude, the mean energy per unit area of the sea 
surface is }pga*. But this is equal to yw, by definition. Therefore, from 


(4.9), a uw dw cos@’ (5.7) 


—— =-—- —@ Loo = —— 
- 1122 dw’ cosé 


(dw cosé’\l2 ws 
and a=a (Z a r (5.8) 
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Equation (5.7) is a generalization of the well-known formula (Burnside 1914) 
for the change in amplitude of waves approaching a coast normally; on 
setting @ = 0’ = 0, we obtain Burnside’s equation (ii). Since, in general 
cos@ > cos6’, we deduce that waves approaching a coast obliquely are 
amplified less than those approaching normally. 

To facilitate numerical calculation, we have plotted in figure 5 graphs of 


’ 


uw 
— = tanhwh (5.9) 
w 
and 
r = (1+2whcosech 2wh)"} (5.10) 
as functions of w’h=o*h g. ‘The amplitude function 
dw \12 
— = (rw/w’)'2 511 
(JE) = Cw) (5.11) 


is also shown. When w’h is large, these all tend to unity. 
When w’h is small, we have 


w’ dw\!2  — (w'h)-14 a 
— ~ (w’h)!?, v~}, —) ~+ S.12 
w asia, (=) V2 ( ) 





(as was first shown by Green, 1838). However, the present linear approxi- 
mation is only valid so long as 


a <tu*h' (5.13) 


(Stokes 1847, Ursell 1953), which implies that 
cos @ 


ee 2 
a«< 3 (x h) h—0"" (5.14) 





If this condition is not satisfied, the waves will depart from the sinusoidal 
form, and may become solitary waves or may break. 


6. NUMERICAL EXAMPLE 
Suppose that waves of mean period 10sec approach the shore initially at 
an angle of 60° to the normal : 
2a - o ns 
o = —sec!, w = — = 0-012 ft}, 0’ = 60°. 
10 g 
When the waves have reached water of depth h = 20ft we have w’h = 0-25 
and so, from figure 5, 
w’ ss dw \\2 " 
— = 0-48, r = 0:54, ——}) = 1-07. 
w dw 
The angle 6 between the wave crests and the shore is given by 


sin@ = 0-48 sin 60° = 0:42, 6 = 24°. 


Let the initial spread in wave-number be one-quarter of the mean wave- 
number, and the initial spread in direction be 15°: 


vy’ = 0-25, y’ = 0-26. 











The refraction of sea waves in shallow water 173 


Then, assuming that q’ = 0, we have from (4.12) 
v = 0-135, y = 0-088, q = 0-68. 


Therefore, the angular spread is reduced by a factor y/y’ = 0-30. 
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Figure 5. Graphs of the functions w’/w, r and (dw/dw’)*/?. 


Taking account of the reduction in mean wavelength(w’/w = 0-48), wesee 
that the mean length of the crests is multiplied by 1-6, i.e. increased by 60%. 
From (3.10), we find for the angle of skewness 
B = 27°, 
and from (5.8), 
a = 0-79a’. 
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Thus the r.m.s. wave amplitude is reduced by 21%, although if the waves 
approached the coast normally it would be increased by 7%. 
Equation (5.14) shows that, for the analysis to be valid, the amplitude a’ 


in deep water must be small compared with 11 ft. 


7. A MIXED SPECTRUM 


So far we have supposed that the distribution of energy is narrowly 
restricted in wave-number and direction. It is interesting, however, to 
consider the changes due to refraction when the sea consists of a long- 
crested swell superposed on short-crested local waves. Under these 
conditions, the swell may be expected to be amplified more than the local 
waves, producing an increase in the length of the crests. 

Let the energy function E consist of two parts, F, and E,, which represent 
respectively a narrow band of swell, of total energy , and a local, fairly 
short-crested, sea of total energy pu. 

For simplicity, we suppose that both systems of waves are approaching 
the coast normally. 

The average number of zero-crossings along a line parallel to the crests 
will be denoted simply by N (so that V = N,,,). It can be shown (see I) 
that 


pitw2[° { eam. (7.1) 


9 


a eee eee 


So, if NV, and \, denote the corresponding numbers for the two component 


systems, 


2 72 72) 
pN? = wy Ni+peN3, | 
(7.2) 
= Py TP, | 
giving 
es Ly N2 + poN2 "4 
NM. See (7.3) 


"pd RL as 


or, if a,, a. are the corresponding r.m.s. wave amplitudes, 


Ne = atN?+ azN2 (7.4) 


The amplitudes and wave-numbers are transformed (approximately) as 


follows: 
, \1/2 
Wow ae, a, = (2) a, | 


dw 


dw, \12 
r oad 3 
No = 9» a, = ( a5, 


dw. ? 


re 
ay + a3 


J 


where w,, &}, W., ws are the values of w and w’ for the two groups. 
These relations enable us to determine N in terms of the initial values of 
Q,, @, N, and Ng. 
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For example, suppose that the mean periods of the two bands are 12sec 
and 6sec, respectively, and their r.m.s. angular widths are 10° and 20°. 
Then we have 

w, = 85 x 10-*ft-}, N’ = 4-7x 10-* ft-}, 
w, = 3-4x 10-* ft“), Nz = 3-8x 10-* ft, 
and, if the wave amplitudes are initially equal, the crest-length is given by 
N’? = }(N2+ N32) = 7:2x 10-*ft-?. 


By the time the waves have reached water of depth 10 ft. the amplitudes of 
the swell and the local sea are increased by 1-35 and 1-01 respectively, and 


we find 
N?2 = 6:3 x 10-6 ft-. 


Hence the mean crest-length is multiplied by 1-07, i.e. increased by 7%. 


8. CONCLUSIONS 
We have shown that, when a band of waves which has a fairly narrow 
spectrum in both wavelength and direction enters shallow water, the 
following changes occur. 


(1) The mean wavelength decreases and the direction of the crests 
becomes more nearly parallel to the shore. 


(2) The crests become ‘staggered’, i.e. the envelope of the waves is 
in a different direction to the direction of the crests. 


(3) The mean crest-length may either increase or decrease, usually the 
former : if the waves are incident normally, the average crest-length remains 
unaltered. ‘The ratio of crest-length to wavelength almost always increases. 


(4) The number of waves in a group is increased, being multiplied by 2 
when the waves reach shallow water. 


(5) The amplitude of the waves first diminishes and then increases, 
being proportional to h-!4 in shallow water. The increase is less for waves 
approaching obliquely. 

If the spectrum of the waves is not narrow, or if more than one narrow 
band is present, we may expect that the longer waves will eventually be 
amplified more than the shorter ones. Since the longer waves usually 
have a much greater crest-length than the shorter waves, this leads to an 
increase in the crest-length. 

The above conclusions suggest two possible reasons for the frequently 
observed change in crest-length as the waves approach the shore. The 
first is the simple effect of refraction, as described in § 4 and §5. The 
second is the relative amplification of a long band of swell, which might be 
hidden in deeper water off-shore. 

Our quantitative analysis is applicable only so long as the waves are 
fairly low (see equation (5.14)), and does not apply to waves in shallow water 
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which have steepened at their crests, still less to breaking waves. For 
these, some of the effects described may still be true qualitatively, but non- 
linear effects such as those pointed out by Jeffreys (1924) must also be taken 
into account. 


APPENDIX 

Some of the effects described in the present paper are clearly illustrated 
by the aerial photograph shown in plate 1, for which I am indebted to 
Mr W. W. Williams. ‘This shows a very regular swell approaching a 
coastline in a gradually diminishing depth of water. Not only do the 
mean wavelength and direction vary, but it will be seen that the wave 
crests arrange themselves in an échelon pa.tern with an angle of skewness 
that becomes very pronounced inshore. Lines of low wave amplitude 
may be traced visually over limited distances, and are not unlike those in 
figure 3 (a). 
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Plate 1. Aerial photograph of regular swell approaching a coastline obliquely. 
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On steady laminar flow with closed streamlines at large 
Reynolds number 


By G. K. BATCHELOR 


Cavendish Laboratory, Cambridge 
(Re Celve d 19 Dece nobe r 1955) 


SUMMARY 

Frictionless flows with finite vorticity are usually made deter 
minate by the imposition of boundary conditions specifying the dis- 
tribution of vorticity ‘at infinity’. No such boundary conditions 
are available in the case of flows with closed streamlines, and the 
velocity distributions in regions where viscous forces are small (the 
Reynolds number of the flow being assumed large) cannot be made 
determinate by considerations of the fluid as inviscid. It is shown 
that if the motion is to be exactly steady there is an integral con- 
dition, arising from the existence of viscous forces, which must be 
satisfied by the vorticity distribution no matter how small the vis- 
cosity may be. ‘This condition states that the contribution from 
viscous forces to the rate of change of circulation round any stream- 
line must be identically zero. (In cases in which the vortex lines 
are also closed, there is a similar condition concerning the circu- 
lation round vortex lines.) 

The inviscid flow equations are then combined with this integral 
condition in cases for which typical streamlines lie entirely in the 
region of small viscous forces. In two-dimensional closed flows, 
the vorticity is found to be uniform in a connected region of small 
viscous forces, with a value which remains to be determined—as is 
done explicitly in one simple case—by the condition that the vis- 
cous boundary layer surrounding this region must also be in 
steady motion. Analogous results are obtained for rotationally 
symmetric flows without azimuthal swirl, and for a certain class of 
flows with swirf having no interior boundary to the streamlines in 
an axial plane, the latter case requiring use of the fact that the vortex 
lines are alsoclosed. In all these cases, the results are such that the 
Bernoulli constant, or ‘total head’, varies linearly with the appro- 
priate stream function, and the effect of viscosity on the rate of 
change of vorticity at any point vanishes identically. 


1. GENERAL REMARKS 
‘The work described herein concerns the steady laminar motion of fluids 
of small viscosity, and is based on the generally accepted premise that, when 
the Reynolds number ofa flow field is very large, viscous forces acting on the 
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fluid are small everywhere, except perhaps in the neighbourhood of certain 
surfaces in the fluid. ‘ Small viscous forces’ here means forces that are small 
compared with unity when made non-dimensional using a length and a 
velocity typical of the flow as a whole (those used in the definition of Reynolds 
number would be suitable). It will usually happen that pressure forces on 
the fluid are of order unity when made non-dimensional in this same way, 
and the above premise is then equivalent to the statement that viscous forces 
are small compared with pressure forces nearly everywhere. 

It is well known, partly as a matter of observation and partly from mathe- 
matical analysis, that, for certain steady flow fields amenable to study, the 
above premise is undoubtedly true. It seems that, if the Reynolds number 
of the flow is allowed to approach infinity, without any other change in the 
conditions of these flow fields, the region of the fluid in which viscous forces 
are not small becomes smaller and smaller, and ultimately reduces, at most, 
to a number of thin layers, usually in the form of boundary layers and wakes. 
In the limit of infinite Reynolds number, the region in which viscous forces 
are not small either disappears altogether or becomes a number of stream- 
surfaces (i.e. surfaces whose tangent planes everywhere contain the local 
velocity vector), which usually coincide with, or are connected with, rigid 
boundaries in the fluid. Across such singular stream-surtaces there may 
be a discontinuity in velocity, as for instance at a rigid boundary where the 
presence of a singular stream-surface ensures that the no-slip condition is 
satistied even in the limit of zero viscosity. In what follows, the above- 
mentioned premise will simply be accepted as valid generally. 

We shall consider those steady flows that take place in a confined region, 
the motion of the fluid being generated by steady tangential motion of the 
surrounding boundaries (which need not all be rigid). It will be assumed 
that the Reynolds number is large enough for the thicknesses of the asso- 
ciated viscous layers to be small compared with the linear dimensions of the 
region of fluid under consideration. It will also be assumed that the motion 
of the fluid is laminar, despite the high Reynolds number (which clearly will 
correspond with reality only if the velocity distribution happens to have 
strong inherent stability). Most of the flows of this type that are capable of 
practical reproduction involve plane rigid boundaries moving in their own 
planes, or rigid surfaces of revolution rotating about their axes. Cases in 
which part of the boundary is rigid and stationary, the remaining part of the 
boundary of the region of closed streamlines being a ‘free boundary layer’ 

for example, the motion in a cavity opening off a plane wall over which 
fluid is streaming—are of interest in a wide range of practical problems in 
aerodynamics and hydraulics. Examples of steady flow in a closed region 
are not common, but they perhaps occur often enough in practice to 
warrant notice of their peculiar features. 

The equations governing the steady laminar motion of a uniform in- 
compressible fluid are 


Va = 0, (1.1) 
u <w—V(p pt 1g?) +vVu = ou ot = 0, (1.2) 
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wherew = V xu, gq = |u|, and the other symbols have their usual meanings. 
When the Reynolds number of the motion is large, viscous forces, according 
to our premise, are small everywhere except in the neighbourhood of certain 
singular surfaces, and equation (1.2) reduces approximately to 


uxw = VH, (1.3) 


Where H = p/p + 3q? is the local ‘total head’ in the fluid. H is here constant 
over stream-vortex surfaces, or ‘ Bernoulli surfaces ’, each such surface con- 
taining the local vectors u and w in its tangent plane everywhere and being 
swept out by a material vortex line. When i 0, stream-vortex surfaces 
will not exist, in general, because uw is then (see (1.2)) not necessarily 
parallel everywhere to the gradient of some scalar function. 

It is well known that equation (1.3), which is approximately valid every- 
where except in the neighbourhood of the singular surfaces, is not suth- 
cient to allow the velocity distribution to be determined from the condition 
of zero normal velocity at specified boundaries in the field. In the case of 
two-dimensional steady motion, the indeterminacy takes a form such that 
the vorticity—which is constant along a streamline—may vary arbitrarily 
from one streamline to another. It often happens that the inviscid flow 
equations can be made sufficient to determine u with the help of additional! 
information, for instance about the variation of w from one streamline to 
another far upstream (most often in the form of a statement that the velocity 
is uniform at infinity). 

However, in the case of flow in a confined region, for which all the stream- 
lines are necessarily closed, the possibility of making equation (1.3) sufficient 
to determine u in the region of small viscous forces, by introducing boundary 
conditions which specify w everywhere in a region ‘far upstream’, no longer 
exists. Other means of making the velocity distribution determinate must 
be found, and it is apparent that there is no further information to be found 
from considerations of the fluid as inviscid ; this is the feature that makes the 
study of flow with closed streamlines novel and interesting. It will be shown 
that the action of viscosity imposes certainly one, and, when the vortex lines 
also are closed, two, integral conditions on the distribution of w, and that in 
the cases of two-dimensional flow and of rotationally symmetric flow (with 
suitable restrictions) these integral conditions render the distributions of w 
and u determinate in the region of small viscous forces. ‘The question of 
what further conditions are needed to make the distributions of w and u 
determinate in general three-dimensional flow with closed streamlines is 
left unresolved. We begin with a derivation of the integral conditions in 
the general three-dimensional case. 

2. INTEGRAL CONDITIONS ARISING FROM THE EFFECT OF VISCOSITY 

To obtain a condition which arises from the effect of viscosity and which 
is valid no matter how small the value of vy may be, we operate on the com- 
plete equation of motion (1.2) in such a way that the contributions from all 
terms other than the term containing v vanish identically. Such an 
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operation is to take the line integral around a closed contour (in space) of 
certain shape, with line element dl, giving 


b(u xw).dl — } ANH — v bv xw).dl = $4 u.dl. (2.1) 


The term on the right-hand side vanishes in view of the steadiness of the 
motion, (It is important to notice that the motion is regarded as exactly 
steady, even though the value of v will later be taken as very small and even 
though steady motion would take a long time to develop from arbitrary 
initial conditions. In other words, we are considering flows subject to the 
double Jimiting operation t —- x, v —- 0, in that order, corresponding in effect, 
to the procedure that would be used in a real observation of the forms of the 
steady flows set up at several different large values of the Reynolds number.) 
‘The second term on the left-hand side vanishes since H is a single-valued 
function of position. We now make the first term on the left-hand side 
zero by choosing the closed contour to coincide with a streamline, of which 
a line element will be donated by ds. ‘Then, since v is non-zero (although it 
will later be assumed to be small), we have the exact integral condition 

(Vv xw).ds = 0, 2) 
to be satisfied for every closed streamline. 

It is natural to enquire if there are any other choices of the closed contour 


for which the integral (u xw).dl vanishes identically. If there existed a 


family of surfaces such that their normals were everywhere parallel to u xw, 
we could make the integral zero by choosing the contour as any closed curve 
on one of these surfaces. However, as already remarked, such a family of 
surfaces does not exist in general ; the surface formed by all the streamlines 
passing through a given vortex line will in general be intersected by other 
vortex lines (except when v = 0), which is irrelevant, since we are seeking an 
integral condition which is exact for small but non-zero values of v). For 
reasons related to the symmetry, surfaces which everywhere contain the 
local vectors u and w exist in cases of two-dimensional motion and of rota- 
tionally symmetric motion without azimuthal swirl when v = 0, but we 
shall see that the condition (2.2) alone is then sufficient to make the distri- 
butions of w and u determinate and no other choice of closed contour is needed. 

There exists the possibility that vortex lines are closed and that the 


integral @ (u xw).dl may be made zero by choosing the contour to coincide 


with a closed vortex line. Vortex lines.may end at a rigid boundary, but 
there is no reason why at least some of the vortex lines in a confined flow 
should not be closed; for such lines we have the additional exact integral 
condition P 


p (v xw).dv = 0, (2.3) 
where dv is an element of a vortex line. It appears to be difficult to decide 
whether vortex lines are, or are not, closed in any given flow field, but in one 
case, described in § 4, this is possible and the condition (2.3) is utilized. 
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Other contours which are closed, and which are such that each line 
element is orthogonal to u xw, may be imagined (one natural choice is a 
contour which consists of four segments, of which the first and third coincide 
with streamlines, and the second and fourth with vortex lines); however, 
their existence depends on the particular properties of the flow, and in any 
case it is not clear how the corresponding integral condition could be made 
use of. 

The integral condition (2.2) (and likewise (2.3), where it is applicable) is 
valid independently of the Reynolds number of the flow. If now the 
Reynolds number is assumed to be large; the equation of inviscid flow, (1.3), 
is approximately valid over nearly all the flow field, and there exists the 
possibility of making use of both equation (1.3) and the condition (2.2). 
Provided the streamline, around which the integral in (2.2) is taken, lies 
entirely in the region of small viscous forces, the integrand in (2.2) may be 
evaluated with the aid of (1.3). It will be shown that, in this way, the form 
of the velocity distribution in the region of small viscous forces may be 
determined in certain cases. 

‘The proviso of the penultimate sentence is equivalent to the requirement 
that the shortest distance from the streamline to any singular surface does 
not tend to zero as v-- 0. Now the velocity in the region of small viscous 
forces will in general be of the same order of magnitude as the tangential 
velocity of the boundaries (as can be verified experimentally in certain 
simple cases, and can in any case be examined a posteriori). Consequently, 
in cases of two-dimensional flow and of rotationally symmetric flow without 
azimuthal swirl, a typical streamline passing through the region of small 
viscous forces, on which the value of the appropriate stream function % 
differs trom that for the outer enclosing boundary or singular surface by a 
finite amount, cannot come close to the boundary of this region without the 
velocity there being infinite. ‘The above proviso is therefore satisfied for 
typical streamlines in these two cases. However, in more general types of 
flow, it is not certain that streamlines lying entirely in the region of small 
viscous forces exist; indeed, there are some closed flows for which all the 


closed streamlines pass through a boundary layer region. 


3. STEADY TWO-DIMENSIONAL FLOW WITH CLOSED STREAMLINES 
When the flow is two-dimensional, we can introduce the stream function 
us, and use (us, €) as orthogonal curvilinear coordinates, the lines € = const. 
being everywhere normal to the streamlines. ‘The displacements corre- 
sponding to increments in & and € are di/s/q and Ad€é, where / is an unknown 


function of % and €. The inviscid flow equation (1.3) can then be written as 
dH (us) 
ws) = ——. 3.1) 
) A ( 


‘This equation will be approximately valid, when the Reynolds number is 


large, everywhere except in the neighbourhood of certain singular curves 
(in the plane of motion) which are themselves members of the family of 
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streamlines. Even when the shapes of outer and inner streamlines (perhaps 
given by the shape of enclosing rigid boundaries) bounding a region in 
which (3.1) holds is given, there will in general be a solution of this equation 
for any choice of function w(), and the flow in the region of small viscous 
forces can be made determinate only with the aid of further conditions. The 
considerations given in §2 supply the integral condition (2.2), which we 
proceed to apply. 

When evaluating the integral in (2.2) for streamlines lying wholly in the 
region of small viscous forces, we can make use of the approximation (3.1), 
whence V xw becomes a vector parallel to the local streamline, and (2.2) 
takes the form 





dw rr 
a sS= (). 3.2 
7h p qds (3.2) 
ys 
Hence wb) = sic ) = w,(const.) (3.3) 
aw 


everywhere in a connected region of small viscous forces. (‘The possibility 
f there being an exception to (3.3) on a streamline for which 4 ds = 0 
—which is possible only if the velocity is zero at all points on the streamline 

with a possible discontinuity in w across such a streamline, is irrelevant 
since the viscous forces would then be large at such a streamline). The 
distribution of velocity in the region of small viscous forces can be deter- 
mined from (3.3) without difficulty when the shape of the streamline 
bounding this region is known*. 

The argument leading to the simple result (3.3) may be summarized as 
follows. In view of the fact that v is small, convection of vorticity will 
dominate viscous diffusion of vorticity when both processes occur, so that 
w is approximately constant along streamlines. But, in exactly steady 
motion, the net viscous diffusion of vorticity across a closed streamline 
must be exactly zero (even when the streamline encloses a rigid boundary), 
and this is then possible only if w is also approximately constant across 
streamlines. It seems that two-dimensional flows with closed streamlines 
cannot be exactly steady until the slow but persistent effect of viscous 
ditfusion of vorticity across streamlines has evened out any variation of 
vorticity that may have been present initially; the time required for this 
isymptotic steady state to be set up will of course increase as v decreases. 


* The notion of a two-dimensional ‘ inviscid’ core with uniform vorticity nas 
already been inferred from arguments rather less general or rigorous than those 
given above, and has been applied to some problems of two-dimensional free con- 
vectionin closed regions. Pillow (ina dissertation submitted for the degree of Ph.D. 
it the University of Cambridge, 1952) and Batchelor (1954) have used it to cal- 
culate the heat transfer across rectangular cavities, and Carrier (1953) has done so 
tor a circular cavity. The temperature in the core was shown to be uniform in all 
these cases from an argument based on symmetry of the streamlines, but in fact 
this result is true generally in a simply-connected *‘ non-conducting’ closed region 


n two dimensions, as may be seen from a proof like that used above 
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Ihe value of the constant wy is undetermined as yet, and since the distri- 
bution (3.3) is such that vV°u is identically zero there is no further inform- 
ation to be gained from considerations of the region in which viscous 
forces are small. Surrounding the region in which (3.3) is valid there is a 
singular streamline, in the neighbourhood of which viscous forces are not 
small, and the value of wy, will presumably be determined by the need for 
steady motion to be possible near this singular streamline. ‘The situation 
can be illustrated by reference to the simple case of flow inside a circular 
cylinder of radius a which rotates steadily with angular velocity }w,, an 
inner sleeve of length 27ac being held stationary (figure 1). It is evident 
here that, when the Reynolds number is large, viscous forces will be small 
everywhere except near the outer circular boundary, so that the region of 
‘inviscid’ flow is circular and in it the fluid rotates as a rigid body with 
angular velocity }w,. ‘The motion in the ‘inviscid’ core can be regarded 
as a standing eddy, which is driven by the motion of part of the outer boun- 
dary, the exact speed of rotation of the eddy being determined by the need 
for steadiness in the viscous boundary layer surrounding the eddy. 


fixed sleeve 


“a ae 
region of inviscid 


motion with uniform annular region in 
cylinder vorticity W, which viscous 
rotating forces are not 
steadily small 





Figure 1. A case of two-dimensional flow in a closed region. 


lhe present paper is concerned with general results rather than solutions 
for particular configurations of the boundaries, but it may be found illumi- 
nating to present some consideration of the boundary layer problem asso- 
ciated with the simple system shown in figure 1. Being unable to solve the 
boundary layer equation exactly, I at first solved the problem by linearizing 
it in the manner of Oseen* (which is equivalent to converting it to a time- 
dependent problem as Rayleigh did for a flat plate). However, it was 

* Professor H. B. Squire has also conceived this simple boundary layer problem 
and has solved it in this same way, in a paper to be published in the Journal of the 
Roval Aeronautical Society. 
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later pointed out to me by Mr W. W. Wood that so far as the relation between 
w, and w, is concerned the problem may be solved exactly by using the 
von Mises form of the boundary layer equation. ‘This form (see Modern 
Developments in Fluid Dynamics, Oxford University Press, 1938, Vol. 1, § 50) 
becomes 
Lc C(« *) 
“4 = ie it (3.4) 
h< € Fe ous? 
in cases in which the velocity outside the boundary layer is uniform, where 
hdé represents a displacement along the streamline ys = const., as before, and 
q is the velocity in this same direction. ‘Then, since the streamlines here are 
closed and q is single-valued, we have 
- 92/72 a2 
(9°) . ( ‘ ; ‘ 
0 p nae gh dé = =; ‘) gh dé, (3. 
( ys z ous" z 
since the mesh parameter / is approximately constant across the boundary 
layer. ‘The total variation of ; across the boundary layer tends to zero as 
v —- 0, so that the solution of (3.5) is effectively 
qth dé = const. (3.6) 


throughout the boundary layer. On evaluating the integral for the stream- 
line at the wall and for one just outside the boundary layer, we find the 
required relation to be 
Wo/w, = (1—a)!?. (3.7 
It is still necessary to rely on some approximate procedure like the Oseen 
linearization for details of the velocity distribution in the boundary layer at 
different values of €, but (3.7) represents the crucial piece of information. 
It should be noted that the basic assumption that the velocity in the region 
of small viscous forces does not tend to zero as v —~ 0 is confirmed in this case. 
Another problem whose examination involves a consideration of the 
exact shape of the boundaries, and which will not be attempted here, con- 
cerns the position of the singular streamlines. The location of these viscous 
layers will sometimes be evident from the nature of the conditions at the 
outer boundary, as in the very simple case represented in figure 1. Cases 
in which their location is not evident will clearly present great difficulties in 
any detailed analysis, akin to those in problems in which a boundary layer 
separates from a rigid wall. ‘he natural assumption that the singular sur- 
faces coincide everywhere with rigid boundaries (except where the contrary 
is evident) needs particularly careful scrutiny. For instance, in a case of 
flow in a region bounded externally by a rigid wall which has a 90° corner, 
a viscous boundary layer would not flow along the wall right up to the corner, 
in general, because there would then be a stagnation point of the inviscid 
flow at the corner; it is probable that a singular surface exists to divide the 
main body of the fluid from a secondary ‘standing eddy’ in the corner, and 
indeed there may be even a whole sequence of such singular surfaces and 


‘standing eddies’ of diminishing size as the corner is approached. 
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Finally, it is worth noting that steady two-dimensional motion of a 
fluid relative to given boundaries is unaffected by steady rotation of the 
whole system about an axis normal to the plane of motion (Taylor 1921). 
Consequently, all the above remarks and results apply not only to two- 
dimensional flows enclosed by outer boundaries whose positions are fixed 
and whose velocities are steady, but also to flows enclosed by boundaries 
whose positions and velocities are steady relative to suitably chosen uniformly 
rotating axes. 


4. STEADY ROTATIONALLY SYMMETRIC FLOW IN A CLOSED REGION 

It is convenient here to introduce orthogonal curvilinear coordinates 
(€,7, 45), where the €-lines (on which y and % are constant) are everywhere 
parallel to the component of u lying in an axial plane, the 7-lines are azimuthal 








axis of symmetry 


Figure 2. Coordinate svstem for rotationally symmetric flow. 


circles, and the ys-lines (4 being the Stokes stream function of the component 
of the motion in an axial plane) lie in axial planes and are orthogonal to the 
é-lines (figure 2). The displacements corresponding to increments in 
- I, « > 
¢,7,W are 6 
al h,dé, rdn, hyd, 
where h,(€, 7) is unknown, r is the distance from the axis of symmetry, and 
the definition of % supplies the relation 
a tone 
h, = (ru,y', (4.1) 

where (u1, Ug, 9) are the velocity components. 

The €-lines, or streamlines of the components of velocity in an axial 
plane, are bounded externally, and perhaps internally also, by a singular 
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curve which is itself a member of the family of €-lines. Provided the velo- 
city component uw, in the region of small viscous forces remains of order 
unity as vy —- 0) —which is evident in the case of flow without swirl, since the 
‘inviscid’ core would not otherwise be responding to the motion of the 
boundaries, but is not necessarily true in the case of flow with swirl since the 
motion generated in the ‘inviscid’ core may here be primarily an azimuthal 
swirl—a typical streamline will lie in the region of small viscous forces 
over the whole of its length. It will thus be possible to combine the exact 
integral condition (2.2) with the approximate inviscid flow equation (1.3). 
However, only in the case in which the azimuthal component of velocity is 
zero are these two equations sufficient to determine the form of the distri- 
butions of w and u in the region of small viscous forces. 


(a) Flow without azimuthal swirl 
The velocity here has components (u,, 0,0), and the components of the 
vorticity are (0, w»., 0) where 


1 o(hyu,) 


Re eS eS «6 
r hh, Cds 
Everywhere in the region of small viscous forces, (1.3) is satisfied approxi- 
mately, whence 1 dH(.) 


Uy Ws = h, a ’ 
ante: We _ dH (i) (4.2) 
r difs 
This last relation describes the known proportionality, in inviscid flow, 
between the vorticity and the length of a material element of a vortex line. 
This is as far as one can go, making use only of the equations for inviscid 
motion, and (4.3) is the counterpart of (3.1) in two-dimensional flow. For 
information about the function H’(ys), which describes the unknown vari- 
ation of vorticity across the streamlines, we must take, some account of the 
action of viscosity, and this will be done by applying the integral condition 
(2.2) to closed streamlines lying entirely within the region of small viscous 
forces. 
The components of V xw are given by 


seas & 2. — 0, i a), (4.3) 


th, ou 











and the condition (2.2) becomes 


—— hy dé _() 


ous rh, . 





or, in view of (4.1) and (4.2), 

" * or 2h 

H"() > Push, dé = -HW) OS dé 
a u ¢ 3 


= 2H (bp cos 6 ds 


© 


=0 (4.4) 


b] 
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where @ is the angle between the tangent to the streamline (€ and s increasing) 


and the axis of symmetry (with regard for the sense). Hence, we require 
H"(s) = 0, A(t) = xb, (4.5) 


where x is a constant (the other constant of integration being absorbed in 
the definition of ys), and 
s/f = If thy = « (4.6) 

everywhere in a connected region of small viscous forces (the possibility of 
a different result applying on certain exceptional streamlines again being 
irrelevant, since viscous forces would not be small at such streamlines). 

The result (4.6) is an obvious counterpart of the result that the vorticity 
is uniform in the region of small viscous forces in two-dimensional flow, and 
it is also true of (4.6) that the net effect of viscosity on the rate of change of 
the vorticity at any point vanishes identically. (But note that in axisym- 
metric motion the vorticity does not satisfy a heat-conduction type of equa- 
tion, and it does not seem possible here to arrive at the result (4.6) by an 
argument in terms of diffusion of vorticity across streamlines.) ‘There is 
also the common linear dependence of H on ¢, although the meanings of ¢ 
in the two cases are not the same. The result established for the region of 
small viscous forces in two-dimensional flow is such that the net viscous 
force on any element of fluid vanishes identically, but this is not true of (4.6). 
The local viscous force per unit mass of fluid is — vV xw, and it is readily seen 
from (4.3) and (4.6) that this is a uniform vector, of magnitude — 2x and 
directed along the axis of symmetry. ‘Thus the local viscous force has a 
simple character, and does not require the velocity distribution to be 
ditferent (by even a small amount) from that for an inviscid fluid, since the 
viscous force can be balanced exactly by a uniform pressure gradient. 

When the shape of the singular surfaces bounding the region of small 
viscous forces is known, the velocity distribution can be found from (4.6). 
‘The constant x, like w,, then remains to be determined from the condition 
that the surrounding viscous boundary layer is steady. 


(b) The general case (flow with swirl) 
The velocity components are now (i), és, 0), and the components of the 
vorticity are (@,, @, 3), Where 











7 1 (rus) 1 c(hyu,) _ 1 e(ruty) 47 
" th, op ” 2 Iyhg ob” 3 th, of wn 
Che approximate equation (1.3) then yields the three scalar relations 
1 oH Uy O(Uy) 
-_———_— Ee  - ,yW 4.8 
h, c& th, o€ ( ' 
1oH u, (rus) 
gee Oh an oe ee, 4.9 
r en rh, o€ ( 
i cH _.: (Ayu) | My O(a) (4.10) 
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From (4.9), we have 
ru, = C(s) (and ws, = 0), (4.11) 
which describes the constancy of circulation around’a material curve in the 
form of a circle about the axis of symmetry, and then (4.8) reduces to 
H = H(s). (4.12) 
‘The last scalar equation, (4.10), can be written as 
Wo dH 1 dC? 413) 
7” db Bap’ aie 
which is the generalization of (4.2). 

‘These relations appropriate to purely inviscid flow are now combined 
with the exact relation (2.2) which arises from the existence of viscous forces. 
This can be done only if there exist streamlines which lie entirely in the 
region of small viscous forces. As before, it follows that such streamlines 
will exist provided u, does not tend to zero, as v-> 0, everywhere in the region 
of small viscous forces. However, whereas the possibility of u, tending to 
zero could be rejected in cases of two-dimensional motion and of rotationally 
symmetric flow without swirl, it cannot immediately be rejected in cases in 
which wu, is not the only component of velocity that may be finite. The 
assumption that wv, does not tend to zero as v—- 0 is here a genuine restriction, 
which places some cases of rotationally symmetric flow with swirl outside 
the scope of the theory. 

We first find that the components of V xw are given by expressions like 
those on the right-hand sides of (4.7), with w, and w, replacing u, and wus. 
All these components are independent of 7, and the integral (2.2) round a 


streamline reduces to an integral round a closed €-line, giving 


' | 
p[(V xo), u, (V eo)» uo] dé 0, (4.14) 


, 1 
When the inviscid relations (4.11) and (4.13) are employed, this condition 


"| ,@H fdC\? dC 1 hu) : ‘a 
2 +C— — —— hk, dé = 0, 4.15 
[ dis? (a) dis hyu, ous viaitiat ( ' 


reduces to 


\ 


which is a more general version of (4.4). 

It is not possible to determine the unknown functions H(i) and C(i) 
from (4.15), and the inviscid flow equations together with the integral con- 
dition (2.2) are not sufficient here to determine the distributions of w and u 
in the region of small viscous forces. However, there is another condition 
that is applicable in this case of rotationally symmetric flow, at any rate 
provided the €-lines are not bounded internally by a solid boundary or a 
singular surface. Since w,—- 0 as v—- 0, the angle between the local compo- 
nents of uandw in the axial plane tends to zero as v—- 0, and the €-lines and 
the vector lines formed from the components of w in an axial plane (the 
latter being referred to as ‘vortex lines in an axial plane’) nearly coincide. 
This approximate coincidence of the two families of curves may take either 
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of two forms. The &-line and ‘vortex line in an axial plane’ that pass 
through any given point may intersect once more at some other point along 
their lengths (or, more generally, may intersect at an odd number of addi- 
tional points) or they may not intersect again (even number of additional 
points.) Inthe former case the ‘ vortex lines in an axial plane’ are necessarily 
closed when v is small, whereas in the latter case they have one end at a point 
exterior to all the €-lines and another end at a point interior to all the €-lines— 
which is possible only if the area intercepted on an axial plane by the region 
of small viscous forces is bounded internally. ‘Thus, provided the €-lines 
have no inner boundary, we conclude that the vortex lines are closed when 
y is small, and hence that, as explained in § 2 


p (V xw).dv = 0. (4.16) 
For small values of v, when w, = 0, this condition becomes 
% ] 
) G XW); Wy (Vv <W)> on | dé = 0. (4.17) 
. WwW 


Combining (4.17) with (4.14), we have 


uty wo\1 Ah 
Us = 0(A,) Ww, W dé wil 
J \uy, ak Op 


which, after some reduction using (4.11) and (4.13), becomes 


-~. dC 

ey iu h,— 7 zr) dé = 0), 

d ¢ c 2 

sal ap r'ujh, ae | = 0. (4.15) 
aus us 


It has already been assumed that the €-lines are not bounded internally, so 
that there will exist an inner limiting ¢ -line which is merely a point. On 
this degenerate €-line, h, = 0 and u,dC/d is finite, in general, so that 
the integral in (4.18) is zero; hence 

dC 

S ( 

qi br? uh, dé = 0. (4.19) 
‘The integral cannot be zero, except perhaps on an isolated €-line (since we 
have already supposed that w, is not zero everywhere, even when v —- ()), and 
(4.19) becomes 


that is 


C = rug = const. (4.20) 
Equation (4.15) now reduces to the form appropriate to flow without swirl, 
and the solution is 


— - = = eon (4.21) 


as before. 
The solution represented by (4.20) and (4.21) is such that the vortex 
lines are circles about the axis of symmetry, as for flow without azimuthal 
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swirl. Moreover, the remarks of the preceding sub-section about the effect 
of viscosity on the rate of change of vorticity at any point, and on the acceler- 
ation of any fluid element, apply here also. When the azimuthal component 
of velocity is allowed to be non-zero, it seems that in truly steady motion 
this component is necessarily an irrotational velocity field corresponding to 
a circulation about the axis (provided, as above, that (a) u, does not tend to 
zero as v —- U0, and (4) there is a degenerate inner €-line). Note, however, 
that the result (4.20) does not apply in the neighbourhood of the axis of 
symmetry, since, in cases in which the region of flow does include part of 
the axis of symmetry, the axis is part of the streamline that bounds the whole 
How (in the axial plane) and this streamline necessarily passes through a 
region in which viscous forces are appreciable. 

In cases in which the area intercepted on an axial plane by the region of 
small viscous forces does have an inner boundary (as, for example, when the 
Huid lies between two anchor rings with a common axis of symmetry, one 
ring enclosing the other), it does not seem possible to deduce the distri- 
butions of w and u in the region of small viscous forces unless the vortex 


lines can first be shown to be closed. 


[1 am grateful to Dr I. Proudman and to Mr W. W. Wood for their 
iseful comments on certain parts of this paper. 
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SUMMARY 

After an extensive survey of mean-velocity profile measurements 
in various two-dimensional incompressible turbulent boundary- 
layer flows, it is proposed to represent the profile by a linear 
combination of two universal functions. One is the well-known 
law of the wall. The other, called the law of the wake, is 
characterized by the profile at a point of separation or reattachment. 
These functions are considered to be established empirically, by 
a study of the mean-velocity profile, without reference to any 
hypothetical mechanism of turbulence. Using the resulting 
complete analytic representation for the mean-velocity field, 
the shearing-stress field for several flows is computed from the 
boundary-layer equations and compared with experimental data. 

The development of a turbulent boundary layer is ultimately 
interpreted in terms of an equivalent wake profile, which supposedly 
represents the large-eddy structure and is a consequence of the 
constraint provided by inertia. This equivalent wake profile is 
modified by the presence of a wall, at which a further constraint is 
provided by viscosity. ‘The wall constraint, although it penetrates 
the entire boundary layer, is manifested chiefly in the sublayer flow 
and in the logarithmic profile near the wall. 

Finally, it is suggested that yawed or three-dimensional flows 
may be usefully represented by the same two universal functions, 
considered as vector rather than scalar quantities. If the wall 
component is defined to be in the direction of the surface shearing 
stress, then the wake component, at least in the few cases studied, 
is found to be very nearly parallel to the gradient of the pressure. 


I. COMPENDIUM 
A. The law of the wall 


Consider a turbulent shear flow which is steady and two-dimensional 
Let u(x, y) and v(x, vy) be the mean velocities in the direction 


on the average. 
Suppose that 


of increasing rectangular coordinates x and y respectively. 
the flow exerts a shearing stress 7,,(~) on a smooth impermeable wall at 


rest at y= 0. For a fluid of constant density, define a friction velocity 


u(x) by 
pus = Ty (1) 
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Experience with turbulent shear flow has shown that, under these conditions 
the mean-velocity profile in a considerable region near the surface is 
described by a relationship called the ‘law of the wall’: 


— = f(%). (2) 
u, Vv 


‘The earliest formulation of the law of the wall was based on observations 
of pipe flow. ‘The current view, however, is one first suggested by Ludwieg 
& Tillmann (1949); the relationship (2) is taken for practical purposes as 
: unique and universal similarity law for every turbulent flow past a smooth 


surface. 
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Figure 1. The law of the wall 


For the special case of steady two-dimensional mean flow of an incom- 
pressible fluid, the form of the universal law is well established. In 
particular, for values of yuw_/v greater than about 50, equation (2) takes the 


form 
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in which « and ¢ are constants to be determined experimentally. On the 
ther hand, the predominance of laminar shear very near the wall requires 
uu, to approach yu, as v approaches zero. ‘To illustrate the presence of 
linear and logarithmic regions in the mean-velocity profile, experimental 
data from several sources are collected in figure 1. <A brief discussion of 
these data will be found in Part II of the present paper. 
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B. The defect law and the equilibrium boundary layer 
The description just given of the mean velocity profile in a turbulent 
shear flow may be summarized in the formula 


De (=) +h(x,¥), (4) 


where the function / is arbitrary except that it is negligibly small in some 
finite region near the wall—say for v/é less than about 0-1, where 4 is the 
thickness of the shear flow. 

For certain special cases frequently encountered (e.g. uniform pipe and 
channel flow and the boundary layer on a flat plate in a uniform stream), 
equation (4) is found experimentally to have the special form 


u (MU) fay ¥ = 
a 1%) +e(165). = 


where II is a parameter which is independent of x and y. Profile similarity 
in terms of the argument y/6 is usually expressed by a relationship known 
as the velocity-defect law, or more properly the momentum-defect law. 
Outside the sublayer, it is an immediate consequence of the logarithmic 
variation of f in equation (5) that 


u,—Uu : \ 
= P(11,%), (6) 
u, 0 
with u = u, at y = 6. 


Quite recently these contributions to the experimental definition of 
turbulent shear flows have been brilliantly extended by F. Clauser (1954), 
who generalized the idea of a defect law by showing experimentally the 
existence of boundary-layer flows for which the similarity laws (5) and (6) 
remain valid although the pressure gradient is positive. 

Finally, the existence of a logarithmic region in the mean-velocity 
profile has been shown by Millikan (1938) to follow directly, without 
reference to any hypothetical mechanism of turbulent mixing, from the 
assumption of the simultaneous validity of the law of the wall in the form 
(2) and the defect law in the form (6). 


C. The law of the wake 

The relationships so far cited, especially the law of the wall, form the 
hard core of empirical knowledge of flow in turbulent boundary layers. 
Once they are accepted, it is natural to look for further inspiration in the 
inverse of Millikan’s problem. ‘That is, given a universal law of the wall 
which is logarithmic outside the sublayer, what other assumptions about 
the motion are necessary in order to establish a defect law? The most 
obvious answer is that the mean-velocity profile must have the form of 
equation (5). ‘The purpose of the present paper, however, is to suggest 
a less superficial answer, one which has several immediately useful con- 
sequences and which at the same time appears to approach somewhat 
closer to the central problem of turbulence itselt. 


F.M. N 
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The essential element in this work is a study of the function h(x, y) 
in the general mean-velocity formula (4). An extensive survey of experi- 
mental data at large Reynolds numbers, reported in Parts II and III, leads 
to the crucial conclusion that this function can be reduced directly to a 
second universal similarity law. Specifically, equation (4) may be written 


in the form 
u yu, I] y 7 
— = — + — / 
u, f( v ) K (3) F (7) 


where II is a profile parameter, as in equation (5), but the function w(y/8) 
is now supposedly common to all two-dimensional turbulent boundary-layer 
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Figure 2. The law of the wake. 


The introduction of a second universal function in the mean-velocity 
profile will be referred to as the wake hypothesis, and the function w(y/5) 
in equation (7) will be referred to as the ‘law of the wake’. Some experimental 
data which support the wake hypothesis are collected in figure 2. The 
measurements in question are included in the survey already mentioned, 
and the development leading to equation (7) and to figure 2 will be found in 
Part III. 

In preparing figure 2 the wake function w(y/5) in (7) has been subjected 


to the normalizing conditions w(0) = 0, w(1) = 2 and | (y/5)dw=1. The 
~ 0 


parameter II is then found to be related to the local friction coefficient 
C, = 2u?/uz by 


5 fi 
7 wIn(™) +04 =, (8) 


u K 


T 


v 
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to the displacement thickness 6* by 








5* u, 
= ( 
K Su. 1+I, (9) 
and to the momentum thickness 6 by 
K? (5* —6\ u2 . ss 


where « and f are constants of order unity. 


D. The equations of mean motion 

To the extent that the similarity laws of the preceding sections are 
empirical, and not based on clear physical principles, these laws cannot 
be extended with confidence to conditions outside the range of observation. 
At the same time, these similarity laws go well beyond the usual limits of 
dimensional analysis. 

In the first place, the thickness 5 and the parameter II are uniquely 
defined by equations (8) and (9), as is easily confirmed by eliminating either 
of these parameters between the equations in question. Inthe second place, 
the mean-velocity field, including the flow in the sublayer, is determined 
whenever the kinematic viscosity v and three of the four parameters u,, u,, 
dand II areknown. Butacomplete analytic description of the mean-velocity 
field implies a complete knowledge (at least within the boundary-layer 
approximation) of the streamline pattern, of the shearing-stress field, and 
of the local rate of transfer of energy from the mean flow to the turbulent 
secondary flow. Some typical calculations based on the formula (7) are 
presented in Part IV. 


E. Physical interpretation 

It is instructive to examine the nature of the wake component in 
equation (7) by letting u, approach zero (i.e. 7,, —> 0) after expressing II 
in terms of 5, 5*, u, and uw, with the aid of equation (9). The result is 


u_ 42 
“ "wer 


which does not involve eitherx or v. It follows, since w(/5) is by hypothesis 
a universal function, that the flow at a point of separation or reattachment 
is locally a pure wake flow. 

This observation suggests a simple physical interpretation for the law 
of the wall and the law of the wake, as illustrated in figure 3. The figure 
shows the profile u(y) for various values of x in a hypothetical boundary 
layer, which is proceeding from separation to separation through a region 
of attached flow. The dashed lines in the figure denote the wake-like 
structure represented in equation (7) by the function w(y/5). The associated 
velocity defect u,—u is given by Ilu,(2—w)/x, and the intercept at y = 0 
of the equivalent wake profile therefore differs from the velocity in the 
external stream by an amount 2I[1u_/«. 

N2 
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According to this interpretation, the law of the wake is to be viewed 
as a manifestation of a large-scale mixing process similar to flow in a wake, 
in that it is constrained primarily by inertia rather than by viscosity. When 
the flow is bounded by a wall, however, it is ultimately necessary to satisfy 
the boundary conditions of vanishing velocity and Newtonian friction at 
the surface. ‘These conditions impose a further viscous constraint on the 
flow whose effect is to modify the mean-velocity distribution as shown by 
the solid lines in figure 3. Near the wall, where the wake mean velocity 
is nearly constant, the constraint provided by viscosity can be observed 
in the sublayer flow and in the similarity relationship known as the law of 


the wall. 
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Figure 3. Mean-velocity field in a hypothet’cal boundary layer, including the 


equivalent wake. 


IJ. ‘THE LAW OF THE WALL 


A. Historical development 

1. Evolution before 1949. ‘The early development of the law of the 
wall, in the hands of Prandtl, von Karmin and others, included a simple 
dimensional argument which has not lost its usefulness. Suppose that 
the mean-velocity profile in a turbulent shear flow is found to be adequately 
represented by a relationship 4(u, y, 4,7,,, “,) = 0, in an obvious notation, 
and that this relationship is found in some region to be independent of the 
characteristic length 6. It follows from the principles of dimensional 
analysis, without any explicit assumptions about the nature of the turbulence, 
that in this region equation (2) must be valid. 

Before the development of the mixing analogy, the function in equation (2) 
was sometimes taken as a power law, for lack of a better representation. 
The sublayer, that is, the region where viscous stress is predominant, was 
treated separately by means of the plausible assumption of a linear velocity 
profile very near the wall. In this approximation 


Cu oy = 2) =F eH tof = u- Vy 


and therefore u/u_ = yu,/v. 
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By 1930 the present general formulation of the law of the wall had 
been achieved, in the sense that no distinction was made in equation (2) 
between the sublayer and the fully turbulent flow. Nikuradse (1930), at 
the suggestion of Prandtl, expressed the law of the wall in the present form 
during an analysis of some measurements in pipe flow. Equation (2) 
also appears in articles by Tollmien and by Schiller in Volume IV of 
Handbuch der Experimentalphysik. 

At about the same time, the mixing analogy of Prandtl (1926) and the 
similarity hypothesis of von Karman (1932) had provided an equation 
Cu(x, y)/ey = u(x)/«ky for the mean velocity in the fully turbulent region, 
with the integral u/u, = (1/«)In[v/yo(x)]+ constant. It is clear, however, 
that the choice of v/u, for the unspecified characteristic length yo(x) is not 
properly a part of the mixing analogy, but rather a part of the dimensional 
argument already mentioned. 


2. Recent developments. Until quite recently, the three most important 
elements in the development of the law of the wall have been, first, the 
dimensional argument leading to equation (2); second, the stipulation that 
the function f is linear at the wall; and third, the recognition that, for 
whatever reason, the function f is very nearly logarithmic in a certain region 
outside the sublayer. 

A fourth important element was added when Ludwieg & ‘Tillmann 
(1949) found experimentally, using an ingenious heat-transfer technique 
for the indirect measurement of surface shearing stress, that for flow in 
a turbulent boundary layer the function f in equation (2) is apparently 
independent of pressure gradient. ‘This result lies at the heart of the present 
study, virtually every part of which depends directly or indirectly on the 
hypothesis that the law of the wall is a unique and universal relationship 
for flow past a smooth surface. 

As several writers in the field have pointed out, equation (2) is an implicit 
equation for uw, (hence for 7,,) when p, », and u(y) are given. ‘The law of 
the wall thus provides a means for accurate determination of the elusive 
wall shearing stress, once the function f in (2) has been established by a 
survey of experimental data on flows for which 7,, is accurately known. 


3. Experimental data. Nikuradse’s classical pipe measurements (1930) 
confirmed the prediction of a logarithmic region in the mean-velocity profile. 
The lowest curve in figure 1 shows the data obtained, including a small 
correction for wall interference, in sixteen surveys at various Reynolds 
numbers of the region between the pipe wall and a value of y/éd of about 
0-15, where 6 is the pipe radius. Nikuradse did not obtain data on the 
sublayer, but inferred the validity of the general law of the wall from the 
argument already given in favour of the limiting form u/u, = yu_/v at 
y =0. This omission was partly repaired by Reichardt’s measurements 
(1940) of mean velocity in the sublayer of a channel flow, and more completely 
repaired by the work of Laufer (1953) and Klebanoff (1954) on pipe flow 


and boundary-layer flow respectively. 
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The latter measurements, shown in figure 1, include data on the sublayer 
obtained in each instance with a hot-wire anemometer. For Laufer’s 
experiments on pipe flow, the wall shearing stress has been derived from 
the observed axial pressure gradient. For Klebanoff’s experiments on 
boundary layers, the wall shearing stress has been derived from the rate of 
momentum loss observed in other experiments with the same model 
(Klebanoff & Diehl 1951; see also Coles 1954). For Reichardt’s experi- 
ments on channel flow, the wall shearing stress was not measured 
independently, as Reichardt’s object was to interpolate experimentally 
in the region of the mean-velocity profile not studied by Nikuradse. 

Finally, some measurements by Sheppard (1947) in a natural wind 
near the ground offer a striking example of a turbulent shear flow for which 
it would seem that no Reynolds number can be explicitly defined. However, 
difficulty in visualizing either an origin of coordinates or a second boundary 
does not prevent the presentation of these data in terms of the law of the wall. 
Sheppard observed the surface shearing stress directly, using the floating- 
element technique, together with the mean velocity at several points up 
to a height of two metres. The result of the measurements is given in 
figure 1. ‘The agreement with wind-tunnel data is surprisingly good when 
it is considered that Sheppard’s measurements were made over a concrete 
surface, and that the vertical temperature gradient may have differed 
significantly from the adiabatic lapse rate associated with neutral stability. 


4. Numerical evaluation. ‘Two empirical constants, « and c, appear 
in equation (3). Throughout the present study, the numerical values 
given to these constants are 

k = 0-40, 

¢ = $1, 
A great variety of other values, especially for «, can be found in the 
experimental literature. However, in practically all cases where equation (3) 
is explicitly taken as a definition, « is found to lie between 0-39 and 0-41. 
Values outside this range are usually the result of operations or assumptions 
which change the definition of x andc. In any event, it is clear in figure 1 
that the data outside the sublayer are well represented by equation (3) 
when « and ¢ are given the values already mentioned. 

Within the sublayer, on the other hand, large fluctuations in velocity 
and cramped quarters for experimentation usually combine to make 
measurements of mean velocity somewhat uncertain. The available 
data in figure 1, therefore, should not be said to establish conclusively the 
uniqueness of the law of the wall in the sublayer, although these data have 
been used elsewhere (Coles 1955) in forming a tentative estimate of the 
function f(yu_/v) and related functions. 


B. Test of the wall law 


1. Momentum balance. Except for a few applications of the floating- 
element technique in flows at constant pressure, values of surface stress in 
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turbulent boundary layers have usually been obtained indirectly from the 
observed pressure gradient and rate of momentum loss, using the momentum- 
integral equation of von Karman (1921). For two-dimensional incom- 
pressible steady mean flow in a boundary layer, this equation is 


- £ uio+ ue, (11) 


where 7,,(x) is the wall shearing stress, u,(x) is the velocity of the free stream 
outside the boundary layer, and 5*(x) and 6(x) are the displacement and 
momentum thicknesses, defined respectively by 


and 


se? u 
= {| —(1-—— 1, 3 
0 fz( “a (13) 


Experience in applying equation (11) to flows with positive pressure 
gradient has usually been that the values obtained for 7,, appear to increase 
rapidly when separation is imminent, contrary to the behaviour which 
is intuitively expected. In fact, the experiments of Ludwieg & Tillmann 
(1949) were originally designed not to study the law of the wall, but rather 
to investigate the validity of the momentum-integral equation (11) by 
providing an independent estimate for 7... 

It is therefore instructive to compare, for various experimental studies 
of flow in turbulent boundary layers, the values of wall shearing stress 
obtained from equation (11) with those inferred from the law of the wall. 
In order to avoid differentiation of measured quantities, equation (11) 
may be integrated with respect to x between two stations x,andx. Replacing 
7 ,, by pu’, the result of the integration can be expressed by the two equations 


uz x §* u? 
(x) = Gy -1+43] Gs a(“), (14) 
Oro “ Lo ) 0 


and 


where uy = U,(X») and 4, = O(xp). 

The dimensionless function ®(x) may be evaluated experimentally 
from equation (14) if 5*, , and uw, are known as functions of x. Quite 
independently, the same function ®(x) may be evaluated from equation (15) 
if the friction velocity u, is obtained by fitting the velocity profile to the 
hypothetical law of the wall. Agreement between the slopes of the two 
functions (x) shows that the flow in question is not inconsistent with the 
logarithmic law of the wall displayed in figure 1. 


2. Presentation of data. ‘The discussion will be limited to data obtained 
at reasonably large Reynolds numbers in experiments at reasonably large 
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scale.* ‘The experimental point of departure is flow at constant pressure, 
for which some mean-velocity measurements of Wieghardt (1943) for a 
free-stream velocity of 33 m/sec are shown in figure 4. Here and in the 
fifteen following figures, the mean-velocity profiles shown are typical of 
the measurements, although some of the profile data have occasionally 
been omitted for reasons of economy in the graphical presentation. 





Figure +. Data of Wieghardt (1943) for constant stream velocity; stations in metres 


All of the profiles, whether shown in the figures or not, have first been 
individually fitted to the logarithmic region of the law of the wall, i.e. to the 
formula u/u, = 5-75 log,)(vu,/v)+5-10, and the resulting values of u,/u, 
have then been smoothed where necessary. In the fitting operation it 
has generally been assumed that the measurements for values of yu_/v less 
than about 200 may be unreliable as a result of large fluctuations in velocity, 
wall interference, poor probe sensitivity at small mean velocities, probe 
position error, or uncertainty in the static pressure. For consistency, 
therefore, the contribution of the sublayer and of the logarithmic region 
to 6* and @ have been computed for the function f of figure 1, i.e. for the 
solid lines in figures 4 to 19, rather than for the actual measurements. 


*This survey of experimental data includes certain material which has not 
previously been reported in detail in the literature. For their courtesy in providing 
this material I am indebted to F. Clauser of the Johns Hopkins University, Baltimore; 
Smith J. DeFrance of the NACA Ames Aeronautical Laboratory, Moffett Field; 
W. Tillmann of the Max-Planck-Institut fiir Strémungsforschung, G6ttingen; 
K. Wieghardt of the Institut fiir Schiffbau, Hamburg; and A. Kuethe of the 
University of Michigan, Ann Arbor. I am also indebted to F. Goddard of the Jet 
Propulsion Laboratory, Pasadena, for making available the services of the JPL 


computing section. 
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Each of figures 4 to 19 includes a sketch showing both the geometry 
of the experiment and the physical extent 5(x) of the shear flow. Finally, 
each of the figures includes a comparison of the two functions ®(x) of 
equations (14) and (15), represented by open points and by a solid line 


respectively. 


Figure 5. Data of Wieghardt (1944) with turbulence grid; stations in metres. 
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Figure 6. Data of Ludwieg and Tillmann (1949) for falling pressure; stations in 
metres. 

3. Unseparated flows. Wieghardt’s measurements (1943) of boundary- 

layer growth for a uniform external stream have already been presented in 

figure 4. Wieghardt (1944) also investigated the effect of increasing the 
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tunnel turbulence level by means of a coarse screen, obtaining the mean- 
velocity profiles shown in figure 5. Using the same channel and instru- 
mentation, Ludwieg & ‘Tillmann (1949) observed a turbulent boundary- 
layer flow in a region of a negative pressure gradient, with the result shown 
in figure 6. ‘These particular profile measurements, and those of figures 13 
and 14 below, are the ones originally cited by Ludwieg & Tillmann in 
their cogent paper on the law of the wall. 
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Figure 9. Data of Bauer (1951) for 60° slope; stations in feet. 
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Some measurements of flow near the lower boundary of a stream 
accelerating under the force of gravity at essentially constant pressure 
have been reported by Bauer (1951). The fluid is water traversing the 
face of a model spillway; data for flow over a smooth plane surface at 
angles of 20°, 40°, and 60° to the horizontal are reproduced in figures 7, 8, 
and 9. In each case the body force per unit volume pg sin f is independent 
of x and y and so plays the same role as a constant pressure force per unit 
volume dp/dx. 
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Figure 10. Data of Schubauer and Klebanoff (1950); stations in feet. 


4. Flows approaching separation. Extensive measurements in a flow 
approaching separation have been made by Schubauer & Klebanoff (1950) 
on a large airfoil at the National Bureau of Standards. ‘Typical mean- 
velocity profiles are shown in figure 10. Upstream of the 7 ft. station, 
the surface shearing stress and free-stream dynamic pressure both increase 








The law of the wake in the turbulent boundary layer 205 


A similar study of a boundary layer near separation is available in some 
work of Newman (1951) for which the mean-velocity data are presented in 
figure 11. Because Newman measured the static-pressure variation within 
the boundary layer, and also made a plausible correction to his observed 
mean-velocity profiles for instrumental errors caused by turbulence, these 
data are among the most accurate and detailed which are available. Newman, 


Figure 13. Data of Ludwieg and Tillmann (1949) for moderately rising pressure; 
stations in metres. 





Figure 14. Data of Ludwieg and Tillmann (1949) for strongly rising pressure; 


stations 1n metres. 
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like Schubauer & Klebanoff, also measured turbulence intensities and 
turbulent shearing stress within the shear flow. 

Boundary-layer flow in diffusers of constant width and rectangular 
section has been studied by Kehl (1943), by Ludwieg & Tillmann (1949), 


and by Clauser (1954), with the results shown in figures 12 to 16. 


In none 


of these five experiments did separation actually occur, and in the last two 


it was deliberately prevented. 
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Figure 15. Data of Clauser (1954) for moderately rising pressure ; 
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Figure 16. Data of Clauser (1954) for strongly rising pressure; stations in inches. 


5. Flows following 


reattachment. 


Turbulent boundary-layer flow 


following reattachment downstream of a separation bubble has been 


investigated by McCullough & Gault (1949). 


Data for the upper surface 


of an NACA 64A006 airfoil section at an angle of attack of 5° are reproduced 
in figure 17. 
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Reattachment of separated flow downstream of a tripping device or 
spoiler, in flow with nominally constant pressure, is illustrated in figures 18 
and 19 with some measurements of Klebanoff & Diehl (1951) and of Tillmann 
(1945). In the experiment of Klebanoff & Diehl transition occurred at 
the spoiler, which was a 1/4-in. diameter rod at the 4 ft. station of the plate. 
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Figure 17. Data of McCullough and Gault (1949) for 5° angle of attack; 
airfoil chord 5 feet. 
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Figure 18. Data of Klebanoff and Diehl (1951) for 1/4-inch rod; stations in fcet. 
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In the experiment of Tillmann the boundary layer was turbulent well 
upstream of the spoiler, which was a rectangular ledge 1-2 cm square at the 
2-02 m station. In both cases it may be noted that the flow far downstream 
has apparently not recovered from the effects of the enforced separation, 


as the mean-velocity profiles do not resemble the profiles for flow at constant 


pressure shown in figure 4. 
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6. Evaluation. ‘The data cited in the previous section obviously do 
not provide a test of the hypothetical law of the wall alone, but provide a 
joint test of the law of the wall together with the turbulent boundary-layer 
approximation and the assumption of two-dimensional mean flow. 

These data confirm that the momentum-integral equation (11) cannot 
be relied upon to give accurate values of surface shearing stress in the 
neighbourhood of separation, for the reason that the left side of (11) is then 
equal to a small difference between two large quantities on the right. 
Consequently, large errors in 7, may be encountered, either as a result 
of inconspicuous departures from two-dimensional mean flow or as a result 
of the omission from the boundary-layer approximation of certain terms 
involving the Reynolds normal stresses and the pressure variation normal 
to the wall. 


Both Ludwieg & Tillmann (1949), using a surface heat-transfer technique, 
and Schubauer & Klebanoff (1950), using extrapolated values of the measured 
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turbulent stress, have found experimentally that the surface shearing stress 
decreases monotonically toward zero on approaching a point of separation. 
Consequently, whether or not the hypothesis of a universal similarity law 
is correct for the flows in question, it is certain from the evidence, for example, 
of figures 10, 13 and 14 that the momentum-integral equation in the form (11) 
is seriously in error. 

Besides the work of Ludwieg & Tillmann, perhaps the most convincing 
evidence for a universal law of the wall is simply that a distinct logarithmic 
region occurs in each of several hundred mean-velocity profiles examined 
here, with very few exceptions; a definite estimate for the wall shearing 
stress is readily obtained, and this estimate is entirely plausible. 

At the same time, it is known (Coles 1955) that a necessary and sufficient 
condition for a universal law of the wall, given the boundary conditions of 
' vanishing velocity and Newtonian friction at the surface, is that the ratio 
u u_is constant on streamlines of the mean flow. ‘The edge of the sublayer, 
as usually defined, is therefore a mean streamline. This result must 
surely be considered in any search for a fundamental order and unity in the 
description of turbulent shear flows. 

It should also be noted that the concept of a universal similarity law 
has recently been reinforced by work of Preston (1954), using an experimental 
technique which depends on the general validity of equation (2) in essentially 
the same way that the use of a Stanton tube depends on the existence of 
a linear profile very near the wall. 

In view of these remarks, the hypothesis of a universal law of the wall 
will be accepted for the purposes of the present paper. In fact, it will 
eventually be suggested that the similarity laws of this and the next section 
may be concepts sufficiently powerful to allow a quantitative treatment not 
only of flows approaching or recovering from separation, but also of yawed 
tlows and flows which are actually separated from the adjacent wall. 


III. ‘THE LAW OF THE WAKE 

A. The defect law 

1. Historical development. A special form of the defect law (6) was 
proposed by Darcy nearly a hundred years ago, and again by Stanton in 
1911, to describe the mean-velocity profile in turbulent pipe flow. The 
defect law in the general form (6) was formulated independently by von 
Karman (1932), who derived an approximate friction law involving two 
empirical constants for the turbulent boundary layer with constant pressure. 
According to experimental evidence from many sources, the defect function 
F(II, v8) in a given flow is insensitive to roughness at the wall, provided 
that the origin for the normal coordinate y is properly chosen. On the 
other hand, it appears from a comparison of figures 4 and 5 that there is a 
small dependence of the defect law on the turbulence level in the external 
stream. 

One well known and useful property of the defect law is that it avoids 
the awkward problem of defining the thickness 6 for a boundary layer. 


O 
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For, if the displacement thickness 6* is computed from equation (12) for 
the profile (6), neglecting the departure of the flow in the sublayer from 
the logarithmic law,* there is obtained 


O* l »\ ’ 
SS ail F(u, 3) a(*) = C(I), 
ou, Jo fo) fa) 


so that 6 is proportional to 6*u,/u,. ‘This implicit notation for the parameter 
II will eventually be replaced by an explicit formula. 





recently been generalized by F. Clauser (1954), who constructed experi- 
mentally two boundary-layer flows with positive pressure gradient, such that 
equations (5) and (6) remained valid. ‘The flows in question have already 
been described in figures 15 and 16. Clauser used the term ‘ equilibrium 
flow’ to denote a flow with a defect law, that is, in the present notation, 
a flow for which the parameter II is constant. 

Examination of figures 7, 8 and 9 suggests that the three spillway flows 
studied by Bauer may be equilibrium flows in the sense of Clauser’s definition. 
So, at least approximately, is the flow with falling pressure studied by 
Ludwieg & Tillmann and reported in figure 6. 


2. The equilibrium boundary layer. ‘The concept of a defect law has 


yu,/8*u, 
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Figure 20. The defect law. 


A comparison of the velocity-defect function F for three equilibrium 
flows can be found in figure 18 of Clauser’s paper, and is repeated here 
in figure 20. ‘This figure, unfortunately, sheds little light on the way in 
which the argument 9/5 and the parameter II are involved in the function 
F(II, y/5) of equation (6), and therefore does not immediately suggest 
any useful generalization of the defect law to non-equilibrium flows. 


* This remark implies that, for example, [fdz is written fz—fzdf, and df is 
replaced by dz/xz as if the function f were logarithmic everywhere. 
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as 


3. The logarithmic region. At this point a further important 
consequence of the defect law in equilibrium flow should be mentioned. 
In the first instance, the mixing analogy of Prandtl or von Karman implies 
a logarithmic variation of the function f(yu_/v) in flows for which 7 does 
not depend on y. However, if the law of the wall is universally valid, 
then the velocity distribution u(y) can be expressed independently of the 
shear distribution 7(y); thus the argument based on the mixing analogy 
in favour of a logarithmic mean-velocity distribution is seriously weakened. 

It is therefore instructive to consider another argument, first proposed 
by Millikan (1938) and based on the wall and defect laws, which also leads 
to a logarithmic function f in equation (2). From the law of the wall, 
u/u, = f(yu,/v), it follows that 


veu vu, f’ yu, 
u,cy yO oe 
From the defect law, (u, —u)/u_ = F(II, y/d), the corresponding expression 
1 : 
Vv Cu vf ¥ 
—~—=--F =}, 
u,cy 3) ) 
where the constant parameter II has been suppressed for an equilibrium 


flow. Now suppose that there is a finite region in which the wall and 
defect laws are simultaneously valid. In this region, the last two equations 


require that 
yu. ,_{ yu, V sol ¥ y eu 1 
LS f eS =-= fk ~ = —— = —— : 
y * v fa) fe) u, oy K(x, V) 


1S 


say. Obviously, «(x,y) is fixed when either of the two variables yu,/v 
or y/6 is specified. But these variables are formally independent of each 
other, since their ratio du,/v may be chosen arbitrarily. It follows that « 
must be a constant. Furthermore, on integrating the expression for 
ou/dy over the region in question, it is found that 


yu, 1, yu, 
f(— } = - In— + const., 
Vv K V 


in agreement with equation (3). 


B. The wake hypothesis 

1. The wake function. ‘The defect law (6) has at least a limited physical 
interpretation, in that the loss of momentum is expressed independently 
of the viscosity. This property, being consistent with the idea of a turbulent 
rather than a viscous transport process, can reasonably be assumed, as in 
Millikan’s argument, to apply everywhere outside the  sublayer. 
Furthermore, the observed sensitivity of the momentum defect to external 
turbulence level and the observed insensitivity to wall roughness are not 
surprising. 

On the other hand, the various mean-velocity profiles so far studied 
are quite systematic in coordinates (i/u,, yu,/v) which involve the viscosity 
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of the fluid. In fact, once the universal law of the wall is accepted, it is 
difficult to escape the conviction that an arbitrarily chosen profile is com- 
pletely determined when the free-stream point (u,/u,,du,/v) is specified. 
To illustrate this remark, figure 21 shows several mean-velocity profiles 
selected from various boundary-layer flows described earlier.* These 
profiles have essentially the same defect law (that is to say, the same value 
of the parameter II) in spite of wide variations in environment. 


U 








Ur | 
y 
0 | 
y | 
4 
fs ee See wee A a, ee 
i 100 1000 100,000 
phic 
Vv 
Figure 21. T'wo-parameter similarity in the mean-velocity profile. 


Taking these properties together, it is something of an anticlimax to 
discover that the puzzle of the defect law is apparently a simple one. The 
key lies not in a study of the defect function F of equation (6), but in a 
study of the original function g(II, v/5) of equation (5), which gives the 
departure of the mean-velocity profile from the logarithmic law of the wall. 
For the three equilibrium flows of figure 20, this departure is shown in 
figure 2, using a linear scale for yu_/5*u,. There is a striking resemblance 
between the three curves, including a common anti-symmetry about a 
midpoint. But this resemblance is obviously not confined to equilibrium 
flows if, as has just been suggested in figure 21, an arbitrary profile in non- 
equilibrium flow coincides in the coordinate system (u/u_, yuo/v) with a 
profile from some member of the one-parameter family of equilibrium flows. 


* Reading from left to right, the profiles are taken from figure 18 (Klebanoff & 
Diehl 1951, station 4.25), figure 16 (Clauser 1954, Series 2, station 108), figure 11 
(Newman 1951, stations D and C), figure 16 (Clauser 1954, Serie 
and figure 14 (Ludwieg & Tillmann 1949, Channel VI b, stations 3. 


2, station 230), 
3 and 3.53). 
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In view of these remarks and especially in view of figure 2, the mean- 
velocity profile may tentatively be written in the form already presented as 


equation (7) 
7 (=) I(x) ( *) 
—_— = f =e a Ww — - 
ut, v K o 


where the function w(y/5), for reasons which will become apparent, will 
be called the law of the wake. If II does not depend on x, then both 
g(II, y/5) in equation (5) and I[w(¥/5) in equation (7) are functions of y/5 
only. This is the property assigned to equilibrium flows by Clauser. 
The present formulation of the mean-velocity profile is, however, more 
general than equation (5), in the sense that the law of the wake, although 
it occurs as a restricted form of (5) which is itself a special form of equation 
(4), is here assumed to apply for non-equilibrium flows. 





2. Normalizing conditions. In order to test the hypothesis of a universal 
wake function in equation (7), it is necessary first to define the thickness 6 
and to specify some normalizing factor for w. ‘To this end, the displacement 
thickness 5* may be computed from the definition, equation (12), for the 
particular profile given by (7). Neglecting the departure of the flow in the 
sublayer from the logarithmic wall law, we obtain 


‘ tr. rs 
é* = 6-(-+-—-—| <=dw), 
u, \K KJg O 


where w, is tentatively defined as the maximum value of w. It is therefore 
convenient to take as a first normalizing condition 


(WU, ay 


5 deo = 1, (16) 





“0 
A second normalizing condition is suggested by the nearly anti-symmetric 
form of the curves in figure 2. The maximum value of w will occur very 
nearly at y/5 = 1, provided that 

ve, = a1) = Z, (17) 


Now w(y/5) is by hypothesis a universal function, so that the boundary- 
layer thickness 6 is uniquely defined in terms of 5* by the integral condition 
(16) and the maximum condition (17). That is, the two relationships 
(8) and (9), together with the identity 


(2) C2)/e@) 


are sufficient to determine all five of the dimensionless parameters u,/w_, 
5*/8, du_/v, d*u,/v, and II (or 5*u,/5u.) when any two are given. For 
example, if the two known quantities are u,/u, and 6*u,/v, as is supposedly 
the case for the data cited in Part II of this study, then equations (8) and (9) 
lead to a simple transcendental equation for II ; 
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3. Profile at separation. Equation (7) as originally written is not in 
a convenient form for use near a point of separation or reattachment. 
However, on multiplication by u,/u, and the elimination of Il by means of 
(9), equation (7) leads to 
u uf ye o* lu.\ a 
cy Cold Oe es ™ Fab 
u, u, \ v ee KU, , 0 


If u_ is put equal to zero in this expression, the result is evidently 


u 5* y ; y 
— — wl —} = dw{ —}). 
uy a) (2) = re) 


Moreover, certain numerical values can be assigned in advance to the ratios 
5*/5, 8/6 and 5*/@ at separation or reattachment. For example, the definition 
of 6 adopted here takes account of the nearly anti-symmetric variation of w 
by requiring that 6*/6 =} at separation. A corollary, anticipating 
numerical values of the next section, is 6/6 = 0-12 (approximately). On 
the other hand, the prediction 6*/6 = 4-2 (approximately) at separation 
or reattachment is a result which is relatively free of preconceptions about 
the form of the wake function w. It should also be noted that this prediction 
is based on examination of the mean-velocity profile in flows which need 
not be close to separation. 


Table 1. The wake function z(€) and related functions. 





Ye w(t) | | Law| ow (1,O w,(II, 2) 
0 0) 0) | 1- 0-06011 1 --0-0001T 0-001 
0-05 0-004 | 0 | 1-0-0021] 
0-10 (0-029 | 0-002 | 1-0-0221 1. 0-02711 +. 0-0001T2 | 
0-15 0-084 | 0-009 | 1-0-0621) | 
0-20 0-168 0-024 | 1—0-119I1 1 0-15411 - 0-008I117 

()-25 (0-272 0-047 | 140-190I1 | 
()-30 0-396 | 0-082 | 1-0-27211 | 1- 0-36011 -0-0441F 
0-35 0-535 | 0-127 | 1+0-36311 | 
0-40 0-685 | 0-183 | 1 0-458I1 1- 0-61411- 00-1271]? | 
0-45 0-838 | 0-248 | 1-0-5521] | 
0-50 | 0-994 0-322 | 1-0-6451] 1- 0-88011 + 0-25711? | 
0-55 | 1-152 | 0-405 | 1-0-7371 | 
(0-60 1-307 | 0-495 | 1-0-8241 | 1+ 1-14311- 0-426112 | 
0-65 | 1°458 | 0-589 1 0-90611 
0-70 | 1-600 | 0-685 1. 0-978I11 | 1- 1-38011- 0-6101F | 
0-75 | 1:729 |. 0°778 1-1-0371] | 
0-80 1-840 | 0-863 1-1-0791] 1+ 1-56111-+0-765I12 | 
| 0-85 1-926 ()-935 1+1-10011 | 
0-90 1-980 | 0-981 | 1-1-0901] 1 - 1-64011 0-82311? | 
0-95 1-999 | 0-999 1+1-051I1 | 
1-00 2-000 | 1-000 1-1-000T11 1- 1-60011— 0-76111* | 





4. Test of the hypothesis. For obvious reasons, the existence and form 
of the hypothetical wake function w(y/5) are most readily investigated in 
flows with a large wake component. Figure 2 shows the function in question, 
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subject to the normalizing conditions (16) and (17), for several mean-velocity 
profiles taken from the present survey. ‘The difference in apparent scatter 
between the non-equilibrium and equilibrium data is the result of applying 


the normalizing conditions to individual profiles in the first case, but to 
For unseparated flows at 


the average of several profiles in the second. 
least, the wake hypothesis appears to be a useful concept, and a tentative 
determination of the wake function w(y/5) is therefore tabulated in table 1 


and plotted in figures 2 and 21. 
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Figure 22. The separating flow of Schubauer and Klebanoff (1950) (measured values 
of shearing stress have been reduced by 31 percent). 


The two-component profile, equation (7), can be made to represent 
quite well the data of Schubauer & Klebanoff, as shown in figure 22; 
of Newman, as shown in figure 23; of Kehl; of Clauser; of Ludwieg & 
Tillmann; of McCullough & Gault at 5° angle of attack, as shown in 
figure 24; and of Tillmann for reattaching flow, as shown in figure 25. 
The solid lines in these figures are computed from equation (7), using values 
for u,, u,, 8, and II which vary smoothly with x. Occasionally it has been 
found desirable to make slight revisions in the original values for the profile 


parameters obtained from consideration of the law of the wall alone. 

A few exceptions to the law of the wake can be found among the numerous 
profiles presented here. For example, it is impossible to find satisfactory 
values for the parameters w,, u., 5, and II in the general formula (7) such 
that the profile at the left in figure 21 can be represented within the apparent 
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Figure 25. The reattaching flow of Tillmann (1945). 
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experimental accuracy. ‘This profile, however, was obtained about two 
boundary-layer thicknesses downstream of reattachment, in a region where 
the real accuracy of the profile measurements is unknown. 

Finally, it may be remarked that Wieghardt’s experiments in figure 5 
indicate that there is a definite change in the shape as well as the amplitude 
of the wake component in flow at constant pressure when the free-stream 
turbulence level is increased. ‘The term ‘universal function’ applied to 
the law of the wake therefore implies that the external turbulence level 
is low, much as the same term applied to the law of the wall implies negligible 
surface roughness. 

IV. Discussion 
A. The equations of mean motion 

1. The turbulent shearing stress. In the past an important aim of 
phenomenological theories of turbulent shear flow has usually been to 
develop a priori some relationship connecting shearing stress and mean 
velocity, by analogy with the relationship 7 = ucu/cy for laminar flow. 
It is now generally accepted that such a relationship probably does not 
exist in any practical sense. However, the same purpose is ultimately 
served by the equations of mean motion, which provide a valid relationship 
between shearing stress and mean velocity. In this respect the present 
analytic representation (7) of the mean-velocity profile may be exploited 
immediately. 

The continuity equation may first be satisfied by introducing a stream 
function u(x, Vv) such that u = c/cy and v = —cys/ex; y& is then constant 
on streamlines of the mean flow. But, if the departure of the profile in 
the sublayer from the logarithmic law of the wall is neglected, equation (7) 
requires that 


Saf udy = e(S  2) eS ae (19) 
, v} . , " ee Ris Gg 6 
A formula for the normal component of velocity may be obtained from 
v = —¢cys/ex, using equation (19) for ¢(x, y), or directly from the continuity 
equation cu/ox+¢czv/cy = 0, using equation (7) for u(x,y). In either case 
it is found that 
ve ydu, yu,wdil 1 d(Iléu.) ¢ “7 ” 1) 
un ude ude tau de Jg3 ™ (2) 
This expression is an exact consequence of equation (7), in the sublayer 
as elsewhere. Its application in the calculation of the shearing-stress 
profile is most easily shown by putting —¢v/éy for cu/éx in the boundary- 
layer momentum equation 


thereby obtaining 
ae _ yam _ "se a(). (21) 
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2. Flow at constant pressure. For flow at constant pressure, i.e. for II 
and u, constant, the expression (21) will be evaluated explicitly for the profile 
given by equation (7). Making the usual approximation in the sublayer, 


the result can be written in the form 


r vy du u \2 u 
2 os ee c— —2i k— lw, + 
: (1 T ) u_ dx | ( = («= Jo 
u u u 
T (« “(x ral — 1) = 2( Ye = W) T 20 | ) (22 


where w, and w, are auxiliary functions which depend on an argument 
y 6 and a parameter II; that is, 


re 


nf... 
w, (II, C) - 1+ — | C' du’, 
S$ JO 
) Mp, M2 
w(,d)=a,++| Olnedw’+—| C(w-—w’)dw’, 
S 0 S =» #9 


where w(¢) is the wake function. ‘The quantities w,(II,¢) and o,(II, ¢) 
are tabulated together with zw in table 1. 

Finally, the derivative du,/dx in equation (22) may be disposed of by 
putting v = 6 to obtain* 


5 du.[T/ u,\? u 
“2 ee («“) 0,-2(« + )Q,+2Q, 
ax tN aU: J ie ’ 


with Q (11) = w, (11,1) = 14+ 11 
and Q.(11) = w,(Tl, 1) = 1+ 1-60011 + 0-761 IT". 


Klebanoff (1954) has recently measured the turbulent shearing-stress 

distribution in a boundary layer with constant pressure such that 
d*u,/v = 9,700, 
6* = 0-400 in. 
The experimental mean-velocity profile for this flow has already been 
shown in figure 1. However, in order to make the present calculation 
independent of Klebanoff’s measurements, it is convenient to estimate 
the parameter II from the data of Wieghardt in figure 4. Using either 
equation (7) and the observed maximum excursion from the logarithmic 
law, or equation (18) and the experimental values for u,/u, and 6*u,/v, 
it is found that II is very nearly 0-55 for flow at constant pressure. ‘Then, 
from equations (9) and (8), with 6*u,/v = 9,700 and II = 0-55, 
du,/v = 2,500, u,/u, = 27-4, 5*/5 = 0-142. 


Furthermore, with 6* = 0-400in., 6 = 2-83 in. 
The profile «/u, computed as a function of yu_/v = (du_/v)(y/5) from 
equation (7) 1s plottedin figure 1. ‘The excellent agreement with Klebanoff’s 


* This expression, like equation (37) of a previous paper (Coles 1954), is a special 


form of the momentum-integral equation (11), and may be used to compute a length 
Reynolds number. The present definition of the thickness 6, however, requires 
the numerical constants 4(1) 7:90, C, = 4:05, and C, = 29-0 of the cited paper to 
be replaced (anticipating the value Il = 0°55) by c+2II/« = 7-85, 2,/« = 3-88, and 
20, x? = 26-4. 
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measurements is to be expected in view of the efforts made at the National 
Bureau of Standards to insure that the flow in question would be typical 
of the fully-developed boundary layer. 

The profile +/q = (7 T )(2u; uz) = —2u'v'/u? computed from equation 
(22) is plotted against y/é in figure 26, together with measured data using 





the present estimate for 6 of 2:83in. Finally, to establish the extent of the 
region of intermittent turbulence in terms of the coordinate y/5 occurring 
in the law of the wake, figure 26 also shows Klebanoff’s measurements of 
the intermittency factor y, defined as the fraction of the time that the flow 
is turbulent. ‘The mean position of the turbulent boundary, i.e. the point 
y = 4, 1s found to be at y/d = 0-825, with a standard deviation of 0-148 
about the mean. 
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Figure 26. Shearing-stress profile in the flow of Klebanoff (1954). 


3. Flow approaching separation. Schubauer & Klebanoff (1950) and 
Newman (1951) have recently carried out hot-wire measurements of turbulent 
shearing-stress in flows approaching separation. ‘The representation of 
the mean-velocity field by the general profile equation (7) has already been 
illustrated in figures 22 and 23 for the flows in question. The corresponding 
shearing stress profiles, computed with the aid of equations (20) and (21), 
are shown in the same figures together with the profiles determined 
experimentally. As noted in the figure, the values of w’v’ reported by 
Schubauer & Klebanoff have been reduced by 31°, in view of the excessively 
large values obtained for 7,, when 7(y) is extrapolated to y = 0. It is likely 
that the hot-wire data are in fact too high, most probably as a result of over- 
compensation. 

Unfortunately, these calculations of shearing stress can be attempted 
only for regions in which the mean-velocity field is in reasonable agreement 
with the momentum-integral condition expressed by equation (11). 
Elsewhere in the flow it is found to be possible to satisfy only one of the two 
boundary conditions 7 = 0 at y = 5 or rt = pu? at y = 0, u, being derived 
from a fit to the law of the wall. The alternative calculations are shown in 
figures 22 and 23 by the dashed lines. Because equation (21) is precisely 
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equivalent to equation (11), the values obtained for 7,, on putting 7 = 0 
at y = 6 are obviously the same as the values implied by the slope of the 
function ®(x«) which is defined by equation (14) and is denoted by the open 
points in figures 10 and 11 for the data of Schubauer & Klebanoff and 
of Newman respectively. ‘The discrepancies encountered in 7,, appear 
to be too large to be caused entirely by the neglect in equation (21) of the 
Reynolds normal stresses and the pressure variation op/éy, although this 
question is still open. 

Another explanation for the failure of the two-dimensional momentum- 
integral equation (11) has been suggested on experimental grounds by 
several writers. ‘The displacement of fluid in a boundary layer is assumed 
to be described by the streamline slope v/u, obtained by integrating the 
continuity equation cu/ox+odv/ey =0. This displacement will certainly 
be affected by lateral convergence or divergence of the general flow, although 
the presence of three-dimensional effects would not necessarily be revealed 
by a study of the mean-velocity profile u(y) measured at various lateral 
stations. 

If z,, is known but the values taken for v/u are even slightly in error, 
the right-hand side of equation (21) will not vanish at the outer edge of the 
boundary layer. Conversely, large errors in 7, may be encountered 


on putting 7 = 0 at y= 6. At one station in each of figures 22 and 23, the 


vu, du “vu YP v 
Vu, du, 

1-— — and — | = dj - 
u> dx So Ce NG 


are plotted separately in order to illustrate the extreme sensitivity of their 
sum 7/7,, to small variations in v/u or in du,/dx. Finally, as a tentative 
correction for three-dimensional flow the streamline slopes v/u have been 
multiplied by suitable constant factors close to unity in order to satisfy 
both boundary conditions on 7. ‘The corrected shearing stress profiles 


are shown by the solid lines in figures 22 and 23. 


quantities 





B. Physical interpetation 

1. Two-dimensional flow. In this paper the term ‘wake function’ 
has consistently been used to denote the function w(y/5) in equation (7). 
The reason for this choice of terminology can be found in some measure- 
ments by Liepmann & Laufer (1947) of a plane half-wake or half-jet; 
that is, in the wedge-shaped region of turbulent mixing between a uniform 
flow and a fluid at rest. 

‘The dimensionless mean-velocity profiles for the fully developed 
half-wake were originally reported by Liepmann & Laufer to satisfy a 
similarity law corresponding to a linear growth of the shear flow. These 
same profiles, after a further translation and change of scale for the 
coordinates,* are compared in figure 2 with the wake function w for flow 


* The present variable € = y/§ is related to the original variable oy/x in figure 13 
of Liepmann & Laufer’s paper by the expression € = 0:505-—+0-33loy/x. When 
€ = 0, oy/x = —1°52; and when € = 1, oy/x = 1:50. Then o8/x = 3-02, and it 
follows that, if o = 12, §/’x = 0:252 rads or 14°4°. 
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in a boundary layer. The residual velocity near y/5 = 0 in the figure is 
the normal component of velocity associated with fluid entrainment, and 
agrees in magnitude with the value v/u, = 0-03 implied at y/5 = 0 by the 
equations of mean motion for the half-wake. 

Although the motion represented by the upper curve in figure 2 is 
bounded on the low-velocity side by fluid at rest, rather than by a solid wall, 
there can be little doubt that the similarity between the various experimental 
curves in the figure is more than accidental, and that substantially the same 
physical phenomenon is involved. The interpretation of the wake 
component in a boundary layer as a large-scale mixing process has already 
been mentioned in Part I in connection with figure 3. 

It goes without saying that if the streamwise mean-velocity distribution 
in a turbulent boundary layer can in fact be expressed as a linear combination 
of wall and wake components, as in equation (7), then so can the normal 
mean-velocity component v, the mean stream function ys, and the displace- 
ment thickness 6*, all of which are obtained from u(x, y) by linear operations. 
This is not so, however, for the flow inclination v/u, the momentum thickness 
4, or the shearing stress 7; and it is most emphatically not so for the turbulent 
fluctuations, except in so far as the wake and wall components might be 
expected to contribute more strongly to the small and large wave number 
regions of the spectrum respectively. On the other hand, a foundation 
has been laid in the present paper for the comparison, at corresponding 
points in various free and bound shear layers, of measurements of inter- 
mittency factor as well as of spectra and intensity of various fluctuating 
quantities. 

Finally, in the event that the parameter II is constant in a turbulent 
boundary layer, a certain balance is implied between the constraints imposed 
by inertia and by viscosity, and between the large-scale and small-scale 
mixing processes. ‘This balance is precisely measured in equation (9) 
by the magnitude of the parameter II or alternatively of the combination 
5*u,/du,. From the point of view adopted in these paragraphs, therefore, 
Clauser’s choice of the term ‘equilibrium flow’ to describe the situation 
when II is constant may well be regarded as inspired. 

It should be noted that Lees & Crocco (1952) have recently attempted 
a qualitative analysis of turbulent shear flows in which they visualize a 
continuous spectrum of mixing processes having both wake-like and 
boundary-layer-like properties. ‘The concept of two flow components, 
one depending on friction and the other on the cumulative effect of pressure 
gradient, has also been independently advanced by Ross & Robertson 
(1951) and by Rotta (1950), both of whom used a term linear in y to represent 
what is called here the wake function w(y/5). ‘These authors did not give 
any interpretation for either of the mean-velocity functions, and were 
therefore able to recommend their formulation of the problem only as 
a useful engineering approach, although Rotta obtained relationships 
which anticipate the present work. As a matter of historical interest, it 
should also be noted that a much earlier attempt by Millikan (1938) to 
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examine the function denoted by A(x, y) in equation (4) was unsuccessful 
because the experimental data studied were not sufficiently precise. ‘These 
efforts have certainly contributed, if not to the specific concept called the 
law of the wake, at least to the atmosphere in which this concept was evolved. 


2. Yawed flow. The present interpretation of the wake and wall 
components in a turbulent boundary layer, as manifestations of the 
constraints provided by inertia and viscosity respectively, is not necessarily 
restricted to two-dimensional flows, provided the two flow components 
are viewed as vector rather than scalar functions of position. Specifically, 
suppose that the general profile equation (7) is rewritten in the form 

q=- 4/74» 


where 


and 





In these expressions, f(yq,/v) and w(y/d) are to be identified with the 
scalar functions previously described for two-dimensional flow, and q, 
is defined as the magnitude of a friction velocity vector q, taken parallel 
to the surface shearing stress T,,, i.e. 

Tw = Pq-4,- 

Furthermore, bothq, and II are assumed to depend on two space coordinates, 
say x and x. Then, if q, and q,, are not parallel vectors, the parameter 
II(x,z) should presumably be interpreted as a linear operator having the 
properties of a square matrix, i.e. as a tensor. It also follows that the 
generalized vector friction law is 

5 2 

qi = q- “() ot ae. 

l K 

This notation is highly tentative, and may have to be revised after more 
data on yawed flows become available. In particular, the present definition 
for qg, is not perfectly consistent with the streamline hypothesis (Coles 1955) 
if t,, is an arbitrary continuous function of two coordinates on the surface. 
However, the notation does allow the interpretation already proposed, 
that the flow near the surface where the wake component is small should 
have the same direction and sense as the surface shearing stress. ‘Thus 
the concept of a constraint provided by friction is a vector concept. 

Given a mean-velocity profile in yawed flow, the vector nature of the 
wall and wake components can be tested in five steps. First, the direction 
of the mean flow near the surface, which is also by assumption the direction 
of the shearing-stress vector t,, and of the wall component q,, is noted. 
Second, the component of mean velocity in this direction is plotted in 
coordinates (q.q,/g=, yg,/v) appropriate to the law of the wall; a fit to the 
function f then yields a value for g,. Third, the thickness 6 is estimated 
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Figure 27. Vector resolution of the yawed flow of Kuethe, McKee, & Curry (1949). 
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from this same plot, for example as twice the value of y for which the profile 
reaches half of its maximum excursion from the law of the wall. Fourth, 
the ‘wall’ vector q, = q, f(5q,/v) is computed and subtracted from the 
free-stream vector q to obtain the ‘wake’ vector q,, = 2IIq,/«. Finally, 
the profile is resolved in oblique coordinates determined by the direction 
of the wall and wake vectors. 

An examination of the experimental data of Gruschwitz (1935) in 
a curved channel suggests that serious errors were introduced in the mean- 
velocity measurements near the surface by the use of a periscope probe. 
In particular, the profiles in the straight portion of the channel differ from 
the concensus of data obtained by other investigators under similar 
conditions. Fortunately, the hot-wire measurements of Kuethe, McKee & 
Curry (1949) on a swept airfoil, although carried out at relatively small 
Reynolds number, involve large angles of yaw within the boundary layer 
and therefore provide a useful test for the concept of vector similarity. 
In treating these data, incidentally, it has been assumed that the indication 
in some of the profiles of a sudden change in flow direction within the 
sublaver is fictitious. 

Ihe airfoil of Kuethe, McKee & Curry was of elliptical planform, 
of 18in. chord and 96-5 in. span, with the major axis swept back at an angle 
of 25. Four profiles obtained near the trailing edge of the airfoil at an angle 
of attack of 14° are plotted in figure 27 in terms of spanwise and chordwise 
components of mean velocity; in terms of streamwise and crossflow 
components; and finally in terms of wall and wake components. 

‘The data in figure 27 can be fairly well represented by the characteristic 
wall and wake functions defined previously for unyawed flows, and it is 
dithcult to say whether or not there is any systematic discrepancy. A more 
stimulating result, which may be coincidental or may illustrate an important 
intrinsic property of strongly wake-like yawed flows, is that the direction 
of the wake component is found in each case to be nearly the same as the 
direction of the gradient of the pressure field over the airfoil surface. ‘Vhat 
is, the final resolution is for practical purposes along the directions detined 
by the two vectors T,, and — grad p, so that (grad p) x (IIt,,) vanishes every- 
where. The inference that the constraint provided by inertia is also a 
vector constraint is potentially useful in investigating the nature of the 
tensor parameter Il when more suitable data become available from 
experiments carried out at larger Reynolds numbers in flow on a larger scale. 
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SUMMARY 

A previous paper drew attention to the collective importance 
tf three physical quantities O, R, S associated with ideal fluid 
flow in a horizontal channel. Invariability of these quantities 
at different cross-sections of the flow implies respectively 
conservation of flow rate, energy and momentum; and their 
values determine a wave-train uniquely. The properties of 
O, R, S are recalled in the present paper to account for the various 
effects of lowering a rigid obstacle into a stream. The conditions 
giving rise to dissimilar types of flow are examined; in particular, 
the circumstances causing stationary waves on the downstream 
side are clearly distinguished from those under which the receding 
stream assumes a uniform ‘ supercritical’ state. A well-known 
result in the theory of the solitary wave is shown to apply to the 
receding stream even when the extreme conditions for the wave 
are exceeded; although it fails to account for a region close to the 
obstacle where the curvature of the streamlines becomes large. 
In passing, a feature of the theory is shown to bear on the practical 
problem of producing a uniform stream. Precise calculations are 
made for the flow under a vertical sluice-gate and under an inclined 
plane. ‘lo account for the region near the bottom edge of the 
sluice-gate, a method based on conformal transformation is used 
whereby an unknown curve in the hodograph plane is approximated 
by an arc of an ellipse. ‘The accuracy of the results is more than 
sufficient for practical purposes, and they compare favourably with 
solutions previously obtained by relaxation methods. A number 
of experiments with water streams are described. 


1. INTRODUCTION 
When the stream in an open horizontal channel is spanned by a fixed 
evlindrical obstacle, the effects of the obstruction may extend to great 
distances both upstream and downstream. Furthermore, each part of the 
stream may take several distinct forms depending on the cross-sections of 
the cylinder, the extent of its immersion, and the nature of the undisturbed 
stream. Kelvin (1886) demonstrated that a train of periodic waves may 
form downstream, whose amplitude increases with the drag on the obstacle 
P2 
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(i.e. with the ‘wave reséstance’). His theoretical treatment was limited to 
cases where the waves are small in both amplitude and length; but Benjamin 
& Lighthill (1954) have shown in general that the wave resistance may 
take any value up to a certain maximum. If the obstacle is lowered 
sufhciently far into the stream, energy losses may occur due to breaking of 
waves, and further immersion may lead to the formation of a fully turbulent 
‘bore or’ ‘hydraulic jump’ behind which waves appear to be absent. 
‘he action of a fixed obstacle can also bring about a transition from a uniform 
subcritical stream (of which the speed is less than , (gdepth)) to a 
supercritical one. ‘This effect is best known through the example of 
a vertical sluice-gate, which is discussed in many text-books on hydraulics. 
‘The conditions under which waves can form ahead or in the rear of a sluice- 
gate have not, up to the present time, been satisfactorily explained. 

The first object of this paper is to present a unified account of these 
various effects, and to make clear the relation amongst them. A general 
method for calculating the form of the receding stream when waves are 
absent will also be outlined, and will be illustrated by applications to the 
tlow under a sluice-gate and under a planing surface. ‘The treatment is 
based on ideal-fluid theory, but a number of experiments have been made 
on water streams in order to test the theoretical results; these experiments 
will be described in the concluding part of the paper. A fact emerging 
incidentally from the discussion appears to have an important bearing on 
the practical problem of producing a uniform horizontal stream of water. 
Such streams are often desired for tests on ship models, and for this reason 
they were made the subject of numerous experiments by Binnie, Davies & 
Orkney (1955). ‘These authors found it possible to produce a satistactory 
subcritical stream without measurable waves, provided the Froude number 
(F speed , (g< depth)) was less than about 0-5. At Froude numbers 
exceeding this value waves always appeared, despite a variety of measures 
aimed at avoiding them. In agreement with this result, the present theory 
indicates that flows converging towards a uniform stream are impossible 
at Froude numbers somewhat less than unity. ‘This fact might be relevant 
to the design of a channel for model tests over a range of speeds. 

In the course of examining the asymptotic properties of a disturbed 
stream, we shall show that a formula known in the theory of the solitary 
wave also applies to flows beyond the extreme conditions for the wave. It 
is subsequently used to calculate the form of the converging stream which 
issues from under a sluice-gate. ‘The flow in the immediate neighbourhood 
of the gate has to be treated in another way; but when suitably combined, 
the two methods account satisfactorily for the whole stream. ‘The results 
of the calculations demonstrate the variation of the contraction ratio with 
Froude number, and hence establish the relation between the sluice-opening 
and the total discharge. ‘They are apparently more accurate than previous 
estimates, notably those due to Pajer (1937), and the single result which 
Southwell & Vaisey (1946) calculated by relaxation methods. ‘The latter 
authors drew attention to the difficulties which confront their methods 
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when applied to flow near points of separation from solid boundaries 
(p. 160 in their paper). ‘These difficulties appear to be responsible for 
most of the discrepancy (about 1-5",,) between their computed value of the 
sluice-opening and the value according to the present theory. From 
a practical point of view there is admittedly little justification for carrying 
the calculations to high orders of approximation, since the theory does not 
apply with precision to real fluids. In the experiments described later, 
for instance, the thickness of the boundary layer on the channel bottom was 
observed to be about one-fifteenth of the sluice-opening, thus indicating 
that frictional effects were far more significant than the final error in the 
ideal-fluid calculations. An accurate form of the theory is nevertheless 
desirable, for at least it gives assurance that the discrepancies observed in 
real flows are due to friction alone. 

‘lhe problem to be considered is indicated in figure 1. We shall assume 
that the motion is steady, and that both far upstream and far downstream 
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Figure 1. Side elevation of sluice-gate. 


the flow is uniform and horizontal, having depth and velocity hy, uy and A, u 
respectively. In the paper already cited, Southwell & Vaisey noted the 
existence of a symmetrical solution in which the free surface regains its 
original level after passing the sluice; but flow of this kind does not occur 
in practice unless there is some severe obstruction downstream. The 
free surface usually falls to a lower level, and Binnie (1952) has shown that 
the approaching stream must then be subcritical (F < 1) and the receding 
stream supercritical. ‘The condition of continuity is 

gts = uh = OQ, (1.1) 
where O is the discharge per unit span; and Bernoulli's equation may be 


written 


gh, + hu; = gh+3u? = gH, (1.2) 
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where H is the total head measured above the bottom of the channel. If 
both O and H are given, Ao, A, uy and u can be calculated from (1.1) and 
(1.2), but there is no simple way of finding the sluice-opening s. 
Alternatively, if s is specified together with either O or H, as is usual in 
a practical example, elementary methods fail to give any of the other 
quantities. In addition to (1.1) and (1.2) a third simple relation may be 
written down equating the reduction in flux of momentum to the resultant 
of horizontal forces, which include the force P on unit span of the gate and 
the pressure forces in the oncoming and receding streams; thus, denoting 
density by p, we have 


p(u-h, — u7h) P — | po(h=—h?). (2.5) 


However, since P is unknown, this equation is of no assistance in the main 
problem. 

At point 4 in figure | the free surface rises to stagnation level, and the 
velocity there is zero. At point B it must slope vertically, since otherwise 
the bounding streamline would make a sharp corner there and the velocity 
would become infinite, which is clearly impossible if the surface is tree. 
One may nevertheless expect that the free surface has infinite curvature at 5, 
for this property has been established in many similar instances when 
gravity is absent (see, for example, Southwell & Vaisey, p. 160). 


2. GENERAL THEORY OF THE FLOW 

The investigation is conveniently begun by considering the ditterent 
types of flow which may exist at points far removed from the obstacle. 
In other words, we shall first examine the possible asymptotic forms of the 
How upstream and downstream. At the outset it is necessary to emphasize 
that the stream cannot be perfect/y uniform at any finite distance from the 
obstacle. ‘This fact becomes obvious when one considers that the hodograph 
variable ¢ = ge’ (where g and # are the magnitude and direction of the 
velocity) is analytic within the stream; therefore, by virtue of an elementary 
property of analytic functions, ¢ cannot be constant throughout any finite 
region unless it is constant everywhere. Having thus dismissed the case 
of a uniform stream, there remain only two types of flow capable of extension 
to indefinitely large distances upstream or downstream. ‘The tirst possibility 
is a train of periodic waves, and the second a steadily converging stream 
resembling the outskirts of a solitary wave. In practice, flow of the latter 
sort may be indistinguishable from a uniform stream except in a region 


close to the obstacle; strictly speaking, however, they become uniform only 
at infinity. Although this distinction is a trivial one in some respects, 
it will assume importance in the subsequent argument. 

It is now desirable to recall some of the ideas put forward by Benjamin & 
Lighthill (1954) in connection with the theory of gravity waves of finite 
amplitude. ‘They pointed out that for any two-dimensional steady flow 
in a horizontal channel, there are three important physical quantities which 
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in the absence of friction and horizontal external forces, have the same value 
at every cross-section of the flow. These are O the volume flow per unit 
span, R the energy per unit mass (i.e. g times the total head as measured 
above the channel bottom), and S the resultant of momentum flux and 
pressure force per unit span divided by the density. A train of long waves 
was shown to be determined uniquety by the values of O, R and S, and 
there is reason for believing that this property is common to all wave-trains 
in parts distant from their origin and termination. ‘This view is supported 
by the recent work of De (1955), who calculated numerical values of the 
quantities in question for a wide range of wavelengths. ‘The former authors 
also demonstrated that the physically realizable combinations of O, Rand S 
are confined within certain limits. For instance, if O and R are fixed, 
then S is restricted between a higher value corresponding to uniform 
subcritical flow and a lower value corresponding to uniform supercritical 
flow. Stationary wave-trains may occur for intermediate values of S; 
but the wave amplitude tends to zero as the upper limit is approached, 
whereas at the lower limit the wavelength becomes infinite, and the only 
possible wave is the solitary wave. 

We proceed by extending a line of argument begun in the paper by 
Benjamin & Lighthill (p. +55). If a rigid obstacle is lowered gradually 
into a slightly subcritical stream, it experiences an increasing wave resistance 
{corresponding to a reduction of S on the downstream side) until the waves 
formed downstream become of great length, and their profile approaches 
that of the solitary wave. A different situation arises, however, if the 
Froude number of the respective supercritical flow (i.e. the flow with the 
same values of O and R as the subcritical flow upstream) is greater than 
the value 1-25, at which, as McCowan (1894) showed, the solitary wave 
takes its extreme sharp-crested form. A gradual approach to the lower 
value of S is now impossible, since the extreme condition for the solitary 
wave is exceeded at the limit. In fact, as S is reduced from its upper limit, 
periodic ‘ waves of maximum height’ occur at a value of S in excess of the 
lower limit, and a further gradual reduction will result in energy losses due 
to breaking of waves. The obstacle must then be lowered by a finite 
amount, bringing S to the minimum value, before steady lossless flow is 
again possible. No wave can form on the supercritical stream under these 
circumstances; but it will be shown presently that the stream has some 
features in common with the solitary wave. Recognition of the lower 
limit to S also shows that any obstruction of a stream initially in a supercritical 
condition is bound to result in energy losses, since a reduction of S is other- 
wise impossible. In fact a bore will be formed upstream. A bore 
sometimes occurs in practice on the stream receding from a sluice-gate. 
It may be caused by some obstruction, or simply by friction if the channel 
is fairly long. The flow assumes a subcritical state behind the bore, but 
with a smaller value of R than the subcritical flow upstream from the gate. 
( he material of this paragraph may be made somewhat clearer by consulting 
figure 1 in the paper by De (1955).) 
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The condition F > 1-25 noted above may be shown, by means of (1.1) 
and (1.2), to correspond to the condition Fy 0-792 for the respective 
subcritical stream. Hence we can conclude that if F, 0-792 a transition 


from subcritical to supercritical uniform flow, in the manner indicated by 
figure 1 (present paper), cannot be achieved with R unchanged, since only 
Wave-trains can result on the downstream side. Although such a transition 
is possible if F,, -- 0-792, the obstacle must be immersed to a considerably 
greater extent than that which causes waves. ‘Thus a clear distinction is 
made between the action of a sluice-gate and that of a slightly submerged 
obstacle causing a wave-train downstream. ‘These considerations seem 
to suppiv a complete answer to the question whether waves can form on the 
downstream side of a sluice-gate. [his matter was raised by Southwell & 
Vaisey (1946), who were led to expect waves in the course of their 
calculations, although none were found. it was also considered by Binnie 
(1952), who presented a tentative argument indicating that waves are 
impossible. ‘The question whether waves can form on the upstream side 
is deferred to a later part of the discussion. 

The special case of critical flow (F 1) deserves separate attention. 
Let us therefore consider the flow in a channel which is supplied from 
a large reservoir, so that R is fixed, and which is initially free from obstruction 
of any sort. Such circumstances are well-known to give rise to critical 
flow, which makes the value of O the maximum possible at the given total 
head. ‘his fact may be established by various methods, a review of which 
was given by Binnie (1949). Suppose now that an obstacle is fixed in slight 
contact with the free surface at a point well down the channel, thereby 
causing a small reduction in S on the downstream side. ‘The disturbed 
flow is clearly unstable; for the change in S necessitates a decrease in Q 
from its maximum value, since R cannot change, and hence an increase 
in depth on the upstream side. As the free surface rises upstream, a greater 
force is exerted against the obstacle, and the downstream value of S further 
decreases. Before steady conditions can be resumed, the free surface 
must rise to stagnation level on the front of the obstacle. ‘he elevation 
of stagnation level above the channel bottom (i.e. the total head H) 1s easily 
shown to be 1-5 times the depth of critical flow. ‘Thus the obstacle is 
‘ wetted ’ toa height }H in the final steady state. ‘The resistance experienced 
by an obstacle this far immersed is much too large to be ascribable to wave 
formation; and the flow is therefore of the kind illustrated in figure 1. 

It is now proposed to examine the profile of the receding stream in the 
case where waves are absent. ‘The work of Benjamin & Lighthill again 
provides a convenient starting point. In a new presentation of ‘ cnoidal 
wave ’ theory, they demonstrated that the condition of constant O and S$ 
leads to an approximate differential equation for the free surface (the condition 
of constant R being superfluous). In terms of coordinates (x, v) with the 
y-axis along the channel bottom, this may be written 


iy\2 
10% =) + evi 2Rv242Sv—O? = 0, (2.1) 


a! - 
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The nature of the approximation on which (2.1) is based is noteworthy. 
The assumption was made that successive derivatives of the slope v’ diminish 
fairly rapidly in order of magnitude, so that, for instance, y?v'‘ can reasonably 
be neglected in comparison with y”. The next stage of approximation 
was observed to contribute terms of the order of v4 to (2.1), and further 
stages would evidently contribute terms of the order of successively higher 
powers of y’*. The trend of the coefficients in the early stages suggests 
that the method of approximation is rapidly convergent, provided v’ is 
fairly small. There is accordingly some justification for assuming that 
an exact solution exists for any flow with constant O, R and S, and that 
(2.1) will provide a close approximation to it in parts where \’ is small. 
As will be explained in the next paragraph, this equation can lead to an 
approximation for the solitary wave. The mathematical existence of this 
wave has been proved by Friedrichs & Hyers (1954), and so the validity of 
(2.1) is confirmed in this one respect. Its validity in the present application, 
however, cannot be established by appeal to the work of Friedrichs & Hyers ; 
for the solitary wave is impossible with the values of O, R and S to be con- 
sidered here. Nevertheless, it may be said that the phy sical aspects of the 
present problem, in which the flow-pattern is supposed to be stationary 
relative to the obstacle, are a good deal clearer than these aspects of the 
solitary wave, which is stationary only with respect to a hypothetical frame 
of reference. 
If O, Rand S take the values 

O = uh, R = gh+4u’, S = hgh? + wh, (2.2) 
which are appropriate to a uniform stream with depth / and velocity u, 
and whose Froude number is therefore F = u/y (gh), the cubic expression 
in (2.1) factorizes and the whole equation can be arranged in the form 


(Hy RGay(my). es 


With F > 1 this is identical with the equation for the solitary wave obtained 
by Rayleigh (1876a). Other than the uniform flow y = 4, its only solution 


i 


2 


2 kis gM (2.4) 


j 1 +(F?—-1)sech? | | 


2Fh 
The solitary wave is thus the only wave which can emerge from a uniform 
supercritical stream without loss of energy or momentum. 

This approximation to the solitary wave is most accurate at the outskirts 
of the wave where the slope and curvature are small. When the wave 
approaches its extreme form (at F = 1-25), the approximation is found to 
be no longer accurate over the whole wave, vet still holds at the outskirts. 
It is clear also that these results apply to the outlying parts of a stream with 
constant O, R and S which converges towards uniform flow at F > 1-25, 
although the solitary wave is then impossible. ‘The streams of concern 


here are clearly included in this category. Accordingly, their free surfaces 
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must satisfy (2.4) approximately up to a region close to the obstacle where 
the slope or curvature becomes large. ‘There is a simple way of testing this 
assertion; for, as Addison (1938) and Binnie (1952) have pointed out, 
the flow under a sluice-gate when the ratio s/H is small resembles a jet 
issuing under pressure in the absence of gravity. ‘The latter example is 
illustrated in figure 2 (a), where the bottom line may be taken to represent 
either a plane wall or the centre-line of asymmetrical jet. Rayleigh (18764) 
found the complete solution of this problem. ‘The free surface satisfies 
the equation 


y os i 
— = 1- -sin#, (2.5) 
d 1 

where # is the angular slope tan '(dy/dx) which ranges between 0 and — }7, 


and d has the same meaning as / in figure 1. In the corresponding gravity 
problem F is very large, and (2.3) consequently gives 
y ] 

SS ee (2.6) 
h V3 

Agreement between the two results is partly attributable to the similarity 
of the numerical factors 1 3 and 2/7 appearing in (2.5) and (2.6), but is 
in fact better than this alone would suggest. “The Cartesian equations 
equivalent to (2.5) and (2.6) are easily found; hence the agreement can 
be checked very readily by plotting the two curves on tracing paper and 


fitting them together. ‘hey are found to fit closely in parts where the 
magnitude of @ is less than about 30, the maximum discrepancy in y/h 
being about 1". ‘he curves diverge significantly in the higher region where 


the curvature of the first of them (from (2.5)) becomes large. Note that 
tor smaller Froude numbers the accuracy of (2.4)is likely to improve, whereas 
(2.5) becomes inaccurate. 

The approximation of (2.4) may be tested in one further instance by 
comparison with the solution for a sluice-gate which Southwell & Vaisey 
(1946) obtained by relaxation methods. In their example the assumed 
conditions were, in the notation of tigure 1, H = 12andh, = 11 length units. 
‘The depth and Froude number of the flow far downstream were not given 
explicitly, but from (1.1) and (1.2) they are easily calculated to be h = 3-854 
length units and F = 2-055. When these values are put into (2.4), the 
curve obtained from this equation is found to coincide with the lower 
part of the profile given by Southwell & Vaisey (p. 151 in their paper), 
but begins to diverge perceptibly from it at a point where the slope is about 
— 30°, a little way below the highest point. 

We turn from the main argument for a moment in order to consider 
another matter explained by the differential equation (2.3). If F < 1 the 
equation has no real solution except y = h, the case of uniform flow. This 
fact demonstrates convincingly that non-uniform flow converging towards 
a uniform stream, as discussed previously, is impossible for values of F less 
than unity yet large enough for a long-wave approximation to apply. Any 


ittempt to produce a uniform stream under these conditions is therefore 
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bound to fail, unless the flow can be made uniform on entry into the channel. 
Clearly, the result of any non-uniformity on entry will be a train of periodic 
waves, the values of O, R and S being differently correlated than for uniform 
How. The occurrence of waves as predicted here was observed experi- 
mentally by Binnie, Davies & Orkney (1955), who mentioned similar 
findings by other investigators. 

This is a suitable place to raise the question of wave formation on the 
upstream side of an immersed obstacle. ‘The situation at small Froude 
numbers, to which a short wavelength approximation can be applied, was 
discussed at length by Lamb (1932, §§ 242-5), whose treatment was restricted, 
however, to infinitesimal disturbances of the free surface. He noted that 
in the absence of dissipative forces the problem is te some extent indeter- 
minate, since waves of small amplitude can always be superimposed on the 
stream without altering certain of the physical conditions (e.g. QO or S). 
There appears to be no simple way of extending Lamb’s argument to account 
for waves of finite size: but the presence or absence of waves seems likely 
to depend only on conditions specified far upstream. Hence one is 
justified in assuming tnat the flow is uniform upstream, as is done in the 
following treatment of the sluice-gate problem. 


3. FLOW UNDER A SLUICE-GATE 

Equation (2.4) has been shown to apply behind, but not too near, any 
bstacle adequately immersed. Other methods are needed, however, 
to deal with the rapidly curving portion of the free surface which is always 
to be expected near a‘ trailing edge ’, and which clearly depends on the shape 
f the obstacle. A sluice-gate and planing surface are the only examples 
treated in this paper, but similar calculations could undoubtedly be carried 
out for other boundary shapes. 

It is reasonable to suppose that in the region of B in figure 1 the flow is 
little atfected by gravity, since the acceleration of the fluid particles is large 
there. Although (2.5) will roughly give this part of the free surface, the 
solution of the problem indicated in figure 2 (4) is rather more accurate.' 
Here the flow upstream from the jet is confined within parallel walls at 
a distance D from the centre-line; and, as a basis for comparison with the 
example in figure 1, we may take D = Hasa useful estimate. This problem 
was solved by Mises (1917), although a formula for the free surface was 
not given in his paper. Straightforward methods show, however, that 


the free surface satisfies the equation 


2d F 
y = d— — tan (m sin 9), (3:5) 
. TH 
2Dd , 3 ' 2tan-'m (3.2) 
with m= ——, and — = 1+——— 3.2 
D? — d?” d =m 


This result can be combined with (2.3) to obtain a rough estimate of the 
complete free surface below the sluice. “The two curves may be joined 
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1 


at a point where 6 = —30,, sav, so that in (2.3) y? = (tan30°)? = 4. To 
calculate the sluice-opening s for a given value of F, which specifies an 
appropriate value of D and also the height of the ‘joining point’ (19, say), 
the two equations (3.2) are used to eliminate k and d from the relation 


2d 


S—Vp) = (tan-! m—tan-! 3m), (3.3) 
; ss 2 


7 
which follows from (3.1). In the example treated by Southwell & Vaisey 
this method leads to a value 0-602 for the contraction ratio h’s, whereas they 








Figure 2. Free streamline problems not affected by gravity 


obtained 0-608.* ‘The method takes no account of velocity changes over 
the upper portion of the free surface, although over all parts the resultant 


velocity g varies according to Bernoulli’s equation 

g" 29(H YY). (3.4) 
As the free surface falls between heights s and yy, the velocity increases by 
a factor , (f/—s)/, (H—.,), which is found to be significantly greater: 


* Their result was misprinted as 0-66; but Miss Vaisev has confirmed that 


the intented value is 0-608 
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than unity in typical examples. We are therefore led to seek a more accurate 
method taking full account of gravity. 

By the method now to be described, a solution is found which satisfies 
the boundary condition (3.4) exactly at the edge B and at a ‘joining point’ 
below which (2.3) can be applied accurately. ‘The method is developed 
from one given by Pajer (1937), and so the present account is made fairly 
brief. Following the usual procedure in dealing with irrotational motions 
in the plane of a complex variable z, use is made of the fact that since the 
velocity potential ¢ and stream function y% are harmonic, the variable 
w = 6+ isa function of zonly.. The first aim of the method is to establish 
an approximate relation between w and the hodograph variable 


¢ = dwdz = qe’. We consider a transformation to the plane of an 
auxilliary variable f(z) defined by 
ail a : 
P a ei =, (3.5) 
¢ q . s 


where p and & are positive real constants. ‘The separation of the real and 
imaginary parts of (3.5) leads to 





cos 4 ; k . 
P -_ = (Ir1- a) cos w, (3.6) 
q ~ 
sin : k ed 
sted = (1 | + A) sin w, (3-7) 
q If] 
vhere w = arg(f). We also make use of the transformation 
f’= fr +1/f?. (3.8) 
Itisnow assumed that w = 6 = — }zatthe edge B,and that |f| = 1 along 


the free surface below this point. According to (3.5), the circle so defined in 
the f-plane is mapped in the hodograph as an ellipse, whose size and eccen- 
tricity depends on p and k. By a suitable choice of these constants, the 
points in the hodograph representing B and the ‘joining point’ may be 
linked by an arc of the ellipse. 


From (3.8) we have that f’ = 2 cos 2w on the free surface below B. At B, 
f’ 2,sincew = —47;andatinfinity downstream f’ = 2, since w = 4 = 0. 
In other words, the free surface is mapped in the f’-plane between — 2 and 2 
along the real axis. On the straight line 4B in figure 1, w = 6 = —!7; 


and g = 0 at A. Thus AB is mapped between infinity and —2 along the 
negative real axis of f’. If the free surface upstream from 4 is taken to be 
horizontal, we have w = @ = 0, and so f’ is real along this surface also. It is 
therefore mapped along the positive real axis between infinity and a value 
f, = f7+1/f2 representing the conditions far upstream. It is a simple 
step to show that the lower half of the f’-plane maps the entire flow-pattern. 

The streamline 4 = 0 extends along the bed of the channel, and the 
streamline ¢ = O(where O is the volume flow rate) forms the upper boundary 
of the flow. As shown above, the latter streamline is mapped along the 
real axis in the f’-plane. The figure in the w-plane is an infinite strip of 
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width O lying on the real (@) axis, and is mapped in the lower half of the 
f'-plane by the Schwarz—Christottel transformation 

: O f’—fe\ > 

w= +m = a loe( 5) , (3.9) 

‘The relation between ¢ and w is determined by (3.5), (3.8) and (3.9) 

The form of the free surface below B can hence be deduced by a standard 

), details of which may be omitted here. We 





method (see Lamb 1932, $7 
tind for this surface 
- a ae 
eS WwW 77 COS Ww 
4 ore He nha 3 
deal | logtan( S 5) ntanh ( - )| (3.10) 
Pog 1b 
CS x m 
\ s— — —— [tan'm+tan-(msinw)], (3.41) 
: mz 1—k 


where n = }(fy+1/fy), m = 2fo (f5—1), and c is the value of y/s for w = 0. 
If the boundary condition (3.4) is to be satisfied at B, we must have, 

by reason of (3.7), 

ai te 


If this condition is also to be satistied at a point below B where v = 1, 


wy), We must also have 


sin 4, 
P : (1 t k)sin wy. 


\ (2e[H — v]}) 
H — vo sin 4, \° - 
AY (= (3.12) 


and Ww 


It follows that 


H- s SIN Wo 

in which w, is determined (from (3.6) and (3.7)) by 
l—- en 
tan (). (3.13) 


tan w, rare 

l+k 

There will be no need to consider further the constant p, which merely 
tixes the scale of the figure in the f-plane. 


introduce a new constant « = tan-'!m (hence 


It is convenient to 
f,), in terms of which (3.11) gives 


2¢ s ] t k 
ee tak 
(in (3.11), that 


Cot 1% 


(3.14) 


cota{x+tan (tan xsinw,)]. 


We also tind, by putting & 
l 21+k eee 
r , (3.15) 





Fs 
Furthermore, if /, is the value of the depth upstream consistent with the 


assumption of its being constant everywhere, the continuity condition (1.2) 


leads to 
h, ' 
(fo —Rk) (1—k) = (cot }a—Atan $x)/(1—R), 


from which, if a substitution for ¢ is made from (3.14), we obtain 


h 2h h 2h : 
h( — - = % cot x—tan ix) = s+ — z~cotx%—cot $x. (3.16) 


s 


Ss 
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‘These results represent a flow which closely resembles the true solution 
up to the cross-section where y = y); beyond this the methods described 
previously can be applied. Our main object is to calculate the contraction 
ratio fs as a function of the ratio s/H, which defines the flow-pattern uniquely 
on the basis of Froude number scaling. It is first necessary to obtain a 
suitable estimate of the quantity A, appearing in (3.16). We have so far 
neglected the small difference between H and h, (see figure 1), but now 
tentatively take account of it by assigning to 4, a value somewhere between 
the two. (Note that no continuity condition is violated by doing this, since 
we are merely seeking a flow particularly resembling the actual flow in the 
region just below the sluice.) Inspection of the flow-pattern obtained by 
Southwell & Vaisey (1946, figure 25) suggests that }(H+2h,) is a fair 
estimate for /,, although the value of A, is not at all critical when s/H is 
fairly small. An initial estimate of A/s is required to calculate Ay from the 
relation 

H(h,—h) = h-+hh+h?, (3.17) 
which follows from (1.2) and (1.3). When 4, is specified, vgs can be found 
by eliminating c, k, x and w, amongst the five equations (3.12) to (3.16). 
The calculation has to be done by numerical methods, but the following 
procedure was found to give fairly quick results. The value of « is first 
guessed. ‘This task is aided by extrapolation after solutions have been 
obtained for some values of s/H, and as a help in this respect we have 
z= for sH=0. Using this estimate of x, k is calculated directly from 
(3.16) ; then ¢ is calculated from (3.15), and w» from (3.13). Finally, vp» 
is obtained both from (3.12) and from (3.14). This sequence is repeated 
with another estimate of «, and it is then possible to assess « more accurately 
by interpolation designed to bring the alternative estimates of yg together. 
Repetitions of the whole process with successively closer approximations 
to x lead rapidly to an accurate estimate of y). ‘The method is helped by 
the fact that the value of y) given by (3.12) varies greatly with x, whereas the 
value from (3.14) varies very little with «. 





Su x | 3 q\ Ye | 4-4 
(ft) | (ft) | (f/s) | (f/s) | (f/s) 
909 | O 1-0 | 16:048 | 16-048 | 0 
75 | 0-0024 | 0-9863 | 16-085 | 16-075 | 0-010 
60 | 0-0197 0-9464 16°186 | 16-155 | 0-031 
45 | 0-0692 | 0-8835 | 16-318 | 16-280 0-038 
35. | -0-1324 | 0-8309 | 16-417 16-383 | 0-034 | 
25 | 0-2368 | 0-7719 | 16-499 | 16-499 | 0 
Table 1 
Forexample, with s/H = 0-2 and 4, = —25° itis found that y/s = 0-7719. 
Coordinates of the free surface in the region where — 90 6 < —25° are 


given in table 1, which also includes values of g given by (3.6) (q,), and the 
correct values according to (3.4); s is taken as 1 ft, and g as 32-19 ft/sec?. 
The largest error in q is seen to be slightly greater than 0-2°,,. 
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After the value of y, has been established in this way, the calculation 
is completed as follows. If (2.3) is applied at the ‘joining point’, it gives 


3 lV 2 = MM dy 2 
mG 1) G i) (5) ewe 


From (1.3) we also have 
H 
h 
Che elimination of F? between these equations leads to 
(Vp —A)?(2H — 2h — yo) — 3(H — h)h? tan*6, = 0. (3.18) 
When divided by s*, this equation becomes a cubic in A s whose coethcients 
are functions of the known quantities H/s, vg/s and 6,; this can be solved 
in a straightforward manner. 
Values of A,s calculated with 6, = —25° are tabulated as a function of 
s H in table 2, which includes the Froude number according to (3.13). 
The table also shows the set of values obtained by Pajer (1937), who used 





Fes h's hs F 
| (Pajer) 
0 0-6110 0-6110 x 
()-] 00-6060 | 5°56 
()-2 0-6022 06046 3-822 
0-3 0-5995 0-6036 3-020 
()-4 05980 00-6043 2-522 
(-5 0-598] 0-6066 2:165 
Table 2 


1 method similar to the first part of the present method but extended over 
the whole free surface. ‘he errors in velocity were then considerably 
larger than at present, and were in fact largest in a region corresponding 
to# = —25 , where, by the present method, the error is deliberately reduced 
to zero. Moreover, no account was taken of the ditference between H and 
h,. ‘The discrepancies between Pajer’s results and the present ones are 
fairly small, although quite significant. A remarkable feature of table 2 
is the smallness of the variation in fs. ‘This effect is mainly due to gravity ; 
for in its absence (the case illustrated in figure 2 (4)) the contraction ratio 
increases quite rapidly as the flow is narrowed upstream, and is always 
greater than the extreme value for the case in figure 2 (a). 
The present method requires slight modification when applied to the 
problem treated by Southwell & Vaisey, in which H, / and F are given and 
s remains to be found. First, an estimate of s'H is made, and then /s is 


calculated as before. Since hs varies quite slowly with s H, this result 
leads to a much closer estimate and the approximation is then further 
improved by repetitions of the process. Starting with h/H = 0-3212 and 
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I? = 4-227, which are the values appropriate to Southwell & Vaisey’s 
example, the method gives finally A's = 0-5991. We recall that the value 
obtained by them was 0-608. ‘The curves derived from the present theory 
and by Southwell & Vaisey are both shown in figure 3. »The two are 
indistinguishable up to a region just below B, but diverge appreciably in this 
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region. ‘The form of Southwell & Vaisey’s curve very close to B seems 
unlikely to be correct, since it requires the streamline to undergo an unduly 
sudden change of direction at B. It should be noted, however, that figure 3 
reproduces only a small part of their complete flow diagram. Even with 
the fine computational net which they employed in the vicinity of B, errors 
of the same order of magnitude as the discrepancy noted here are to be 
expected (the differences between their recorded values of (H—y) anc 
q/\ (2g) are of this order). ‘he present method thus appears to be superior 
(in this one respect) to the relaxation method, even when the latter is applied 
with the high degree of precision achieved by Southwell & Vaisey 


4. FLOW UNDER A PLANING SURFACE (HYDROPLANE) 

We now turn to the problem indicated in figure +, where the obstacl 
is a plane inclined at an angle 8(< 90°) to the horizontal. ‘The plane is 
assumed to have been lowered into the stream to an extent sufficient to 
avoid standing waves on the downstream side. ‘This problem has received 
much attention in the past due to its bearing on hydroplanes and sea-planc 
floats. ‘The scope of existing theories was discussed by Southwell & Vaisey 
(loc. cit.), who also treated an example by relaxation methods. ‘They took 
8 = 30°, and assumed the same initial conditions as in their sluice-gaté 


( 
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example. ‘lLheir solution shows that the free surface curves very sharply 
just below the trailing edge B, but beyond this region follows the curve 
obtained in the previous example (figure 3). 

To account for the region near B, we may use the solution of the problem 
indicated in figure 2(c). ‘This was included in the series of jet problems 
solved by Mises (1917); but, as in the example previously mentioned, he 
did not give an expression for the free surface. It can be shown without 
difficulty, however, that the intrinsic equation of the surface is 


2d 6+ 8 
P= ra logsec oY) s F 


ay 


where / denotes the length of the arc. Hence, for 8 = 30, we find in the 


¢ 2sin@ 
ee 1 + | 3 tanh ( =n a sind |, (4.1) 
d T1LVd Var, 


5 


ae 


usual way that 





1-229 for @ 


‘The curves obtained from (2.3) and (4.1) are suitably linked at a point where 
#= —20. Inthe example treated by Southwell & Vaisey, this procedure 
leads to the value 0-777 for the contraction ratio. ‘lhe same value is found 
from the coordinates given in their figure 26. ‘The value in the absence of 
gravity is, from (4.1), 1/1-229 = 0-814. For a refinement of the method, 
the solution of the problem in figure 2(d) may be used instead of (4.1). 


5. EXPERIMENTAL TESTS 

While in no respect being intended as a complete investigation, the 
following experiments serve to illustrate some of the matters considered in 
the preceding sections of this paper. ‘The experiments were carried out 
in a horizontal channel previously used in the work of Binnie, Davies & 
Orkney (1955), whose paper can be consulted for full details of the equipment. 
After the completion of these authors’ experiments, the wooden contraction 
connecting the channel with the supply tank was modified, with the result 
that uniform subcritical flow could be obtained at higher Froude numbers 
than before. ‘The channel was 14in. wide and 8 ft. long, with glass walls 
about 12in. high. When necessary, the stream could be held up by means 
of an adjustable sharp-crested weir at the downstream end; but otherwise 
it fell freely into a sump after leaving the channel. A movable gantry 
supported a point gauge and Pitot tube, which provided accurate measure- 
ments of the stream depth and total head. Froude numbers were calculated 
from the formula F? = 2(total head/depth —1). The water supply was 
controlled in such a way that the flow rate O was independent of any 
obstruction in the channel. ‘This feature assumes significance when the 
theory of $2 is brought to bear on the experimental observations. The 
description of the experiments is suitably divided into three parts, as follows, 
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(i) Experiments in which the obstacle was fixed 


A rectangular brass plate was contructed with a width slightly less than 
that of the channel. It was suspended from a horizontal axle mounted 
above the channel at about 1 ft. from the upstream end, so that it could be 
swung down broadside to the stream. A sharp edge was machined on the 
bottom of the plate, and the small gaps between the sides and the walls of 
the channel were made water-tight with rubber sleeving. Struts attached 
to the rear face of the plate held it in position at any desired inclination. 

In the first experiment, the position of the plate was gradually adjusted 
until the bottom edge was brought into contact with the surface of a uniform 
subcritical stream. A large disturbance was always observed as soon as 
contact was made, provided the Froude number of the undisturbed stream 
was not too small. An appreciable increase in depth occurred on the 
upstream side, and a zone of eddying water was formed immediately in 
front on the obstacle. The latter effect evidently prevented the formation 
of a stagnation point in the manner indicated by figure +. When steady 
conditions were resumed, the obstacle was well immersed in the stream 
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Figure +. Side elevation of planing surface. 


The result on the downstream side varied in character according to the 
initial Froude number and the inclination of the plate. If the Froude 
number was fairly small, say less than 0-3, and the plate sloped backwards 
as in figure 4, a stationary wave-train appeared, its surface being everywhere 
smooth. At rather higher Froude numbers the waves were larger both in 
amplitude and length, and the leading wave broke. At a Froude number 
about 0-7 and with the plate inclined slightly forward, a wave-train appeared 
immediately after contact but was then swept away, leaving the stream 
in a uniform supercritical state. The wave-train could be held in place, 
however, by slightly raising the weir at the channel exit. 
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‘lhe ettects of lowering the plate further into the stream are most easily 
described in reverse, that is, we recount the results of gradually lifting the 
plate backwards from a vertical position as in figure 1. ‘The adjustments 
to the plate were made slowly enough for the flow to approximate the steady 
state appropriate to each stage, and the weir was kept out of action. As 
the plate was raised, the water level fell upstream and rose downstream, 
so that the respective Froude numbers increased and decreased. Eventually 
a Wave appeared near the end of the channel and moved upstream, thus 
causing other waves to form behind it. ‘lhe leading wave then came to 
rest just behind the obstacle, and broke at its crest. ‘This exemplified the 
type of weak undular bore (hydraulic jump) discussed theoretically by 
Benjamin & Lighthill (1954). By comparing the results of several experi- 
ments like this one, the Froude number at which the bore first occurred 
was found to vary with the depth of the stream. ‘The critical value was 
about 1-5 for a depth of about 2-5 in., and increased for shallower streams. 
When the weir was not in use, the bore clearly resulted from the action of 
friction alone By raising the weir it could, of course, be made to form at 
much higher Froude numbers. If the plate was lifted a little further, the 
turbulent zone at the crest of the leading wave became smaller, and in a few 
of the tests it vanished entirely. At this stage, however, some sort of 
instability evidently occurred, for the water level upstream began to rise 
and fall periodically, while the wave-train downstream began to oscillate 
longitudinally. ‘lhe oscillations usually increased in amplitude until the 
surface of the water suddenly fell away from the bottom of the plate, 
whereupon the flow became uniform. 

According to the theory of §2, a uniform supercritical stream cannot 
be formed by the action of an obstacle unless F > 1-25, since a wave-train 
is bound to result from any attempt to reduce F below this value. ‘lhe 

xperiments indicated that, owing to the formation of bores which is not 
considered in the theory, the value of of F at which waves first appear is 
in practice somewhat higher than the theoretical limit. Nevertheless, 
when waves appear for / > 1-25, a significant amount of energy must be 
dissipated in the bore; whereas in the range 1 < F < 1-25,a minute amount 
of dissipation, however small, is sufficient to precipitate a wave-train (see 
Benjamin & Lighthill, p. 452). It seems likely, therefore, that if measures 
were taken to reduce frictional effects and so inhibit the formation of a bore 
(by boundary layer suction, say), the critical value of F could be lowered, 
but could never be reduced below 1-25. 


(11) Experiments in which the wave resistance was specified 


‘To obviate the confusing oscillations observed in the foregoing 
experiments, the following device was employed. As before, a planing 
plate was supported from a horizontal axle broadside to the stream. ‘This 
time, however, the plate was delicately counter-balanced, and a small 
clearance was lett open between each of its sides and the adjacent channel 
wall. ‘lhe whole apparatus was of light construction, and the friction of 
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the support was extremely small; consequently a light finger-touch was 
sufficient to bring the bottom of the plate into contact with the stream. 
The plate could be immersed to any desired extent by loading the plate 
with suitable weights ; hence it was a simple matter to estimate the horizontal 
force exerted against the stream. Disturbances arising from the sides of 
the plate tended to spoil the two-dimensional character of the flow, but 
were small enough to be unimportant. 

By allowing the obstacle freedom of movement in this way, the resistance 
which it experienced on immersion was made independent of small changes 
in the level of the stream. ‘This property was, of course, not possessed by 
the fixed obstacle in the previous experiments, and its absence then was 
evidently responsible for the instability accompanying slight immersions. 
No oscillation was observed with the second device in use. 

A number of experiments were made in which the resistance was increased 
gradually from zero by the addition of weights to the plate. ‘Che waves 
first to appear were very small in amplitude, and their length A was found 
to be in good agreement with the usual formula 


F2 = —— tanh > (5.1) 
*” oe »”? 


where f, and F, are the depth and Froude number upstream. ‘lhe 
amplitude of the waves and the wavelength both became larger as the 
resistance was increased, although the proportional increase in the wave- 
length was quite small. Eventually the leading wave broke. 

The reduction in the quantity S just sufficient to produce breaking waves 
was estimated for various values of Fy. ‘These estimates were found to 
agree fairly well with the theoretical results of De (1955), the method of 
comparison being as follows. If the total drag on the obstacle is M grams 
weight, the reduction in S is obviously S* = Mg/(pb), where 6 is the breadth 
vf the channel. In presenting his numerical results, De considers a dimen- 
sionless quantity s = S/S,, where S, is the value of S for a critical stream 
with the given flow rate Q. It can easily be shown that 

S, = ght FG". (5.2) 
Hence we have, for the reduction in s, 
2M F548 

3ph2b 

After specifying Fy, the required theoretical value of S* can be measured 
on figure 2 in De’s paper. Advantage was taken, however, of an enlarged 
and more detailed version of this figure kindly made available to the author. 
A typical set of measurements, made when the leading wave was on the 
point of breaking, was hy = 14:7 cm, Fy = 0-701 and M = 52 gm. In 
addition we had 6 = 35-6 cm and p = 1 gm/cm’, and with these values (5.3) 
gives s* = 0-0072. The value of s* estimated from De’s chart was 0-0079. 
The slight deficiency of the experimental estimate, which is a common feature 
in all similar tests, may have been due to the well-known fact that in practice 
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stationary waves break before their height reaches the theoretical limit. 
For the conditions quoted here the wavelength of breaking waves was 
estimated to be roughly 50cm. The wave-number 2z/,/A was therefore 
about 1-8, which agrees fairly closely with the value indicated on the chart. 


(iti) Measurements of the contraction ratio for a sluice-gate 

‘These measurements were taken in order to check the theoretical values 
derived in §3. ‘The brass plate used in the first group of experiments was 
fixed in a vertical position, and careful adjustments were made to insure 
that the bottom edge was strictly horizontal. The opening s was 3-571in. 
for one series of tests, and 1-02in. for another. In each series, readings 
of the total head and depth were taken in mid-stream at a distance about 
4-5s downstream from the sluice. ‘The value of H was varied over as 
wide a range as possible by adjusting the flow rate. The maximum useful 
value of H was that at which the water level in front of the sluice rose to the 
top of the channel walls. At the minimum value a bore formed on the 
downstream side. 

Experimental values of the contraction ratio A/s are plotted as a function 
of s/H in figure 5, which includes the theoretical curve according to $3. 














-—— 
0 0-2 0-4 S 6 
Figure 5. Comparison between theoretical values of contraction ratio (——————) 
and experimental values for sluice-openings of 1-02 in. (——-o——) and 
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Both series of tests specified above are represented in the figure. The two 
experimental curves are seen to lie above the theoretical curve, the 
discrepancy being larger for the smaller value of s. ‘The three curves appear 
to converge for small values of s/H, that is, when the downstream Froude 
number is large. At first sight figure 5 tends to exaggerate the influence 
of friction on the flow; but the following interpetation casts a more favourable 
light on the comparison between theory and experiment. 

Following the usual procedure in hydraulic experiments, traces of 
colouring matter were placed in the stream and their movements observed. 
In this way the flow-pattern was seen to be approximately as the theory 
predicts, except in an eddying region near the stagnation point and in the 
boundary layer on the channel bottom—which was apparently always in 
a laminar state in the vicinity of the sluice. ‘The boundary layer on the face 
of the sluice appeared to have only local importance, probably because the 
velocity is small over most of this region. Let us therefore consider that 
the discharge under the sluice is affected by friction only in as far as the 
effective opening is reduced to a value s—6, where 6 is a quantity roughly 
the same as the thickness of the boundary layer. The depth downstream 
is similarly reduced; hence the discrepancy 4 in the contraction ratio is 
approximately proportional to 6/s. If we now tentatively assume that 
s determines the length scale to be associated with the growth of the boundary 
layer, it follows that 4 is proportional to the Reynolds number based on the 
length s. Consequently, for two experiments at the same Froude number, 
we have 4,/4, = (s/s,)°4. Withs, = 1-02 ands, = 3-57, the values in inches 
for the tests in question, this ratio is 2-6. ‘The discrepancy in s/H varies 
in the same way. Examination of figure 5 shows that the experimental 
values are roughly in agreement with this result over most of the range of 
s/H. ‘The extent of the agreement is in fact quite surprising in view of the 
very rough and ready nature of the above calculation. By observing the 
drop in total head as the Pitot tube was brought near the channel bottom, 
the thickness of the boundary layer near the wider sluice (s = 3-57 in.) was 
estimated to be of the order of }in. Note that the corrections necessary 
to bring the experimental results into agreement with the theory are also 
of this order. 

These three groups of experiments amply confirm the theory’s usefulness 
in qualitatively assessing the flow. ‘They suggest moreover that if allowance 
could be made for the boundary layer on the channel bottom, the theory 
would give accurate predictions of the drag on an obstacle and other physical 
quantities. The eddying zone in front of the obstacle was the most noticeable 
effect of friction in the experiments, but it appeared to have little influence 
on the flow downstream. Very complex phenomena, due to surface 
tension and shedding of the boundary layer, were to be expected in the 
region where the free surface first springs clear of the obstacle; but they 
did not appear to affect the overall characteristics of the flow. 

Further experimental work is desirable, especially on the boundary 
layer. For instance, each of a series of measurements as in (iii) might 
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usefully be accompanied by careful estimates of the boundary layer thickness. 
‘This procedure might establish a close agreement with the theory of §2. 
Further, the method used in (ii) could probably be refined a good deal, 
and might profitably be extended to measurements of wave resistance over 
the entire range of De’s theoretical results. 
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